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Developments in the physics of 2D electron systems during the last decade have revealed a new 
class of nonequilibrium phenomena in the presence of a moderately strong magnetic field. The 
hallmark of these phenomena is magnetoresistance oscillations generated by the external forces 
that drive the electron system out of equilibrium. The rich set of dramatic phenomena of this 
kind, discovered in high mobility semiconductor nanostructures, includes, in particular, microwave 
radiation-induced resistance oscillations and zero-resistance states, as well as Hall field-induced 
resistance oscillations and associated zero-differential resistance states. We review the experimen- 
tal manifestations of these phenomena and the unified theoretical framework for describing them 
in terms of a quantum kinetic equation. The survey contains also a thorough discussion of the 
magnctotransport properties of 2D electrons in the linear response regime, as well as an outlook 
on future directions, including related nonequilibrium phenomena in other 2D electron systems. 
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I. INTRODUCTION 



Fundamental research on two-dimensional electron 
gases (2DEGs), as well as on quantum wires and quantum 
dots made on the base of 2D structures, has largely deter- 
mined the developme nt of condensed matter physics in 
the last half-century (iBeenakker and van Houtenl . [l99lt 
iFerrv and Goodnickl . 11997 ). Continuous advances in the 
fabrication of semiconductor structures and the devel- 
opment of experimental techniques and theoretical ap- 
proaches keep bringing up fascinating new areas of re- 
search. On the side of applications, these structures, 
serving as a basis for the planar semiconductor technol- 
ogy, are of enor mous importance for modern nano- and 
optoelectronics (|Alferovl . T200ll ; iKroemeii l200ll ). 

A large body of research findings concerning the 
character of electronic states, mechanisms of scatter- 
ing, activated transport in the insulating regime, the 
cyclotron resonance, quantum magnetooscillations, and 
many other properties of a 2DEG has been summa- 
rized i n the famou s revie w article by Ando, Fowler, and 
Stern ( Ando et all ll982T) . The discovery of the integer 
and fractional quantum Hall (QH) effects in 1980-1982 
largely shifted the focus of research towards the range 
of strong magnetic fie l ds B where the QH physics is 
observed ( von Klitzin 1 119861: iLauehlinl . Il999t IStormerl . 



M Il999t iTsuil . Il999l) . More recently, the behavior of very- 
low-density 2DEGs at zero B has attracted a great deal 
of attention: the physics of these systems is governed 
by the interplay of quantum localization and strong in- 
ter actions, which leads to a metal-insulator transition 
(|Abrahams et aZ.L l200lh . 

During the last two decades evidence has accumulated 
that the range of moderately strong B, where the cy- 
clotron dynamics already sets in whereas the quantum 
localization effects and thus the Hall quantization are 
not developed yet, also reveals a rich variety of deep 
and important physical phenomena. These include ge- 
ometric resonances in pe riodically modulated structures 
( Weiss et all I1989L Il99ll ) and magnetoresistance due to 
quasiclassical memory effects and to interaction-induced 
qua ntum corrections to t he conductivity [for a review, 
see (|Dmitriev et aU [2008h ]. 

A particularly prominent research arena was opened 
by the experimental discovery of t he micr o wave- 
induced resi s tance oscillations (MIRO) (|Ye et al. , 2001; 
IZudov et a/.l . |2001al) and the zero- resis tance states (ZRS) 
( Mani etaU , 120021 : IZudov et all l2003t ). It became clear 
soon after these discoveries that a new field of nonequilib- 
rium physics in high Landau levels (LLs) emerged. Most 
importantly, the observed phenomena demonstrated that 
the combined effect of the weak Landau quantization 
and the relatively weak microwave radiation can give 
rise to very strong changes of the transport properties 
of a 2DEG. These novel nonequilibrium phenomena re- 
vealed unexpected and conceptually interesting physics 
which has proven to be of interest in many areas of con- 
densed matter physics. In addition to MIRO and ZRS, 
a number of other closely related phenomena have been 
observed, such as the Hall field-ind uced resistance oscil- 
lations (HIRO) (|Yang et aLl . l2002a|). the phono n- induce d 
resistance oscillations (PIRO) (fZudov et all l2001bl) . 
in particular, u nder nonequi li brium conditions created 
by a dc field (|Zhang et all l2008bl) , the microwave- 
induce d -B-periodic oscillatio ns in stronger magnetic 



fields (Kukushkin et all 2004). the zero-differential re- 



sistanc e states (ZdRS) (|Bvkov et all 120071: IZhang et al\ , 
2008bh. as we l l as the fractional MIRO (jDorozhkinl 1200 



Will ett et all 12004 IZudov et al. , 2003). 



In the decade that has passed since the above dis- 
coveries, intensive experimental and theoretical work by 
many groups of researchers has greatly advanced the un- 
derstanding of the physics of nonequilibrium phenom- 
ena in high LLs. These advances have motivated us 
to write this review. While several short reviews are 
available (mainly written on a m uch earlier stage of de- 
velopm ent of the research field) dDmitriev et all 1 2008 1: 
i Durstj. [2001 iDurst and Girvinl. |2004 iFitzgerall booa 
iLeil . l2010t iLvapilin and Patrakovj 12004: IVitkalovl . l2009h 
there is a clear need in a full-scale review article on the 
subject. In addition to summarizing the past achieve- 
ments, we also review important open questions as well 
as perspectives for further developments (including re- 
lated nonequilibrium phenomena in other 2D systems). 



3 



Particularly interesting in this respect is the recent dis- 
covery of photoinduced dissipa tionless transport in a 2D 
electron sy stem on liquid He ( Konstantinov and Konol 
2009. 120101) . 



II. LINEAR TRANSPORT OF 2D ELECTRONS IN A 
MODERATELY STRONG MAGNETIC FIELD 

We begin by reviewing the transport properties of a 
2DEG in the linear-response regime, with emphasis on 
the effects that are qualitatively sensitive to the specific 
nature of disorder. 



A. 2D electron gas: Types of disorder 

Almost all work on the novel nonequilibrium magneto- 
transport phenomena at low magnetic field B has relied 
on ultra-high mobility selectively doped GaAs/AlGaAs 
heterostructures — with the electron mobility /i at the 
level of /i ~ 10 7 cm 2 /V s (the mean free path / ~ 10 2 fim). 
Therefore, we focus below on the scattering mechanisms 
that are specific to these realizations of a 2DEG. Note 
that it is on the structures of this very design and this 
sample quality that the most spectacular advances in the 
study of the fractional QH effect have been made in the 
opposite limit of strong B. Typically, the 2DEG in high 
mobility structures is confined to a single quantum well 
(QW) in undoped ("pure") GaAs. In the single- interface 
setup, a triangular-shaped QW is created by the con- 
duction band offset at the interface between AlGaAs and 
GaAs on one side and by the electric field produced by the 
charge of dopants, and of the 2DEG itself, on the other. 
The stability of the 2DEG charge confinement near the 
interface is supported by the chemical potential fixed by 
the overall distribution of charges. The double-interface 
design with a GaAs QW squeezed between AlGaAs layers 
is also frequently employed. 



1. Remote donors 

The key idea behind the suppression of impurity scat- 
tering in the heterostructures is a spatial separation of 
the layer to which the 2DEG is confined and the ionized 
donor impurities by an undoped "spacer" with a typical 
width d ~ 10 2 nm. Since the first experiments on th e 
modulation-doped heterostructures (jDingle et all ll97S|) . 
the optimization of the structure design has led to the 
growth of achieved mobilities from ab out 10 4 cm 2 /Vs to 
the values in excess of 3 x 10 cm /Vs dPfeiffer and Westl . 
120031 : lUmanskv etail [2009). Much of this progress has 
been related to the strong dependence of the momentum 
relaxation rate 1/r = e/m/j, on d for scattering off remote 
charged impurities (here and below h = 1): 



(1) 



where m is the electron mass, kp is the Fermi wavevector 
of the 2DEG, m — Jdz cn(z)(d/\z\) 3 , and c R (z) is the 
three-dimensional density of remote charged impurities 
located at \z\ > d (z is counted from the interface). 

The transport scattering rate at zero temperature T is 
obtained as an integral over the scattering angle </>: 



in 



2n 



(1 - COS (f>)W q=2kF sin^/2) , (2) 



where W q is the Fourier component of the correlation 
function of the screened impurity potential in the plane 
of the 2DEG 



W{r) = (V(Q)V(r)) , r = (x,y) 



(3) 



at the transferred in-plane momentum q. In the limit 
kpd 3> 1 (typically kpd ~ 10), scattering is predomi- 
nantly on small angles and Eq. ((5J) reduces to 



1,3 
F 



J dq 
2^ 



q 2 W q 



(4) 



The function W q that describes remote charged impuri- 
ties and upon substitution in Eq. dU yields Eq. (JTJ) is 



W q =[- 



,-2qd 



(5) 



Equation (|SJ) assumes that (i) q -1 is much larger than 
the 2DEG thickness w (typically w ~ 10 nm); (ii) the 
only source of screening of the impurity potential is the 
2DEG and the dielectric function e(q) is given by the 
2D random phase approximation (RPA): e(q) — en(l + 
2/qag), where eo is the lattice dielectric constant at the 
interface and = eo/me 2 is the Bohr radius; (iii) q^ 1 
is much larger than as — 10 nm. Conditions (i) and (iii) 
for the characteristic values of q^ 1 ~ d are reasonably 
well satisfied in the high-mobility structures. The use of 
the RPA in (ii) relies on the smallness of the Wigner- 
Seitz parameter r s — [i:n e a 2 B )^ 1 /' 2 , where n e — k 2 F /2ir is 
the density of the 2DEG. For the typical value of n e = 
3 x 10 11 cm -2 , r s ~ 1 and the conventional use of the 
RPA is only marginally justified. Fluctuations of the 
random potential with the correlation function (|SJ) can 
be accurately described by Gaussian statistics in view of 
the typically large parameter riid 2 ^> 1. 
Note that the total scattering rate 



2n 



W, 



q=2k F sinO/2) 



(6) 



[which determines, e.g., the LL broadening induced by 
disorder, see Sec. III. C] is, in the case of remote impurities, 
much larger than 1/r. In the limit kpd 3> 1, l/^q is 
rewritten, similarly to Eq. ([!]). as 



1 _ 2m r°°dq 
T q k F Jo 2tt q ' 



(7) 



8m(k F d) 3 



which for W q from Eq. ((SJ) gives r/r q = (2kpd) 2 
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In an important case of "delta-doping" , charged impu- 
rities are concentrated in a thin layer whose thickness is 
much smaller than d (down to a few lattice constants), 
so that cr(z) can be approximated as riibiz + d) with 
rii being then the sheet density of these impurities. The 
(5-layer is grown either by directly adding impurities in 
AlGaAs or by doping into a short- period AlAs/GaAs su- 
perlattice embedded i n the alloy (jFriedland et all Il996t 
lUmanskv et al. , 2009). In the latter case, X-valley elec- 
trons in AlAs yield additional screening of disorder and 
a further reduction of 1/r. In modern high-mobility 
GaAs/ AlGaAs structures, the growth sequence design 
typically includes one or two Si-doped (5-layers (these 
may be followed by other layers with Si dopants which 
compensate for surface states at the boundary with the 
vacuum) , which are the main supply of electrons to the 
2DEG and are thought to be the main source of scatter- 
ing in the 2DEG as far as the intentionally doped layers 
are concerned. A setup with two <5-doped layers, one on 
each side of the 2DEG, has the advantage that it yields 
larger n e for given d. The effect of the increase of n e on 
1/r is stronger than the enhancement of the strength of 
disorder. In the case of remote impurities with n e = n,, 
the time r for a QW symmetrically doped from both sides 
is, according to Eq. JT]), a factor of larger than for the 
single-side doping. Moreover, the double-side doping has 
been instrumental in attaining the highest reported mo- 
bilities for which the scattering off remote impurities is 
likely to be of little importance, see Sec. III.A.2I 

Equation (|TJ) implies that charged impurities are 
randomly distributed in the 2D plane. This as- 
sumpt io n appears to overestimate 1 / r dBuks et all . 



2. Background impurities 



Electrostatics of the heterostructure ( Sternl. 19831) dic- 
tates that n e is lowered with increasing d ( Heiblum et all 
Il984t lUmanskv et all Il997t ). Details of the relation be- 
tween n e and d depend on the concrete doping design of 
the heterostructure; however, as long as the electric field 
in the spacer layer is mainly determined by the 2DEG 
charge, n e scales as 1/d, which gives for 
cording to Eq. ©, oc d 5 / 2 (jLee et all Il983l) . The 



growth of /i with increasing spacer thickness, observed 
at smaller d, typically sto ps in the ultra-high m obility 
samples at d ~ 60-70 nm (lUmanskv et all Il997l) . This 
behavior is usually argued to be associated with a com- 
petition between scattering off remote charged impuri- 
ties [Eq. ([]])] and scattering off "background impurities" 
which are present in a small concentration in the spacer 
layer and also in the GaAs layer. Being distributed in 
close vicinity of the 2DEG, the background impurities 
can lead to large-angle scattering, thus giving the main 
contribution to 1/r at sufficiently large d. 

The momentum relaxation rate due to scattering off 
background impurities uniformly distributed with the 
three-dimensional densities c B on one side of the interface 



7TC B 



mk F a B 



— In (min{fcp, w 1 }a B ) 



(8) 



where c B 



and the corresponding correlation 



function W q for qui <C 1 is given by 



W q = (c B /2q)[27re 2 /e(q)q} 



19 

(C 



Coleridgi 
tures with d 



19911 11997ft measured in some heterostruc- 
The reason why t he measured 



< 50 nm. 



mobility is higher than expected was argued (jBuks et all 



)ility 

ll994aT lbf) to be related to the fact that the Si-dopants in 
AlGaAs (in a range of Al content which includes the typi- 
cal values of 30-40%) may exist in two configurations: as 
shallow donors (which become positively charged after 
having supplied electrons for the 2DEG) and as DX cen- 
ters (whose energies lie deeper in t he forbidde n gap a nd 
which can be negatively charged). iBuks et a/1 ( 1994al lbh 
interpreted their experimental results to prove that the 
DX centers are negatively charged in the ground state 
and demonstrated that correlations in the spatial distri- 
bution of the positively and negatively charged impuri- 
ties may lead to a strong reduction of the sc attering rate , 
simil ar to the correlations considered earlier ( Efros et all 
Il990h for the case when not all donor impurities are ion- 
ized and some remain neutral. Control of the ratio of 
the impurity densities in the shallow and D X configura- 
tions i nd eed provided exp erimental evidence (jBuks et all 
Il994allbl : [Colerideel . Il997ft for the correlations capable of 
substantially reducing 1/r (by a factor of up to about 6), 
depending on the parameters of the technological process, 
and a similar beha vior was also demonstrated in the mea - 
surements of l/r q (|Coleridgel Il997t IShikler et ail . Il997l) . 



(9) 

1994alfbl: IColeridgt^ 1997t iHeiblum et aLl . [l984 ) and l/r q Equation ((Sj> is accurate provided both kp 1 and w are 

much smaller than a B . In reality, both spatial scales in 
high-mobility GaAs/ AlGaAs structures are about a B — 
10 nm, so that 1/r due to scattering off background im- 
purities should then be sensitive to the exact shape of the 
electron density profile in the direction across the 2DEG 
plane. A comparison of Eqs. (TTJ) and ([5]) shows that at 
rii/d 3 ~ c B /a B (up to the logarithmic factor) the two 
mechanisms of scattering give equal contributions to the 
total scattering rate and at larger d the scattering on 
background impurities dominates. 

The dependence of fi on n e varied in situ by a con- 
trolled illumination in structures with d larger than 85 
and up to 200 nm, exhibited — when fitted to a power law 
over about one decade in n e within the interval between 

10 10 cm -2 and 2 x 10 11 cm -2 — the relation a cx rii with x 
i 1 1 — — — 1 1— 1 1 

in the range 0.6-0.7 (Pfciffcr et aLI.119 89: Shavcgan et all 
119881 : lUmanskv et all Il997fh iPfeiffer and Westl (|2003f ) 
also emphasized that the power law with x ~ 0.7 holds 
over two decades in n e for a series of samples with dif- 
ferent n e for /i > 10 6 cm 2 /Vs. These values of x are 
substantially smaller than 3/2 expected for scattering 
off remote impurities [Eq. ([T])]. which is commonly re- 
garded as evidence pointing toward the dominant role 
of background impurities. Note that Eq. ((5]) for back- 
ground impurities also predicts the exponent 3/2 for 
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sufficiently large n e . However, in the low-density case 
n e < 10 n cm~ 2 , the dependen ce of a on n P — accordin; 
to the numerical calculations ( Andd . Il982t Goldl . Il98 
using variational electron wavefunctions to treat the ef- 
fect of a finite thickness of the 2DEG — is indeed substan- 

3/2 

tially weaker for backgrou nd impurities than /i oc n e 
[the results bv lAndol ( 1982h for background impurities are 
obtainable from those for remote impurities by putting 
d — > 0] . On the other hand, the nume rical variational cal- 
culation bv lWalukiewicz et all (1984) for n e > 10 11 cm~ 2 
gave very similar dependences of /i on n e for remote and 
background impurities. 

A systematic experimental study ( Umanskv et all 
Il997l) of the relative contributions of the scattering mech- 
anisms was performed in single-interface heterostructures 
grown in a wide range of d and n, and mobilities as 
high as [i ~ 1.4 x 10 7 cm 2 /Vs. By measuring the mo- 
bility at given d and n e as a function of rii (see also 
lUmanskv et all . 12009^ . experimental evidence was pre- 
sented that for d larger than 70 nm the contribution of 
remote impurities to the total 1/r was smaller (in the 
"best samples" — about 10-15%) than that of background 
charges. These results indicate that /i in the presently 
available ultra-high mobility structures of various designs 
is likely to be mainly, or to a large extent, limited by 
background impurities. 

In both Eqs. ([T]) and (JSJ), 1/r was obtained within the 
Born approximation. While this level of approximation 
is well justified for remote impurities, for kpciB ~ 1 back- 
ground impurities that are within the distance as from 
the 2DEG and give the main contribution (as far as the 
background charges are concerned) to the scattering rate 
are not really weak scatterers. The momentum relax- 
ation time r for the background impurities sitting effec- 
tively right at the interface, is given simply by the time 
of flight l/vpCBa% (where vf is the Fermi velocity) be- 
tween two of them along the straight line. The presence 
of the strong scatterers at the interface will have impor- 
tant ramifications in the subsequent sections. 



3. Surface roughness 

Interface roughness is also an important source 
of disorder in high-mobility GaAs/AlGaAs structures 
and is argued (IMarkus et all 11994 ISaku et all Il996t 
lUmanskv et a/.I . Tl997 ) to be probed by the anisotropy 
displayed by [i when the current is measured along 
different directions in the 2D plane. The anisotropy 
of interface roughness is viewed as an inherent prop- 
erty of the kinetics of growth of GaAs-based struc- 
tures. Grown, as is usual in high-mobility structures, 
on the (001) plane, they exhibit atomic-scale terracing 
at the interface, with randomly shaped terraces being 
typically elongated in the [110] direction. Anisotropy 
of fi was reported to be in some h igh-mobility sam- 
ples as large as 40% (IMarkus et all. 1994 ;_ Saku et all , 
ll996tlTokura et aUll992t lUmanskv et aUll997|) . increas- 



ing with electron density. The experimental results in- 
dicate that the strength of interface roughness varies de- 
pending on the g rowth conditions: f or two ultra-h i gh mo - 
bility samples in lSaku et all (jl996l) ; lTokura et all (|l992[ ). 
both having almost the same [i ~ 10 7 cm 2 /Vs along the 
[110] direction, the anisotropy was about 40% in one of 
them and about 10% in the othe r . In the "best struc- 
tures" studied by lUmanskv et all ( 19971) the anisotropy 
was also reported to be reduced to 5-10%. The exper- 
iments suggest that the characteristic correlation radius 
of the interface in homogeneities 1 that are relevant to the 
mobility is 10 nm (|Markus et al. L 11994 ISaku et al. 1. 119961: 
iTokura et q"Zl . ll992l) . 

Modeling the interface as a hard wall, the problem of 
scattering by surface roughness reduces to that of a per- 
turbation imposed on the boundary condition. Spatial 
fluctuations of the height of the surface h(r) along the z 
axis create then a 2D random potential whos e correlation 
function at the transferred m omentum q is ()Ando et all 
119821: iPranee and Ned . fl968h 



(10) 



where F — J ' dz (dU j 'dz)p^(z) — (dx/dz) 2 \ z=+0 /2m, 
U(z) is the unperturbed confining potential, pd(z), 
nonzero for z > 0, is the unperturbed electron density 
profile (normalized to unity) in the direction normal to 

1/2 

the 2DEG plane, xi z ) — Pd ( z )' anc ^ (^)q ^ s tne cor " 
relation function of the surface corrugations. Equation 
(fTTj| is valid for the characteristic amplitude of h(r) much 
smaller t han th e 2DEG thickness. Estima t es based on 
Eq. (1TU1) (lAndol. Il982t iMarkus et all Il994t ISaku et all 
Il996t ITokura et all 11992) show that the experimentally 
observed strong anisotropy can be explained in terms of 
scattering by interface roughness for realistic parameters 
of the latter. If the charge of the depletion layer can be 
neglected compared to the 2DEG charge, then F oc n 
and the dependence of the mobility on n e , as it follows 
from Eq. (fT0|). is /i oc n~ 2 , i.e., with increasing n e this 
mechanism of scattering becomes stronger — in contrast 
to scattering off cha rged impurities. The calculation by 
IMarkus et all (|l994l) shows that the n~ 2 scaling of fx (for 

— 1/2 

F oc n e ) changes to n e if scattering by interface rough- 
ness is predominately on small angles. 

Although the term "interface roughness" commonly 
refers to the long-range fluctuations of h(r), there are 
also inherent imperfections in the interface related to al- 
loy disorder in the AlGaAs part of the heterostructure. 



1 The issue of anisotropic interface roughness in ultra-high mobil- 
ity GaAs/ AlGaAs structures has attracted attention also in the 
context of co operative effects in electron transpo r t at h alf-filling 
of high LLs l lCooper et all |2001| ; IWillett et all l200lf) . These 
experiments probed surface roughness at the interface with the 
vacuum and were focussed on long-range statistically anisotropic 
fluctuations of the surface height with a much larger correlation 
radius on the scale of 1 fira. 
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These are correlated on the scale of the lattice constant 
ag and are described for qa <C 1 by 



W q = A 2 c all, 



(11) 



where A c is the energy scale of the fluctuations and the 
integral / = fdz p 2 (z), with p(z) being the electron den- 
sity profile (normalized to unity) across the 2DEG, is 
taken over t he region where alloy scattering is present 
(|Andol . Il982t) . The integral is estimated as / ~ £, 5 /w 6 , 
where £ ~ ao is the tunneling length under the barrier on 
the AlGaAs side of the interface and w is the total thick- 
ness of the 2DEG. Although the values of A e and I can be 
calculated only numerically, early estimates ( Andol . Il982l) 
show that for n e = 3 x 10 11 cm -2 and the Al fraction in 
the range 30-40% alloy scattering at the interface estab- 
lishes the mobility limit at the level of about 10 7 cm 2 /V s. 
This suggests that alloy scattering, commonly deemed ir- 
relevant at the GaAs/ AlGaAs interface, may in fact be 
important in the modern ultra-high mobility heterostruc- 
tures. 



B. Magnetotransport beyond the Drude theory 

At the most basic level, dc transport in a magnetic 
field B is described by the Drude formula according to 
which the resistivity tensor p(B) in two dimensions reads 



P(B) 



1/t W c 
— LU C 1/t 



(12) 



where uj c — eB/mc is the cyclotron frequency. Equa- 
tion (IT2|) can also be derived within the formalism of the 
Boltzmann kinetic equation which has been very success- 
ful i n explaining the magnetotransport p roperties of met- 
als (jLifshits et ad Il973l iPippardl . Il989l) . the behavior of 
p(B) in which is known to depend in an essential way 
on the shape and topology of the Fermi surface. For an 
isotropic 2DEG, which is a virtually perfect approxima- 
tion for the 2DEG in GaAs as far as the band structure 
is concerned (note, however, the possible importance of 
anisotropic disorder discussed at the end of Sec. III. A[) . 
this theory predicts the simple Drude result (fl"2"|) . A strik- 
ing feature of Eq. (|12l) is that the longitudinal resistivity 
Pxx(B) — pr> = m/e 2 n e T is independent of B, so that 
the deviation 



Ap xx = Pxx(B) - p xx (0) , 



(13) 



termed a positive/negative magnetoresistance (MR), de- 
pending on the sign of Ap xx , is exactly zero. 



1. Quantum magnetoresistance 

The MR in 2D systems has been studied extensively — 
in the last three decades particularly in regard to In- 
dependent contributions coming from the influence of 
a magnetic field on the quantum corrections, small in 



the parameter 1/kpl < 1, to the transport coefficients 
( Altshuler and Aronovl Il985l ). One quantum contribu- 
tion to the MR is related to the suppression of the 
weak- localization correction to pr>, which occurs in a 
"classically weak" (in the sense of uj c t <C 1) magnetic 
field and results in a small sharp spike of p xx {B) at 
B = 0. Another quantum contribution is due to the 
interaction between degenerate electrons multiply scat- 
tered by disorder. In the diffusive regime (Tt <C 1), 
the Coulomb interaction-induced MR obeys (for r s «l 
and uj c t > 1) Ap xx /pu = — (uj 2 t 2 /irkpl) lnfl/IY), irre- 
spect i ye of the nature of d i sorder (lAltshuler and Aronovl 
119851 : iGirvin et all Il982t iHouehton et all Il982ri . The 
theory of this MR was exten ded to arbitrary Tt by 
iGornvi and Mirlinl (|2003l . l2004h . In the ballistic regime 
(Tt 3> 1), particularly relevant to high-mobility struc- 
tures, the MR caused by electron-electron interactions 
depends in an essential way on the type of disorder, es- 
pecially for uj c <C T, where it is strongly suppressed 
for smooth disorder. For lj c 3> T, if 1/t is limited by 
smooth disorder, the Coulomb interaction-indu ced MR 
for r s < 1 reads (|Gornvi and Mirlinl . l2003l 120041 ) 



Ap xx = 3C(3/2) 1 {lu c t) 2 
pu 32tt 3 / 2 k F l (Tt) 1 / 2 



(14) 



for Tt <C r~ 2 , while for Tt ^> r~ 2 it retains the same 
amplitude as in Eq. (fl"4"|) but changes sign. If, by contrast, 
1/r is limited by white- noise disorder but the smooth 
component of disorder is still strong enough to produce 
a contribution to the momentum relaxation rate 1 /tl 3> 
uic/kpVp, the MR is enhanced compared to Eq. (fT4"]) by 



a factor 4(t l /t) 1 / 2 (|Gornvi and Mirlinl . [200 



A crossover in the temperature behavior of the 
quadratic-in-S MR between the diffusive and ballis- 
tic regimes was observed in the intermediate range of 
Tt ~ 1 in modera te-mobility GaAs /AlGaA s structures 
bv iLi et all (l2003f) [se e also earlier works (|Choi et ali 
11986 ; Paalanen et al. |, I1983T) ]. in close agreement with 
the theory for smooth disorder. The behavior of 
the MR in the cross over region was al s o stu died in 
Si/SiGe structures by lOlshanetskv et al. | (|2006l) . The 
suppression of th e MR as T increases, observed by 
iGalaktionov et all ( 20061 ) in samples similar to those in 
( Li et al. . l2003h in a wider range of Tt > 1, agrees 
favorably with the T~ x / 2 scaling characteristic of the 
ballistic regime. The enhancement of the MR com- 
pared to Eq. (fT4|) . which can be largely attributed 
to the interplay of the short- and long-range compo- 
nents o f disorder, was reported in a number of exper- 
iments dBockhorn et al. . 1201 ll : iGalakt lonov et all 120061: 
lOlshanetskv et alT 2006), with a particularly s trong en- 
hancement observed by Bockhorn et al\ ( 201ll ) in ultra- 
high mobility GaAs/ AlGaAs structures. 
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2. Classical magnetoresistance 

Quite apart from the quantum MR, it was half a cen- 
tury ago that it was realized that long-range inhomo- 
geneities may have a profound effect on magnetotrans- 
port; specifically, that even weak inhomogeneities with a 
spatial scale l arger than the mean free path I can yield 
a strong MR ( Dreizin and Dykhnel . Il973t iHerringl . Il96dt 
llsichenkol . Il992j) . In fact, even in macroscopically ho- 
mogeneous (on the scale of I) electron systems, disor- 
der with the correlation radius d S> kp 1 can produce a 
strong quasiclassical MR (QCMR). The QCMR was dis- 
cussed in a variety of systems: in the 3D Coulomb plasma 
dMurzinl. fl98l iPolvakovi |l986h. in the 2D Lorentz-gas 
model (|Baskin et all \197&L iBobvlev et all [l995h . and in 
the 2DEG with smooth (iMirlin et aUll999h and "mixed" 
( Mirlin et al .1 120011 : iPolvakov et all boOlT disorder. The 
essence of this phenomenon is that transport retains sig- 
natures of the underlying quasiclassical dynamics of elec- 
trons, which are not captured by the Boltzmann-Drude 
kinetic theory. Specifically, the QCMR is due to correla- 
tions in the otherwise random multiple scattering process 
at the points where the quasiclassical paths self-intersect. 
The strength of these "memory effects" — neglected in the 
collision-integral formalism — grows as a power of d/l. Be- 
low, we briefly discuss the results for the QCMR in a 
2DEG. 

The QCMR depends in an essential manner (even its 
sign does) on the correlation properties of disorder. Let 
us first review the case of Gaussian (in the sense of fluc- 
tuation statistics) disorder with the co rrelation radius 



d > fcp 1 . The M R is in fact nonzero dKhveshchenkol 



119961: IMirlin et aLl . ll998h within the collision-integral for- 
malism because of the cyclotron bending of electron 
paths on the scale of d. This gives a small nega- 
tive MR, Ap xx /pD ~ -(d/R c ) 2 , where R c = v f /lj c 
is the cyclotron radius. Correlations of diffusive paths 
at the points of self-intersection give rise to a much 
stronger posi t ive M R: for W q from Eq. (|5|) and lj c t S> 1 
(IMirlin et aLl . ll999h 



Ap xx = 2C(3/2) (d 
Pd n \l 



(w c r) 



9/2 



(15) 



The characteristic momentum transfer q in the scattering 
processes leading to Eq. (fT5j) is q ~ 1/ A <C l/d, where 



(16) 



is the mean-square shift of the guiding center of a cy- 
clotron orbit after one revolution. The ratio Ap xx /pn in 
Eq. (p~5|) is much larger than that related to the effect of 
B on the collision integral for uj c t ^ (l/d) 2 / 5 (note that 
l/d ~ 10 3 in ultra- high mobility structures) and becomes 
of order unity at lo c t ~ (l/d) 2 / 3 . The latter condition 
corresponds to A ~ d. As B is increased further, the 
strong positive MR is foll owed by an exponential falloff 
of p xx (jFogler et all . Il997t) : 



which is due to the increasing adiabaticity of the electron 
dynamics and the related classical localization. The MR 
in Eq. (|T5j) . induced by rare (A ^> d) self-intersections of 
diffusive trajectories, may be considered as a precursor 
of the strong (adiabatic) localization. Numerical simula- 
tions (jMirlin et all Il999h confirm the strong QCMR for 
the case of smooth disorder. The QCMR of composite 
fermions scattered by effec tive random magnetic field was 
argued (jEvers et all Il999h to explain th e dependenc e of 
p xx (B) around a half-filling of LLs (e.g.. ISmetl . Il998l ). 

Rare strong scatterers randomly distributed on top 
of smooth weak potential fluctuations can profoundly 
change the magnetotransport propertie s of a 2DEG 
(|Mirlin et aill200ll [PoTvakov et all l200lh . As discussed 
in Sec. III. At the combination of these two types of 
disorder — of comparable strength in what concerns the 
momentum relaxation rate — is likely to be an adequate 
description of a random potential in high mobility struc- 
tures. The Drude formula (fT2")) totally fails at B ^ 
when only strong short-ra nge scatterers are p resent: in 
an important early work bv lBaskin et al. (1978). a strong 
negative MR was shown in the classical 2D Lorentz-gas 
model. In this model, electrons are scattered by impen- 
etrable hard disks ("voids"). In the limit of the den- 
sity of the voids n$ — > oo with the momentum relax- 
ation time ts held fi xed, the model is exactly solvable 
( Bobvlev et all [T995T) for the resistivity tensor p(B)- in 
particular, p xx /pu — 97t/8w c ts for oj c ts 3> 1 and the 
MR is exponentially small in the opposite limit. At finite 
ns, the model shows a classical metal-insulator transition 
at a critical value o f R c ~ ^g 1 ^ 2 ( Baskin et all Il978t 
IBobvlev et all 1 19951) and a quadratic MR in the low-B 
limit ([Cheianov et all [2004) . 

When both types of disorder are present, the momen- 
tum relaxation rate at zero £? is 1/t ~ 1/tl+1/ts, where 
1/ t l,S stands for the contributions of smooth disorder 
(L) and strong scatterers (S). Consider the relevant case 
of ts <C tl but smooth disorder being strong enough to 
produce A [Eq. (|T6)) with r —> tl] larger than both cor- 
relation radii of disorder d and a, where a is the effect ive 
radius of strong scatterers. Then (| Mirlin et HI. 12001b 



Ap xx /p u = -(uj c /u) ) 2 



(18) 



(17) 



for w c < wo, where 

lo = (2Trn s ) 1/2 v F (2T 2 /T L T 0S ) 1 / 4 , (19) 

with tqs ~ ts being the total scattering time for strong 
scatterers. The ratio (lj c /ujo) 2 gives the fraction of the 
area "explored" twice because of the self-intersections of 
electron paths, which reduces the exploration rate and 
leads to a longer time between collisions with different 
strong scatterers, hence the negative sign of the MR. This 
should be contrasted with the positive MR (fT5|) for one- 
scale smooth disorder, where the passages through the 
same area increase the scattering rate. For lj c 3> ljq, 
the scattering rate is strongly suppressed by the memory 
effects, which gives the B~ 4 falloff as B increases for 



8 









11 


\ v ^^^^^ C 






\ V b 


T=0.3 K 


• 


\. a 


1 — \ — : — : , i , , 



0.0 



0.5 1.0 
5(kG) 



1.5 



FIG. 1 Magnetic-field dependence of the resistivity p xx (in 
units of p X x at B — 0) of three different samples with the 
mobility of (a) 3 x 10 7 cm 2 /Vs, (b) 1.2 x 10 7 cm 2 /Vs, and 
(c) 8.6 x 10 6 cm 2 /Vs at T = 0.3 K. Adapted from IDai eiail 
\20ld ). 



n s R 2 c > 1, (jMirlin et all . l200lt IPoTvakov et al. 1. 1200111 

pxx/pD ~ (T S /T L )(n s R 2 c ) 2 , (20) 

and a plateau of p xx for nsR 2 <C 1 with 

Pxx/pu = t s /t l , (21) 

i.e., the magnetic field "switches off' short-range scatter- 
ers which give the main contribution to 1/t at B = 0. 
This behavior of the QCMR has bee n observed in numer- 
ical simulations ( Mirlin et all l200lf ) . 

A negative MR with a pronounced resistance minimum 
at about 0.5 kG was observed in early works on ultra-hig h 
mobility structures ( Smetl . [T"997t lUmanskv et al. . 1997). 
Recently, a strong negative MR in the same range of B 
was reported in s imilar structures by IDai et al. (|2010h 
and lHatke et all ( 2011bl). In partic ular, in one of the 
samples studied by IDai et all (2010), namely sample A 
in Fig. [TJ the resistivity was shown to decrease by a fac- 
tor of about 50 between B = and 0.5 kG and exhibit 
a wide plateau for larger B before the onset of magne- 
tooscillations, in qualitative agreement with the picture 
of the QCMR for the model of two-component disorder. 
With decreasing mobility (samples B and C in Fig.[TJ, the 
plateau disappears but the negative MR persi sts. Anal- 
ysis of the measured dependence of p X x{B) led lDai et all 
(2010) to conclude that the observed phenomenon can be 
co nsistently described in terms of the QCMR theorized 
bv lMirlin etUll (|200ll ); iPolvakov et all (|200lh . 

The behavior of p xx (B) in the model of two-component 
disorder proves to be substantially more intricate when, 
as B increases, A becomes smaller than 2d (in high- 
mobility structures d^> a) and scattering by the smooth 
rando m potential acquires t he character of an adiabatic 
drift (iPolvakov et all . l200l[ ). Although the adiabatic- 
ity of scattering suppresses p xx when only smooth dis- 
order is present [Eq. (|TT|) ]. it can also give rise to the 
growth of p xx with increasing B in the presence of short- 
range disorder. This behavior can be most clearly seen in 
the "hydrodynamic" model where short-range disorder is 



characterized by white-noise correlations and the prob- 
lem is mapped on to that of advection-diffusion transport 
( Isichenkol . Il992f) . In the limit of large B, the hydrody- 
namic model predicts p xx /pr> ~ (TjvFd/TLi? 2 ) 5 / 13 . Note 
that for realistic parameters the full-fledged QCMR in 
the adiabatic regime competes with the effects of Landau 
quantization (Sec. IlLCj) . A g rowth of p x x with increas- 
ing B in the experimen ts by IDai et all (|2010f ) (Fig. [T]) 
and lHatke et all (|2011bh was observed in a range of B 
well befor e the onset of oscillations [see also the earlier re- 
sults (e.g.. lGalaktionov et aZ.l . l2006h on a positive MR ob- 
served at sufficiently high T in moderate-mobi lity struc- 
tures]. T he huge negative M R observed by IDai et all 
(|2010h and lHatke efcdl (|2011bJ) . as well as the growth of 
p xx at higher B, warrants further study. In this respect, 
it is interesting to note that, experimentally, there ap- 
pears to be a connection between the huge negative MR 
and a pronounced anomaly in the photoresponse, for de- 
tails see Sec. IVII.A.2I 



3. Classical magnetooscillations of the ac conductivity 

Within the Boltzmann kinetic theory, the dissipative 
diagonal ac conductivity cr(uj) = <j+{uj) + er_(w), which 
determines the absorption rate for linearly polarized elec- 
tromagnetic waves, is given by the Drude formula 



1 



2m 1 + (uj c ± w) 2 t 2 



(22) 



The quasiclassical correlations in the dynamics of mul- 
tiple scattering, which lead to the MR discussed in 
Sec. III. B. 21 can also strongly modify the ac response of 
a 2DEG at B ^ 0. In particular, they yield periodic-in- 
lu/l) c oscillations in cr(u>) which are of essentially clas- 
sical nature. For the model of two-component dis- 
order (Sec. III.F3.2[) . the oscillatory classical correction 
Aa^(u>) — (u>) — ar,(uj) in units of the Drude con- 
ductivity <7d(w) takes, in the large- B limit, the form of 
sharp resonant features at the harmonics of the cyclotron 
resonance (CR) at ui — M ui r . with \M \ = 1,2,..., the am- 



plitude of which is given (|Dmitriev et all |2004) by 

Aa^(Miu c ) = a^TL) 1 / 2 
<j u {Mlu c ) 



37r|Af|A 



(23) 



where a — 1/nsvpTs, and the width Tm — 3A/ 2 /2tx + 
1/tqs- Equation (|23[) is valid for Tm *C lo c . In the oppo- 
site limit of strongly d amped oscillations, Aer^ (w) reads 
()Dmitriev et all . |2004 



a 2nuj 

■—rrn cos — ex P 



3nuj 2 

UJ 3 T L 



(24) 



The behavior of (uj) is illustrated in Fig. [21 where also 
the quantum oscillations related to Landau quantization 
(Sec. III.C.3[) . are shown. One sees that the classical os- 
cillations can be much more pronounced than the quan- 
tum ones for r q <ti tx (which is the case in high-mobility 
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FIG. 2 Classical (c) and quantum (q) oscillatory ac conduc- 
tivity [normalized to the Drude conductivity (Td(w)] vs ui c /u 
for u/2-k = 100 GHz, r h = 0.6ns, T L /T a = 50, t s /tl = 0.1 , 
and a/A = 0.25 at oj c /lj = 1/2. From iDmitriev et al\ (|2004T ). 



structures). Resonant features, of similar origin, at the 
harmonics of the CR we re also discussed for t he case of 
a random antidot array (|Polvakov et aUl2002t) . The CR 
in the opposit e limit of a smooth ra n dom p otential was 
considered by iFogler and Shklovskiil ( 1998() who found, 
in particular, a sharp jump, as B is increased, in the CR 
broadening from the Drude width 1/r to the width of 
separated LLs {lOc/t) 1 / 2 at the crossover to the adiabatic 
localization regime. 



C. Quantum magnetooscillations 

In our discussion so far, we have neglected the effects 
related to Landau quantization by assuming that the 
electron DOS v{e) at e equal to the Fermi energy ef 
does not depend on B. In this section, we review the os- 
cillatory behavior of the transport coefficients of a 2DEG 
with varying £p/oj c and uj/lo c — which arises from Lan- 
dau quantization and is collectively known as "quantum 
magnetoos c illatio ns" — in the linear- response regime [see 
IShoenberd (|l984ft for a review of the basic concepts and 
of the resul t s for magnetooscillations in 3D metals, and 
lAndo et all (|l982h for an early review of magnetooscilla- 
tions in a 2DEG]. 



1. Magnetooscillations of the density of states 

In a clean noninteracting 2DEG the DOS per spin is a 
sum of delta-functions at energies equal to those of LLs: 



ie) = (2nX 2 )- 1 J2S(e-e N ) 



(25) 



N 



where A = (mw c ) -1 / 2 is the magnetic length and En = 
(N + 1/2)lu c . As throughout the review, we concen- 
trate on the limit of moderately strong magnetic fields in 



which the number of occupied LLs is large. For N ^> 1, 
the broadening of the Nth peak in Eq. (|23|) by disorder 
changes only weakly when N changes by unity. The DOS 
of the Nth LL i/W(e) is then written as 



Im ■ 



2tt 2 A 2 e - e n - E(e) 



(26) 



where S(e) is the self-energy of the disorder-averaged re- 
tarded Green's function in the LL representation. The 
approximation is that S(e) is treated as independent of 
N. The singular character of Eq. (|25f) implies that 
cannot be calculated at the level of the Born approxi- 
mation when LLs are separated, i.e., when disorder is so 
weak that — Im£(£jv) < uj c . 

The mathematically simplest scheme to calculate 
E(g) is the self-con sis tent B o rn approximation (SCBA 

(lAndo and Uemura . 1974 ILaikhtman and Altshule: . 

ll994tlRaikh and Shahbazvani fl993l ) in which T,(e) obeys 



E(e) 



1 



27TT n 4^ 



A' 



e N - E(e) 



(27) 



where l/r q is the total (quantum) scattering rate [Eq. ©] 
at B = 0. The SCBA gives in the limit w c r q 3> 1, where 
LLs are well separated: 



/WW = 



/ 2uj r 



(e) = v r q J — - - (e- e N f 

V nT i 



(28) 



if the argument of the square root is positive and zero 
otherwise , with z/n = m/2it being the DOS per spin at 
B = 0. iBenedict and Chalked (jl986l) and ICarra et al\ 



(|l989h showed that the SCBA— and thus the semicir- 
cle shape of the broadened LLs in Eq. (|28p — are exact 
for white-noise disorder in the limit N —¥ oo. For a 
nonzero correlation length d of Gaussian (in the sense 
of fluctuation statistics ) diso rder, it was recognized by 
iRaikh and Shahbazvani (|l993f ) that the SCBA is a para- 
metrically accurate approximation for N 3> 1 provided 
d <C A. The hard gaps in between are an artefact of the 
SCBA, but the exponential tails of the DOS which fill 
the gaps smear on ly slightly (for N ^> 1) the boundaries 
of th e semicircles ( Benedict] . Il987t lEfetov and Marikhinl . 
Il989l). The gaps predicted by the SCBA close at uj c t = 
tt/2 (|Andol . Il974bt ILaikhtman and Altshulerl . 11994) . 

Apart from the conditions N 3> 1 and d -C A, the 
SCBA also assumes that disorder is weak; specifically, 
that at zero B the conventional Born approximation is 
valid, which for Gaussian disorder means 2 d <C vpT q . 
Moreover, the SCBA crucially assumes that the non- 
Gaussian component of fluctuations of the random poten- 
tial can be neglected. The Gaussian character of disorder 



2 According to Eq. (0 with W q from Eq. J5]l and n e = m, this con- 
dition can only be marginally satisfied for randomly distributed 
remote impurities. 
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implies that individual impurities are weak and their con- 
centration rii obeys rij max{d 2 , S} ^> 1, where S is the 
area over which disorder is averaged because of the quan- 
tum uncertainty of electron paths. For separated LLs, 
the characteristic S is given by the area 2ttX 2 occupied by 
one electronic state, which is assumed within the SCBA 
to be much larger than d 2 . The DOS for disorder with 
non-Gaussian statistic s of fluctuati ons wa s considered 
beyon d the SCBA by lAndol (Il974alibh by iBrezin et all 
for N = 0, and by iBenedict and Chalkerl (|l986h 



for separated LLs in the limit iV 3> 1. In particular, it 
was found in these works that within a model of delta- 
function impurities there develops a singularity in the 
DOS at e — en for 27m; A 2 < 1, independently of the 
strength of the impurities. 

If uj c T q <C 1, Landau quantization only leads to a weak 
modulation of the DOS and the SCBA predicts that 
v^ N ' (e) in this limit is described by a Lorentzian with 
ImS(e) = — l/2r q . The DOS is then represented, by 
means of the Poisson summation formula, as a sum over 
harmonics in the following form: 



2nke 



cos ■ 



fe=i 



where the Dingle factor (jShoenbergl . Il984j ) 
<5 = exp(-7r/w c T q ) . 



(29) 



(30) 



The first harmonic — the least damped — yields the lead- 
ing oscillatory correction Av(e) to the DOS: 



Av(e) = — 2vqS cos(2tts/uj c ) . 



(31) 



For d> A, the SCBA fails and instead the quasiclas- 
sical approximation in which the DOS is represented as 



^ N \e)^(2nX 2 )- 1 (5[e 



V(v))) (32) 



is valid ( Raikh and Shahbazvani Il993l) . Here 



rZTT in 

V(r)= —V{x + R c cos0,y + R c sm6) (33) 
Jo 2?r 

is the effective random potential acting on the guiding 
center of the cyclotron orbit. For the case of Gaus- 
sian statistics of fluctuations of V (r), the shape of 



z/_j (e) is also Gaussian and i s given (jMirlin et al. 
iRaikh and Shahbazvani Il993l) by 3 



lapeot 
jl ll996t 



AN) 



2vrA 2 



1/2 



exp 



(34) 



for A < d < R c and by 



1 



2ttA 2 



[27rW(0)] 1/2 exp 



(e - s N ) 2 
2W(0) 



(35) 



for d > R c [in which limit V(r) ~ V(r)]. In Eq. ([3"5]>. 
W(0) is the variance of V(r) [Eq. ©]. Note that the 
broadening grows with increasing B and saturates at 
d^> R c . Equations (|34| and ([35]) describe in effect inho- 
mogeneous broadening and are valid for both separated 
and overlapping LLs. 

The sum of (e) from Eq. AMI) over N gives the sum 
over harmonics 



l + 2^(-l) fe J fc2 cos 



2irke 



fc=i 



(36) 



Note that the only difference between Eqs. ([29]) and (|36|) 
is in the power of the factor 6: it is k in the former case 
and k 2 in the latter. It follows that the damping of the 
leading (k = 1) oscillatory correction in the limit of over- 
lapping LLs is given by Eq. (|3"T|) in both cases, despite the 
shape of LLs being quite different (Lorentzian vs Gaus- 
sian). If d 3> R c , the exponent of the damping factor 
for the fcth harmonic is proportional to fc 2 , similarly to 
Eq. (|36[) , with the k — 1 term in vie) given by 

Ai/(e) = -2^ exp[-27r 2 LT(0)/o; 2 ] cos(27re/o; c ) . (37) 

The Dingle plot [the term conventionally used to describe 
the behavior of In Ai/(er) plotted in the low-T limit as a 
function of 1/B] is seen to be quadratic in Eq. (|3T|) . in 
contrast to Eq. (|31l) where it is linear. The quadratic 
Dingle plot 4 in the DOS can thus be used to probe the 
presence of long-range (with the correlation radius d > 
R c ) inhomogeneities in the sample. 



2. Shubnikov-de Haas oscillations 

The magnetooscillations of the DOS give rise to the 
oscillatory behavior of the dc transport coefficients as 
Sf/^c is varied — t he Shubnikov-de Haas (S dH) effect. 
Within the SCBA (jAndo and Uemural . 11974 . the con- 
ductivity tensor for N ^> 1 was obtained for short- 
range (d <C kp 1 ) disorder by lAndo (Il974b ) f or the di- 
agonal conductivity a xx and by Ando etaU ( 19751 ) for 



the Hall conductivity a xy . Note that the use of the 
SCBA to treat also long-range dis order with d A in 
the early wor ks (|Ando et all Il982f) is not justified since 



this condition violates the applicability of the SCBA 



3 Equation H34II was also obtained by Gcrhardt^ dl975T ) within a 
cumulant expansion of the self-energy J26j in the time domain 
to first order in W q . However, the control parameters of the 
approxim ation were not cl early specified; in fact, Eq. d 341) was 
argued by Gcrhardts ( 1975) to describe short-range disorder with 
d <C A, contrary to the actual condition which is the opposite. 



4 If the damping of the oscillations in the DOS is solely due to re- 
mote impurities, a finite range of B within which LLs are not sep- 
arated and the Dingle plot is quadratic exists provided disorder 
is sufficiently strong; specifically, if the "out-scattering length" 
VFTq <C d. In ultra-high mobility samples, vpr q is typically of 
the order of d. 
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(jRaikh and Shahbazvanl Il993l ). A generalization of the 
SCBA approach which provides a framework for studying 
transport in the case of disorder with a n arbitrary cor- 
relatio n length cl -C A was developed by iDmitriev et all 
(2003). An essential ingredient of the theory describing 
long-range disorder is the inclusion of vertex corrections 
in averaging over disorder. The result is that the the- 
ory can be formulated solely in terms of the oscillatory 
DOS v(e) and the transport scattering rate tb(e) which 
is renormalized by Landau quantization as follows: 

1 - lv ® (38) 



r B (e) 



Expressed in terms of v(e ) andrgfe), o xx (neg lecting the 
Zeeman splitting) reads (jDmitriev et aH. 120031) 



ds 



de 



(39) 



l + ^ri( £ ) ' 

where is the thermal distribution function. Equation 
(l3"9")l describes both separated and overlapping LLs and 
is valid also in the crossover in between. Note that, com- 
pared to the Drude formula ctd = (e 2 n e r/m)/(l + w 2 t 2 ), 
Landau quantization manifests itself in Eq. (1391) in two 
ways: in the renormalization r — > tb(e) and in the ap- 
pearance of the additional factor v(e)/i>q in the numera- 
tor [where, in view of the relation (|3"8"]) . the DOS oscilla- 
tions cancel out]. As B is varied, a xx at zero T is repre- 
sented in the limit of separated LLs as a series of peaks 
the height of which is given by cr"^ ax — CTD(2oj c T q /7r) >• 
<td- For overlapping LLs, the oscillatory correction Act 
to ctd reads 

^,^,^7^ , (40) 

CTD 1 + WJT Z Vq \ UJ C J 

where Az/(e) is given by Eq. (|3"Tj) and the factor F(x) = 
x/ sinhx describes the thermal averaging of the oscilla- 
tions. Importantly, in the case of long-range disorder 
(k F d ^> 1), the magnetic field may be "classically strong" 
in the sense of w c r ^> 1 and at the same time, if o; c r q <C 1, 
lead to only weak SdH oscillations even at T = 0. Note 
also that for T ^> lu c ^> t^^t -1 the thermal averaging 
over the contributions of separated LLs in Eq. (|39|) leads 
to a nonoscillatory contribution to the MR 

pxx/pD = (8/3^)(2w c r q /^) 1 /2 . (41) 
Extending Eq. (|3"8"]l to the Hall conductivity yields 

df?\ He) 



en e c e 2 v 2 F . 

a xy = ?r" H ' de 

B uj c 



de J l+u, 2 r 2 ( £ ) 



(42) 



where the first term describes collisionless drift of elec- 
trons in crossed electric and magnetic fields. The oscil- 
latory MR for overlapping LLs is then given by 5 



Ap _ 2 &v{e F ) ^ItPT 



(43) 



for arbitrary u> c t <C T/r q . It is worth noting that the 
main contribution to the oscillatory MR (|4"3")l comes from 
the oscillations of a xx for u> c t 3> 1 and from the oscil- 
lations of a xy for uj c t C 1, In the case of short-range 



disorder ( t = r q ) , Eqs. (|39|) and (|421) agree with those 
derive d bv lAndol (|l974b . Eq. (2.12)) and bv lAndo et~al\ 
(1975, Eq. (3.26)), respectively. 6 In the case of long- 
range disorder, t hey confirm the f orm o f SdH oscillations 
hypothesized bv lColeridge et al\ ( 19891) . 

It is worth noting that the Fermi energy ef for a 
fixed electron density n e also exhibits magnetooscilla- 
tions. These are qualitatively important for the depen- 
dence of the MR on B in the case of separated LLs at 
zero T: it is because of the oscillations of ep that p xx as 
a function of B does not have gaps similar to those in the 
dependence of v(e). For overlapping LLs, the oscillations 
of Ef only yield a subleading (in lo c /ef *C 1) oscillatory 
term in the MR. The thermal averaging for T 3> w c sup- 
presses the oscillations of ef exponentially, similarly to 
Eq. (00]). 

As discussed in Sec. III. Al disorder in high- mobility het- 
erostructures is a mix of long- and short-range compo- 
nents characterized by vastly different correlation radii. 
Because of the large ratio T/r q for scattering off the long- 
range component of disorder, it is possible that the mo- 
bility fi and the Dingle factor 5 are determined by differ- 
ent sources of disorder. Specifically, it is likely that, in 
ultra-high mobility samples, p is limited by background 
impurities and interface roughness, whereas the damping 
of the magnetooscillations is mainly due to scattering off 
remote impurities. 

Note also that, in real samples, there may exist ultra- 
long range inhomogeneities with a correlation radius 
larger than R c in the regime of SdH oscillations. These 
inhomogeneities may not affect the mobility but will pro- 
vide an additional damping of the oscillations of A.p/p in 
Eq. (|4"3"1) . which is described by the exponential factor in 
Eq. d37j). If W(0) > w c /r 9 , the damping of the SdH os- 
cillations will be mainly determined by these macroscopic 
inhomogeneities and the Dingle plot for the MR will be 
quadratic in 1/B, a s was indeed reported for some SdH 
measurements (e.g., IColeridgel Il99ll Il997l) . 

So far, the magnetooscillations of the MR and the DOS 
have been directly related to each other, being character- 
ized by the same Dingle factors. One important reason 
why this is not true in general is the suppression of the 
transport coefficients by localization effects. These come 
in two varieties: classical and quantum. An example of 
classical localization is the adiabatic drift along equipo- 



5 Inelastic electron-electron scattering (or, for that matter, any 



scattering whose strength is proportional to e 2 + 7r 2 T 2 ) does not 
lead to an a dditional exponential damping of the SdH oscillations 

for r>u c jAdamov et all |2006| ; iMartin et aH. |200J). 

The e xpansion of a xy to first order in 8 obtained by Ando et al. 
lll975f) was also derived by means of a diff erent representation 
of the Kubo formula for Hall transport by llsihara and Smrckal 
jl986t) . 
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tential lines of a smooth rando m potential, which results 
in the exponential suppression ( Fogler et q?J .|l997) of a xx 
with increasing B [Eq. (fTTj) ] . Oscillations of the MR are 
then associated with the oscillatory DOS for only a subset 
of electron paths; specifically, for those giving the main 
contribution to a xx ■ The statistical properties of disorder 
present along the "conducting" paths do not necessarily 
coincide with those on average, which leads in general to 
different Dingle facto rs for the DOS an d the MR. The 
difference was shown ( Evers et all Il999t) to be very pro- 
nounced for the case of a spatially random magnetic field. 



Quantum localization becomes strong when a xx drops 
with increasing B to a value of the order of the conduc- 
tance quantum e 2 /2tt, at which field the QH effect sets 
in. 7 Importantly, if remote impurities, distributed ran- 
domly with the density rii = n e , are the only source of 
disorder, the cyclotron frequency for the crossover field 
Bioc obeys w c T q ~ (fc^d) -1 / 2 <C 1, i.e., strong localiza- 
tion develops when LLs are still well overlapped and the 
DOS does not exhibit any gaps. If, however, the short- 
range component of disorder substantially decreases r 
while not affecting r q (which is the likely situation in 
ultra- high mobility structures, see above), B\ oc shifts up- 
ward. It follows that it is only because of the admixture 
of sufficiently strong short-range disorder that developed 
SdH oscillations can be observed in high mobility samples 
before the crossover to the QH regime. 



At this point, it is also worth emphasizing that quan- 
tum magnetooscillations of the MR may not be related 
to those of the DOS at all. Even if the latter are ne- 
glected, the quantum interference of diffusive waves for 
&xx e 2 /27r gives rise to the QH oscillations Aa xx cx 



xxi, where 

Qxx,xy IS C XXl xy in units 



cos(2itq x v ) exp(-27 rg_ 

of e 2 /27r (lPruiskerJ.ll990h . These are the oscillations into 
which the QH effect transforms as B is decreased — not 
into the SdH oscillations, contrary to the common belief. 
Note that the period of the QH oscillations corresponds 
to one flux quantum through the area min{/ 2 , R 2 }. That 
is, the -B-dependent interference occurs within the area 
the size of which is not related to (and is smaller than) the 
localization length, so that the observation of these oscil- 
lations for a xx e 2 /2ir does not require e xponentially 
low T — in contrast to the QH plateaus. See lMirlin et al\ 
(|l998l) and lEvers et al\ ( 19991 ) for a discussion of the in- 
terplay of the SdH and QH oscillations in the case of a 
random magnetic field. 



3. Quantum magnetooscillations of the ac conductivity 

The oscillations of the DOS induced by Landau quan- 
tization have important consequences also for ac trans- 
port. Specifically, the Drude formula in Eq. (|22|) for 
a±(iv), which describes a single peak of the CR, and 
Eq. (1591 for the oscill atory a xx at lu = generalize 8 to 
ilDinitricv et all 120031) 



J„,2 



<T±(U>) 



2lo 



x de ■ 



! (e)+T B 2 (e + cu)}/2+{uj±uj c y 



(44) 



At zero T, for the Fermi energy lying between two sepa- 
rated LLs, Eq. (j44|) gives a CR peak whose height 9 



<7-{u c ) = {e 2 n e T/m)(2/iruj c T q ) 1/2 



(45) 



is much smaller than that following from the Drude 
theory [Eq. (|22p]. In the case of long-range disorder 
(t ^ r q ), the CR lineshape from Eq. (|4"4")l is given by 



a (uj) _ f ^ .j, 



*CR 



where 



and 



/(*) = 1 - 



2-y/T 



An 



A C r = (2 Wc r q / 7 r) 1 / 2 /r 



(46) 



(47) 



(48) 



The width of the CR for separated LLs Acr is seen to 
be much larger than 1/r. At the same time, in the case 
of long-range disorder, the CR as a function of w is much 
narrower than the peaks in the DOS. 

Landau quantization in combination with disorder 
leads to the emergen ce of CR ha rmonics at ui = Mu> c 
with |M| = 2,3,... (|Andol . [l975h . For separated LLs, 
Eq. P4]) gives for the height of the Mth peak 



o-(Mw c ) 



4 e 2 n e 
37r muj r 



J c'<i 



1/2 



M 



f2 



1 



r (M 2 - l) 2 



(49) 



which is seen to be smaller by the factor of T q /r than the 
result (jAndol . ll975T) for the case of short-range disorder. 



7 Strictly speaking, single-particle states of a 2DEG at B ^ are 
localized except those at a discrete set of critical energies — this 
localization constitutes the essence of the QH effect. However, 
the temperature needed to observe localization is exponentially 
small if <j xx 2> e 2 /2n in the absence of localization, so that in 
reality the QH effect becomes pronounced only when <r X x per 
spin drops down to the values of a few conductance quanta. 



Note that Eq. 144li differs i n an e s sentia l way from the expres- 
sion for <j± (w) obtained by lAndol lll975h (for the case of short- 
range disorder). In particular, the e-dependent factors in the 
integrand of Eq. 1144 \ do not reduce near the CR to the simple 
product (/jT — fT+,.,) y[ e ) u (£ + t^) m * ne limit of separated LLs, 
as suggested bv lAndol U975I') for the C R lineshape. 
The peak height found by Ando (1975) for short range disorder 
is larger than that given by Eq. i 14">l) at r = r q by a factor of 4/3. 
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CO / co c 

FIG. 3 Magnetooscillations of the ac conductivity of a 2DEG 
with smooth disorder (r/r q = 10) in the limit T 2> Td, in 
units of the dc conductivity, as a function of uj/u c for fixed 
WcTq/TT = 3.25 (solid line) and 1 (dashed). Inset: the con- 
ductivity vs uj c /u for fixed uit q /2tt = 1 in units of the ac 
conductivity at B — 0. From [Pmitriev et all l[2003l ). 



In contrast to the CR peak, the width of the higher- 
harmonic peaks is given by the width of the peaks in the 
DOS. Thus, if one imagines that the correlation radius 
of disorder is increased at fixed r, the CR peak becomes 
higher and narrower, while the tendency for the higher- 
harmonic peaks is the opposite. 

For strongly overlapping LLs, the dynamic response 
of the 2DEG is given by the CR, which is described for 
w c r q <C 1 by the Drude formula (|2"2"1) . and by the in- 
dependent analog of the SdH oscillations (|4"0"[) . To first or- 
der in 5, the oscilla tory correction A<7+ (u i) to the Drude 
conductivity reads ( Pmitriev et all 12003) 



Acr ± (u) 
CD,±M 



Aj/(e F )^_/27r 2 T 
u \ Ul c 



Xi 



where 



Xi 



n 2 ± 



1 



TTUJ 



27TCJ 1 



sin ■ 



3i^_|_ Q 

sin^ 



(50) 



TTUJ , 

— 51 



and f2± = (lu±lu c )t. In the dc limit, Eq. ([SI?)) reduces to 
Eq. ([401) . If T is much larger than the Dingle temperature 
Tjj = l/27TT q , the main source of damping of the oscil- 
lations with varying ep/ui c is the thermal averaging. In 
ultra-high mobility heterostructures, Tjj ~ 10 2 mK and 
for the typical measurement temperature T ~ 1 K the 
linear- m-S oscillations (f50[) are completely washed out. 

Unlike the oscillations with £f/uj c , which constitute 
the essence of the SdH effect in the dc case, the oscilla- 
tions with lo/ui c survive the thermal averaging. Specifi- 
cally, they survive in the even-order terms in the expan- 
sion of a± (ui) in powers of S, in which the intermodulation 
of the DOS oscillations at energies separated by u) docs 
not go away upon averaging over energy. For strongly 
overlapping LLs, the oscillatory behavior of <t±(w) as 
Ld/uJc is varied at T > Tjj comes from the correction 
Aa±(ui) of order 5 2 , which in the high-T limit is given 



by ([Dmitriev et all . l2003f) 

Aer±(w)/er D ,±M - 



(52) 



with 

X2 



n 2 ± 



2ttuj 



(O^ - 3) cos + 



3^4 - 1 



2ttuj 



■ sin ■ 



(53) 



As noted already in the discussion of the SdH effect, in 
the case of long-range disorder, for a classically strong 
magnetic field with uj c t 1, the modulation of the DOS 
is still small as long as w c r q -C 1. Then, away from the 



CR (for |0±| > 1) Eq. ([52]) reduces to the simple form 

A<t±(w)/c7d,±M ~ 2<5 2 cos(27rw/o; c ) . (54) 

The overall behavior of the T-independent oscillations of 
<r(u>) is illustrated in Fig. [3J 

Intimately related to the independence of Eqs. (|52p 
([54]) on T is the fact that the oscillations with lj/oj c 
arc not sensitive to the presence of macroscopic inho- 
mogeneities (which may complicate the determination of 
T q from SdH experiments, as discussed at the end of 
Sec. III. C. 21) . The measurement of the damping of the 
u) /w c -oscillations as a function of B in the high-T regime 
may thus be the most reliable means of extracting r q from 
the magneto-oscillatory behavior of the transport coef- 
ficients. Note, however, that the CR harmonics — and 
the related periodic modulation of cr±{ui) with ui/ui c — 
may be associated not only with Landau quantization 
but also with the quasiclassical memory effects discussed 
in Sec. III. B. 31 In fact, the latter can produce strong os- 
cillations of <j±(ijS) with an amplitude comparable to the 
Drude conductivity ctd,±(w) [see Fig. [5J for more details 
of the comparison between these two sources of the mag- 
netooscillations in the ac conductivity and the dc photo- 
conductivity see iDmitriev et ~al. I (|2004l) ]. 

The oscillations of u±(uj) as ui/ui c is varied on the 
low-i? side of the CR, as well as the oscillations of 
the dynamic conductivity with ef lw c , have been ob - 



served in the early studies by lAbstreiter et all ( 1976[) 
in the far-infrared light transmission experiments on a 
low-mobility 2DEG in Si inversion layers. However, 
only recently both these types of magnetooscillations 
have also been observed in microwave absorption exper- 
i ments on high-mobilit y GaAs/AlGaAs heterostructures 
(|Fedorvch et all , l2010h . One reason why the measure- 
ments of <t±(ll>) in the absorption/transmission experi- 
ments may be substantially complicated is that the con- 
ductivity expresses the current as a response to the total 
(screened) electric field whereas what is probed in this 
type of experiments is a response to the (unscreened) 
field of the incident electromagnetic wave. As a result, 
the behavior of u±{ui) in the absorption/transmission 
may b e masked by t he ex c itation of magne t oplas- 
mons dFedorvch et all l2010l : IStudenikin et all , l2007t 



IWirthmann et all , 12007) . 
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FIG. 4 Absorption coefficient A, differentiated twice with 
respect to B, as a function of uj c /uj for fixed to. Three lower 
curves, marked by different radiation frequencies ui/2n, were 
taken at T = 2 K in a sample with n e — 3.6 x 10 11 cm" 2 
and fi ~ 5 x 10 6 cm 2 /Vs. The upper curve is the calcu- 
lated dependence of —d 2 A/dB 2 on oj c /lj at fixed uj according 
to Eqs. Q54p and (|55|) with r q = 9.1 ps determined from the 
oscillatory-photocon ductivity experiment. The curves are off- 
set for clarity. From iFedorvch et al\ (|2010l ') . 



the limit where the Drude formula for g± can be applied, 
2 7 r 7T 



A D .± = 



1 



T D ,± = 



(l+ 7 r) 2 +^ 



(56) 



and the dynamical-screening-induced broadening 

7 = 2ire 2 n e /mcn r . (57) 

In ultra- high mobility GaAs/AlGaAs heterostructures, 
the product 7 t may be as large as a few tens. A sub- 
stantially enhanced — compared to the transport scatter- 
ing rate — broadening of the CR in the absorption (for 
magnetic fields in which the LLs were not much sep- 
arated) was reported in a high-mobility structure by 
IStudenikin et al. (2005), similar to the very strong broad- 
ening of the CR in the transm i ssion through ultra-hig h 
mobility structures (|Smet et aLl . l2005tlTung et alL\20(M . 
The oscillatory behavior of the absorption coeffi c ient a s 
a function of U)/u) c , observed by IFedorvch et al] (|2010h . 
was related therein to the oscillations of the conductiv- 
ity tensor through Eq. (|5"5")) . The \ow-B portion of the 
data shows good agreement with the simple asymptotic 
behavior resulting from Eq. ([54|) . see Fig. |4j 



One more point — specific to a high-mobility 2DEG 
with \a xx (tu) ± i<j xy (u)\ ^ cn r , where c is the speed 
of light in the vacuum and n r is the refractive in- 
dex of the medium on the sides of the 2DEG — is a 
strong reflection of the electromagnetic wav e caused by 
the dynamical screening inside the 2 DEG (jChiu et all 
119761 : iFal'ko and Khmel'nitskil Il989f ). In a magnetic 
field, the enhanced reflection leads to a suppression and 
an additional br oadening of the CR ((Ch iu et al ., 1976; 
iMikhailovl I2004J ) in the absorption/transmission coeffi- 
cient. For a linearly-polarized wave normally incident on 
the 2DEG, the absorption coefficient A = A + + A_ and 
the transmission coefficient T = T + + T_ are given by 10 



4. Intersubband magnetooscillations 



A± = Rc.g±/|l + .g±| 2 , T± = 1/2|1 + g±\ 



(55) 



where g± = 2Tr(a xx ±ia xy )/cn r (neglecting the difference 
in the dielectric properties of GaAs and AlGaAs). 11 In 



10 Note that Eqs. (J35J describe scattering by the 2DEG itself. To 
relate them to the absorption and transmission coefficients of 
the electromagnetic wave incident on the sample, one should 
take into account multiple Fresnel reflections from the bound- 
aries of the latter. Changes of the refractive index on scales 
smaller than the wavelength of the radiation are averaged. If 
the wavelength inside the sample 2irc/um r is much larger than 
the thickness of the sample, one can neglect the polarization of 
the medium around the 2DEG and put n r — > 1 in g± and 7. 
When the distance from the 2DEG to the top surface of the het- 
erostructure is much smaller than the wavelength (which is a 
typical situation), sca ttering by the 2DEG and the top surface 
l|Fedorvch e t al, 2010) is described by Eqs. l|55|l with the change 
n r — ¥ (1 + n r )/2 in g± with a simultaneous multiplication of A± 
and T± by 2n r /(l +n r ) and 4rir/(l + « r ) 2 , respectively. 

11 To avoid notational confusion: the abbreviation cr± is used 



The mechanism in Sec. III.C.3I that relates the oj/lj c - 
oscillations of <j(uj) — which survive the thermal averag- 
ing [Eq. (|52p] — to the intermodulation of the DOS os- 
cillations at energies separated by cj can also lead to 
T-independent magnetooscillations in dc transport. For 
that, one needs two (or more) parallel 2D electron gases, 
allowing for the exchange of electrons between them, with 
the bottoms of their energy bands located at different en- 
ergies. One system of this type is a QW wi th two lowest 
subbands o ccupied by ele c trons , for which iPolvanovskvl 
(|l988l ) and lLeadlev et al] ( 1992t ) pointed out that there 
should exist oscillations of the MR periodic in A/w c , 
where A is the subband spacing, not suppressed by the 
thermal averaging. An SCBA theory of SdH oscillations 
in the two-subband conductivity t ensor for an ar bitrary 
type of disorder was developed by iRaichevI (|2008h . 

The correction Acr m i so to the dissipative Drude con- 
ductivity, which describes the magnetointersubband os- 
cillations (MISO), is a sum of two terms Acr m i SOiQ with 
a = 1,2, coming from electrons in the ath subband. 
For uj c T a -c 1, the partial contr ibutions Aa m i s? a are 
given for whit e-noise disorder by ( Averkiev et txL 120011 ; 
IRaichevI . l2008h 

2e 2 n a rl ( 2r a \ 2ttA 
Acr miso , Q = 1 did 2 cos , (58) 

TOT12 V r 12 / u c 



in this review to denote Re (a XJ: ± i<r yx )/2 (as, e.g., in 
iDmitrie v et all |200 3ft , whe reas in a number of relevant papers 
(e.g., Fcdorych et al., 2010) <j± denotes cr xx ± icr yx . 
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where n a and 1 /t q are the electron density and the zero- 
B scattering rate in the ath subband, l/r a includes both 
intra- and intersubband scattering, l/ri2 is the zero-i? in- 
tersubband scattering rate, and S a is given by Eq. (|30[) 
with r a substituted for r q . Equation (|58|) follows 12 di- 
rectly from Eq. (|39|). Importantly, Eq. (|38|) does not 
hold for the relation between the oscillatory DOS v a (e) 
and the oscillatory scattering rate rs. Q (e) in the ath sub- 
band, since l/rB jQ (e) contains the rate of intersubband 
scattering, whereas v a {e) does not (in the leading ap- 
proximation in 1/ 'ef.oTyi-, where ep, a is the Fermi energy 
in the ath subband counted from its bottom). 13 This 
is the reason why the preexponential factor in Eq. (|58|) 
does not vanish at B — > 0, in contrast to Eq. (|52|) at 
uj = 0. As a result, the MISO contribution to the MR 
/p D — — Ao- miso( o D at co c T a <C 1 is mainly given 
by the oscillatory dissipative component of the conduc- 
tivity tensor (in contrast to the oscillatory dc MR in the 
single-subband case, Sec. III.CT2j) . 

Within the SCBA, Eqs. (J2S1) and (g^J) allow one to de- 
scribe MISO for long-range disorder as well. Similar to 
Sec. III.C.21 the scattering rate in the exponent of S a re- 
mains then the quantum scattering rate, while 1/Ts ja (e) 
becomes the momentum relaxation rate. However, in 
general, if scattering is not isotropic, the zevo-B momen- 
tum relaxation rate in the ath subband 1 /t q is related to 
the intra- an d intersubband sc attering rates in a nonlin- 
ear manner (|Ando et all Il982l) . It is only in one simple 
case, when the intersubband scattering is isotropic, that 
1 jr a is still given by a sum of the momentum relaxation 
rate for (not necessarily isotropic) scattering within the 
subband and the intersubband scattering rate l/ri2- In 
this case, 14 Eq. (|S"5|) for uj c T a <C 1 is reproduced also for 
long-range disorder. 

The cumbersome expressions for the MISO amplitudes 
for th e general case of anisotropic scattering 



(Raichcv, 2008) simplify significantly in the limit of a 
classically strong magnetic field (w c t q ^> 1) but overlap- 
ping LLs (6 a <C 1) — which conditions are only compati- 



12 The MISO amplitude res ulting from Eq. (139 I t and the one ob- 
tained by iRaichevI {2008) agree wi th each other f or arb itrary 
u) c T a , whereas the one obtained bv lAverkiev et al. I (l200lh coin- 
cides with them only for u) c r a — > 0. 

13 At order 0(6162), there are two contributions to MISO for 
u c T a <C 1: one comes from the thermal averaging of the cross 
term between the oscillations in v a (e) and those in Tg Q (e), the 
other — from the averaging of tb a( e ) alone. In Eq. H58H . the for- 
mer brings —1 to the expression i n the brackets, the latter bring s 
2(1— T a /ri2). In an early work by Raikh and Shahbazyan (1994), 
where the SCBA approach to describe SdH oscillations in tjxx 
was generalized for short-range disorder to the case of two sub- 
bands [see also an extension of the the ory to essentially similar 
oscill ations in biased bilayer graphene (Mkhitaryan and Raikh, 
l201lf) 1. only the former contribution was taken into account. 

14 The isotropy condition for intersubband transitions requires that 
the Fermi wavelength in each of the subbands be much larger 
than min{«), d}. As discussed in Sec. III. Al this condition is only 
marginally satisfied in high- mobility GaAs/AlGaAs structures. 



ble in the case of long-range disorder. In this limit, 

2e 2 (n 1 + n 2 ) x 2vrA 
A<T miso = 7TL c>i o 2 cos (59) 

and, correspondingly, the MISO term 15 in the MR 

Apmiso 2 n 1 Ti+n 2 T 2 z S . 2ttA 

= 0id 2 cos , (60) 

Pd T12 ni + n 2 lo c 

where l/fi2 is given by an integral over the scattering 
angle <p [cf. Eq. ©]: 

^=-f ?(l-^^cos^W-), (61) 
T12 J 2vr V n x +n 2 V q V ' 

(12) 

with Wq being the Fourier component of the correla- 
tion function of the intersubband matrix element of the 
random potential at the transferred in-plane momentum 
q = [2to(£f,i + £f,2 - 2 % /£f,i£f,2C0S (f))] 1 ^ 2 - 

Oscillations similar to MISO occur also in a double QW 
with the lowest subband states split by tunneling or, for 
that matter, in any other system whose energy spectrum 
exhibits two or more series of LLs, which are offset with 
respect to each other, in the presence of scattering be- 
tween the series. Experimentally, MISO were studied in 



2010a; 



1992, 



two-subband sin gle quantum wells (iBvkov et al. , 
I ColeridgeL [19901 iGoran et a/.l. l2009t iLeadlev et al. 
Il989l; iRowe et all 120011: ISander et all Il998l) . in dou- 
ble dBvkovl 12008b!; iMamani et all 120081) and triple 
( Wiedmann et al. , 2009bl) Q Ws, and in wide QWs spli t 
in two layers electrostatically ( Wiedmann et all l2010ah . 
The behavior of MISO in nonequilibrium conditions is 
discussed in Sees. IIII.CT41 and IVTU1 



III. MICROWAVE-INDUCED RESISTANCE 
OSCILLATIONS (MIRO) 

A. MIRO: Experimental discovery and basic properties 

When a sufficiently high-mobility 2DEG is subject to a 
weak magnetic field and illuminated by microwave radi- 
ation, the longitud inal magnetoresistivity p (B) exhibits 
giant oscillations (jZudov et aU ll997. 200llJ, termed the 
microwave-induced resistance oscillations (MIRO), see 
Fig. [5] In a Corbino disk-shaped 2DEG, experiments 
revealed cor responding microw ave- induced conductance 
oscillations ( Yang et all . 120031 ) . Most commonly, MIRO 
are observed in dc measurements, with radiation be- 
ing delivered to the 2DEG via oversize waveguides (al- 
though other means, such as planar microw ave trans- 
mission lines patterned on top of the 2DEG ( Ye et all 



5 The expression for Ap m i so in terms of the thermall y averaged 
produ ct v(e)v(e+ A) yielding Eq. J 60 1) was proposed bv lColeridgd 
ijl990T l. although it was erroneously also suggested there that 
MISO arc suppressed by the thermal averaging if T 3> lo c . 
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FIG. 5 Magnetoresistivity with (solid lines) and without 
(dashed) microwave irradiation for different frequencies / = 
45, 72, and 94 GHz, normalized to its value at B — 0. The 
traces for different / are vertically offset for clarity. Integers 
show the order of MIRO peaks. The arrows mark the magne- 
toplasmon resonance. The data were obtained at T ~ 1.7K in 
a 200 )im wide Hall bar sample with n e — 2.0 x 10 11 cm" 2 a nd 
(i~3.0x 10 B cm 2 /Vs. Adapted from lZudov et al\ i|200ibl ). 



l200lh and dipole antennas (|Willett et all l2004) have 
also been succes sfully implem e nted). Rece ntly, they were 
also o bserved bv lBvkov et all ( 2010bl ) and lAndreev et al\ 
in contactless measurements in a capacitively- 
coupled 2DEG. MIRO have been studied not only 
in high-mobility single-subband 2DEGs but in a vari- 
ety of other 2D systems, inclu ding 2DEGs patte rned 



with a triangular antidot lattice (lYuan et all 2006) 
Sec.IIII.C5l two-subba nd 2DEGs ilBvkuv ct al. 



sec 



2010d: 



IWiedmann et 120081) . see Sec. IIILC.5l and hole sys- 
tems based on C-doped GaAs/AlGaAs QWs (|Du et all 
2004). Phenomenologically similar oscillations were re- 
cently discovered in a nondegenerate electron system on 
surface of liquid 3 He ( Konstantinov and Konoll2009l ). see 
Sec. IVIIB.lt 

Period and phase. The radiation-induced oscillatory 
part 5puj of the photoresistivity p u at a radiation fre- 
quency uj oscillates with the ratio 



(62) 
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FIG. 6 (a) Magnetoresistivity oscillations with varying B at 
fixed / = 81 GHz. (b) Microwave-induced correction Spu, to 
the resistivity as a function of e ac , obtained by subtracting the 
slowly varying background. The data were obtained at T ~ 
1.5 K in a 100 pm wide Hall bar sample with n e ~ 3.0 x 10 11 



cm" 2 and p ~ 1.2 x 10 7 cm 2 /Vs. 



i and occur at 



n=F ^ ac , n = 1,2,3, 



(63) 



In accordance with theoretical predictions for the regime 
of overlapping LLs (Sec. lIII.B.l.a[) . the observations sug- 
gest that Spu oc — sin27re ac for e ac > 2, i.e., the phase 
ip ac in Eq. (J63J) is 1/4. This value of the p hase was re- 
ported as un iv ersal i n many exp e rimen ts ( Manil. 2004 , 
2004 l2007a|| b|. 120081: iMani et all l2010l 120091 l2004aUb . 



2002, 20 04dldD. O th er expe ri ments (IDai et al 
Hatke et a/.I.T2008allbl. feoilbllcl [Studenikin et al 



2010: 



2005 



20071: iZhang et all l2007al IzTidovl 12004 IZudov et al , 
200630)" found that the lower-order maxima and min- 
ima are pushed towards the harmonics of the CR 
and are best described by <y9 ac which is considerably 
smaller than 1/4. Such a phase r eduction is ex - 
pected in the regime of separated LLs (|Dmitriev et al. 



2005: St udenikin et al. , 20051 IVavilov and Aleinerl . 12004 
Zudovl . l2004D . see Sec. lHLRT Taj Moreover, the phase de- 
creases with increasing radiation power, as discussed in 
the next paragraph. 



see Fig.[6p. In contrast to PIRO and HIRO (SecEJ), 
the MIRO period is not sensitive to the carrier density 
n e . MIRO maxima (e ac ) and minima (e ac ) are roughly 
symmetrically offset from the harmonics of the CR at 



Early experiments <|Ye et all l200ll : IZudov et all . Il997l . l2001bf ) 

loosely associated the low-order oscillation maxima with integer 
£ac = n. Now it is established that at e ac = n there are zero- 
response nodes where the photoresistivity Sp^ = 0. 
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FIG. 7 (a) Amplitude and (b) phase <p a c measured at 
the second MIRO maximum vs microwave power P. Straight 
lines in part (a) represent, on the log- log scale, the linear and 
square-root dependences on P. (c)-(e) Photoresistivity <5p„ as 
a function of ui c at fixed lj/2-k = 33 GHz for different P. The 
traces are labeled according to the attenuation levels. The 
data were obtained at T ~ 1.5 K in a 200 /im wide Hall bar 
sample with n £ ~ 2.9 x 10 1 1 cm" 2 and fi ~ 2.4 x 10 7 cm 2 /V s. 
Adapted from lHatke et all (j2011al V 



Power dependence. At sufficiently low microwave 
power P, one expects that the amplitude [as de- 
fined in Eq. (|64p below] of the oscillatory photoresistiv- 
ity Spu is linear in P. While the linear dependence o f 
A^ on P wa s clear ly observed by lHatke et all ( 2011rJ k 
IZudov et all ( 20031 ). a number of experiments found a 
strongly sublinea r dependence dMani et al. ,_ 201(i 20041 ; 



Studenikin et all . 12004 IWillett et all , 12004 lYe et al. 
20011 ). The apparently conflicting experimental reports 
can be reconciled by the existence of two distinct regimes 
of weak and strong power, characterized by the linear 
and sublinear scaling of A u with P, respectively. The 
experimental data in Fig. ~~ show that the linear depen- 
dence of A u oc P crosses over with increasing P into 
Au oc P 1 / 2 , accompanied by a strong reduction of the 
phase ip ac compared to its value (p ac — 1/4 in the linear 
regime. Such behavior is expected for the inelastic mech- 
anism of MIRO, see Sec. IIII.B~2l For the displacement 
mechanism, a similar effect is expected at the crossover to 
the multiphoton regime of the photoresponse, see Sec. lVIl 
Temperature dependence. Most experiments show that 
MIRO are best observed at temperature T ~ 0.5-1 K, get 
strongly suppressed with increasing T, and become al- 
most invisible at T ~ 4-7 K. The overall behavior in the 
regime of weak oscillations, illustrated in Fig-03 is consis- 
tent with the theoretical predictions, see Fig.ll3b.c and 
the discussion in Sees. IIII.B.l.cllTli.B.l.fl and |III.B.l.g| 
(the T dependence of MIRO in the regime of ZRS is dis- 
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FIG. 8 (a) Evolution of MIRO with varying T and (b) the 
T dependence of 27r/ii>Tq extracted fro m these data by fitting 
them to Eq. (J64j| l|Hatke et all l2009d). (c)-(d) The T depen - 
dence o f X [Eq. (16411 obt a ined by (c) lstudenikin et all ||2007D 
and (d) IWiedmann et all (|2010tJ ). 



cussed in Sec. HV.A]) . At low B and P, the B dependence 
of MIRO almost perfectly follows 



8p u = -A w e ac sin(27re ac )exp(-27r/w c T*) 



(64) 



where A M is nearly ^-independent for e ac > 2, see, e.g., 
Fig.OSK- Such fits of the B-traces of MIRO at different 
T generally reveal that both r* in the factor e _27r ^ cTc i 
and A^ depend on T. Specifically, the exponential factor 
in Eq. (|64[) is well described by a quadratic-in-T depen- 
dence of 1/t*(T) = l/r q + aT 2 , see Fig.[Sb, which is 
attributed to the LL broadening induced by electron- 
electron interactions (Sec. lIII.B.l.fj) . A similar expo- 
nential dependence on T was also reported for HIRO 
(lHatke et all l2009bl) PIRO (lHatke et all l2009al) and 



MISO dMamani et all 120081). [Studenikin e/ al 
I2007D ; IWiedmann et all (|2010bl ) found A~ 



(2005. 



oc T- z at 

T < 4 K, see Fig.~~:,d, consistent with the inelastic mech- 
anism of MIRO (Sec. IIIIB.l~c~ . The saturation of the 
depend ence of A u on T observed by IWiedmann et all 
(l2010bl) at T > 4 K, see Fig.~~i, and the results of 
lHatke et all ( 2009cf ). where A u almost independent of 
T was reported, indicate that the displacement contri- 
bution to MIRO [which produces A^ — const (T) and 
strongly depends on the correlation properties of disor- 
der, see Sec. lIII.B.l.dj may remain relevant down to low 
T. 

Frequency dependence. Over the past decade MIRO 
have been observed in a wide range of r adiation fre- 
quencies, from 3 GHz (jWillett et all |2004| ) to 1.5 THz 
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( Wirthmann et all , 120071 ). However, the majority of ex- 
periments employed frequencies from 30 GHz to 150 GHz, 
which appears to be the optimum range for a typical 
2DEG. At lower frequencies, MIRO are shifted towards 
weaker B and are therefore suppressed by the Dingle fac- 
tor. At highe r frequencies, the oscillation amplitude de- 
cays as wel l (IStudenikin et all l2007t iTung et all l2009t 
lYang et all 120031) . as illustrated in Fig. [HI for the pho- 
tocon ductance measured in a Corbino disk (jYang et all 
2003). At fixed both e ac > 2 and the microwave inten- 
sity, the observed decay is consistent with Sp^ oc w~ 4 , in 
agreement with Eqs. ([HI), ([75]), <[~~|), and ([Ml) below. 




FIG. 



9 Conductance G a 
„-2 



of a Corbino sample (n e = 3.55 x 
— 1.28 x 10 7 cm 2 /Vs) as a function of uj/oj c 
for different fixed cj. The traces, labeled accord ing to / = 
2ttuj, a re offset in steps of 0.25 mS. Adapted from Ya ng et al\ 
(120031 ) . 



B. Microscopic mechanisms of MIRO 

Most experimental findings concerning MIRO can be 
explained as a combined effect of Landau quantization 
and external fields either on the momentum relaxation 
(displacement mechanism) or on the energy distribu- 
tion of electrons within disorder-broadened LLs (inelas- 
tic mechanism). Section IIII.B.ll describes both effects 
in terms of quasiclassical kinetics of the guiding cen- 
ters of cyclotron orbits. Section IIII.B.2I deals with the 
main mechanism of MIRO saturation which establishes 
the conditions for the observation of ZRS and fixes the 
amplitude of a dc electric field in spontaneously formed 
domains. A quantum kinetic equation formalism is de- 
scribed in Sec. IIII.B~3l This formalism provides a solid 
foundation for the approach used in Sees. Sec. lIILB~T1 and 
IIII.B.2I and serves as a basic tool for the theoretical de- 
scription of other nonequilibrium phenomena throughout 
the review. Additional quadrupole and photovoltaic con- 
tributions to MIRO (these govern, in particular, magne- 
tooscillations of the Hall part of the conductivity) are dis- 
cussed in Sec. IIII.B~4l Alternative mechanisms of MIRO, 
not related to the Landau quantization, are shortly dis- 



cussed in Sec. IIII.B~~l 



1. Inelastic and displacement mechanisms 



Initially MIRO were at tributed (jDurst et all 120031: 
IVavilov and Aleinerl . |2004) to the displacement mech- 
anism which accounts for spatial displacements of 
quasiclassical electron orbits due to radiation-assisted 
scattering off disorder. Due to Landau quantization, 
which leads to a periodic modulation in the DOS 
v(e) ~ v(e + uj c ) (Sec. III.C~~) . the preferred direction 
of these displacements with respect to the symmetry- 
breaking dc field oscillates with lj/uj c . This results in 
MIRO with a phase and a period which agree with those 
observed in experiment. Photoconductivity oscillations 
governed by the di splacement mechanism were in f act 
predicted long ago ( Rvzhiil . Il970t iRvzhii et all Il986l ) in 
the limit of separated LLs and a strong dc field. This 
mechanism was further studied using various approaches 
and approximations in a n umber o f theoretical works 
(lAnderson and Brinkmanj. 20031 lAuerbach and Pal 



20071: iDmitriey et all |2009~I 

2003; Kashuba, 200 6a||bJ: iKhodas and Vavilov 



2007a; 



Durst et at 



Lee and Lcinaas, 2004; 



Parkl 



2008; 



2004; 



. Lei and Liul. 200 
Rvzhii et all 12004 IRvzhii and Surisl. 120031: IShi and Xie 
20031; iTorres and Kunoldl. l2005t IVavilov and Aleiner 



12004 IVolkov and Takhtamirovl 120071 ). Soon after the 
experi mental observations of MIRO, iDmitriev et all 
(|2003l ) proposed that the dominant contribution to 
MIRO was due to the inelastic mechanism associated 
with radiation-induced changes in the occupation num- 
bers of ele c tron s tates. Similar ideas were discussed by 
iDorozhkinl ((2003); however, the calculation there was 
not directly applicable to the experimentally relevant 
systems, in particular, due to an unrealistic model of 
inelastic relaxation. Later studies confirmed that the 
inelastic mechanism ge n erally dominates the observed 
MIRO (iDmitriev et all 120051 ) at low T, while the 
displacement mechanism can be relevant at higher T 
and only if a sufficient amount of short-range i mpuri- 
ties is present in t h e sys tem ( Dmitriev et all [2009b; 
IKhodas and Vavilovl 120081). or else, i n the limit of a 
strong dc field ( Khodas and Vavilovl 120081 ) or high 
microwave power ( Dmitriev et aLl . l2007al ). 

In this section, we formulate a description of both 
mechanisms in terms of migration of the guiding cen- 
ters of cyclotron orbits. A systematic way to obtain 
the same results within the quantum kinetic approach 
is sketched in Sec. lHLB~~l Quasiclassically, each scatter- 
ing event leads to a shift of the guiding center Ai?i 2 = 
R c e z x (n ipi - n V2 ), where n Vk = (cos tpk, sin ip k ) with 
k = 1,2 are the unit vectors in the direction of motion 
before and after the collision, see Fig. [T~J 

For a macroscopically homogeneous 2DEG subjected 
to a dc electric field E = e x E, the dissipative current 



j d = 2v Q e / dx 1 / dx 2 (W Xl 



W X2 



(65) 
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FIG. 10 (a) Shift of the guiding center of a cyclotron orbit 
due to quasielastic scattering off disorder, (b) Schematics of 
the correlations between the direction of the cyclotron-orbit 
shifts in the process of photon absorption and the sign of the 
detuning 8^ — uj/ui c — 2 for the second harmonic of the CR. 
The yellow stripes mark the DOS maxima e n — (n + l/2)u c 
in LLs tilted by a dc field. 



rates r'°^ and r' ph ) need to be specified. The calculation 
in Sec. lIII.B. 1 .el below shows that the inelastic scattering 
integral can be approximated in the form 

Sti„{A} = (jf - fe) , (69) 

which describes thermalization to the local Fermi distri- 
bution /J with the effective rate 1/T m . The scattering 
rates rS_ and for transitions ip ip' are given by 
(|Khodas and Vavilovl . [2008h 



r (ph) = 



<p+(p' • 2 f ' 

sin — 



<P 



It 



p(el) = 



1 



- 2r 



(Ph) 



The last term in T 



(el) 

<pvj 



<p<p 
■ tp—tp- 

(70) 

describes the microwave-induced 



modification of the elastic scattering. For the isotropic 
2DEG, the disorder-induced scattering rate r~ i 1 _ ¥ , 2 is ex- 
pressible in the most general case as a series in angular 
harmonics 



is expressed in terms of the probabilities W x 



of 



the guiding-center shifts x\ — > x 2 along the axis 
AX VlV2 — e x ■ Ai?i2 = R c (sm(pi — sin ip-£). The proba- 
bilities are given by the integrals over the initial and final 
energies in the scattering event 



W x 



= ( dei de 2 M ElE2 S(xi - x 2 + AX Vlip2 ) 



5(As\f±u;) 



(66) 



where the angular brackets denote averaging over the ini- 
tial and final angles tpi.2- The delta functions in the 
square brackets express conservation of the total elec- 



tron energy Aer^ = £ i 
elastic channel (oc v\p^ V2 



eEAX x 



in the 



£2 -r u.j^^j^v yi( p 2 

and its change by ±w in the 
photon- assisted scattering channel (oc T^p 2 ). Quantum 
magnetooscillations originate from the factor 



V e V e '!e (1 - fe') 



(67) 



where v € — v(e)/vq with vq = m/2ir is the dimensionless 
DOS in disorder-broadened LLs and f £ is the nonequi- 
librium distribution function in the steady state. In the 
homogeneous case, both v £ and f e are functions of the 
local kinetic energy e. Since the disorder-induced broad- 
ening of LLs is determined by r q <C t (Sec lILCTTj) , the 
effect of the external fields on v e is negligible in the rel- 
evant range of the dc and microwave field strength. By 
contrast, the modification of the distribution function is 
crucially important. The nonequilibrium occupation of 
electron states is governed by the kinetic equation 

{K^Y t JM.*±»-M,± u e) 

+ {K^'iMee' - M e , £ ))^ = St in {/ £ } , (68) 

where e 1 — e + eEAX vip i . To close the set of equations, 
the type of inelastic scattering and the expressions for the 



• <pi— <p 2 



-1 in(<pi-<p 2 ) 



n— — oc 



(71) 



For the model of two-component disorder, the coefficients 
T n are discussed in Sec. lIII.B.l.dl For the microwave field 
(screened by the 2DEG, see Sec. lIII.B.l.cj) of the form 



E ul {t)=E L 



• E± Re ( s ± e ± eiWt ) 



(72) 



where 2 1 / 2 e± — e x ± ie y and (s + ,s_) is the complex 
vector of unit length which determines the polarization 
of the field, the dimensionless power Pg in Eq. ([70|) is 
written as 

P e =V- 2Rc (£+£*_ e w ) , (73) 

V = \£ + \ 2 + \£-\ 2 , (74) 

£± = s±cvfE uj uj^ 1 (uj ±ui c )^ 1 . (75) 

We will see later that the golden-rule approach for- 
mulated above describes a rich variety of phenomena in 
strong dc and microwave fields. We first discuss it in 
the case of MIRO and calculate the linear direct current 
at high temperature, namely for 2tt T/lo c 3> 1, to order 
EE 2 ,. The current reads 



3d 



cr D (^ 2 ) e + a dis E + a ln E 



(76) 



Here, the first term — describing the linear dark 
conductivity — follows directly from Eq. (|39p and (. . .) e 
denotes the energy averaging over the period u) c . The 
two other terms are microwave-induced corrections, at 
this order in the external fields completely independent 
of each other. The displacement contribution, propor- 
tional to cr dls , originates from the photon-assisted dis- 
placements in Eqs. (|65|) . (|66|) if one substitutes the equi- 
librium function f £ for f e . In the limit 2tt 2 T/uj c 1, 
a straightforward calculation yields 

<J dis = a D ^- [V- Re (£+£*_)] (Ki-H 3 ) , (77) 

TZi = udu (i> s P s+u ) e , (78) 
K 3 = {i> 2 e -v s v e+ j) e , (79) 
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where t^ 1 is expressed in terms of the partial contribu- 
tion (|7T|) to the disorder-induced scattering rate as fol- 
lows (the angle brackets denote averaging over 9): 

t- 1 = 2 (t 9 \1 - cosfl) 2 ), = 3T- 1 - At- 1 + r" 1 . (80) 

The inelastic contribution to jd, proportional to cr ln , ac- 
counts for the microwave-induced change 8f e = f e — fj 
of the distribution function. For E — > 0, the last term in 
Eq. ([68]) vanishes and e = e'. To first order in 



(81) 



The approximation in the last line is valid for T>u. 

At order EE 2 in jd, the correction (|5T|) of order E 2 
should be substituted into Eqs. (|65|) - (|67|) taken at first 
order in E, which gives 



a ln = -(7 D / de v 2 e d £ 5f e = a D -^TK 2 , (82) 



At' 



K 2 = U>d u (v 2 e (v e+u] + Ve-u)). 



(83) 



Below we summarize the most important properties of 
MIRO that follow from the above results. 



a. Period and phase. These are determined by the factors 
Tli defined in Eqs. dZHJ), ((79]), and (JHSl). In the limit of 
overlapping LLs, the DOS is given by Eq. (|3"Tj) . v(e) — 
1 — 2(5cos(27re/cj c ) with 5 — exp(— 7r/oj c r q ) <C 1, so that 



fti = K 2 /A 



-AS 



, TTU! 



2-kuj 



Tl 3 



a c2 • 2 ^ 

AS 2 sin — 



(84) 



We see that, for both the displacement and inelastic 
mechanisms, MIRO are proportional to sin(27ra;/w c ) with 
a negative coefficient in front of it in accord with the ex- 
perimental findings (for the displacement mechanism, the 
phase of MIRO also agrees with the observed one for large 
iruj/u) c , when 72-3 can be neglected compared to IZi). 

The phase of the oscillations can be qualitatively un- 
derstood as follows. Owing to the oscillatory behav- 
ior of v e v e ±u+eEAX , which enters Eq. (1551) via the fac- 
tor ([6"T]) , the average displacement of a cyclotron orbit 
AX (Fig. ITOb ') is positive (negative) for a small posi- 
tive (negative) detuning 8 U — lo/lu c — N from the iVth 
harmonic of the CR. The uphill drift for fe > pro- 
duces a contribution to the dissipative current (1651) which 
is directed against the electric field, i.e., er dls < 0. For 
the inelastic mechanism, a small positive detuning 5 U in- 
creases (decreases) the occupation of states lying right 
above (below) the maxima of v{e) compared to equilib- 
rium (Fig. [11]) , while the negative detuning leads to the 
opposite sign of the nonequlibirum correction to the dis- 
tribution function. The phase difference between the os- 
cillations of u(e) and those of d £ 8f e in Eq. (|82l) is such 
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FIG. 11 (a) Emission/absorption of microwave quanta for 
the detuning 5 U = ijj/oJc — 2 = 1/4 and (b) the DOS v(e) 
and the resulting oscillatory distribution function f £ for T = 
cj c , compared to the thermal distribution /J. Adapted from 
iDmitriev et al\ (|2005l '). 



that the sign of (—d e Sf e ) at the maxima of vie) coincides 
with the sign of <r iD (Fig.[TT]), i-e., S u > yields <r in < 0. 

In the limit of separated LLs, w c r q 3> 1, the 
DOS is a sequence of semicircles ([2"81 of width 2r = 
2(2w c Ar q ) 1 /2 « Wc7 i.e., j}( £ ) = Tq Re^ - (fe) 2 , 
where fe is the detuning from the center of the nearest 
LL. In this limit, one has 



fti=fto£$^sgn(fi n )«2(|n n |) 



(85) 



Alolo, 



K 2 = -ft 0T ^^sgn(O0$ 2 (|r! n |), (86) 



i-£fti(|n„| 



(87) 



where IZo [which also enters Eq. ([76]) ] reads 

K = (v 2 ) £ = 16w c /37r 2 r . 

The parameterless functions of fl n = (lj — nu c )/T are 
nonzero at 0< |fl n | <2, where they are given by 



Hx{x) 
H 2 (x) 



(2 + x) 
3x[(2 



' (A + x 2 )E(X) -AxK(X)} 
■x)E(X)-AK(X)}/8, 



(89) 
(90) 



Air$ 2 (x) = 3ir arccos(a; - 1) - x(l + x)y/x(2 - x) . (91) 

Here X = (2 — x) 2 / (2+x) 2 and the functions E and K are 
the complete ellip tic integrals of the first a nd second kind, 
respectively [see ([Dmitriev et all l2007bh for a graphical 
representation of the functions ([89 ]) -([91 ]) ]. A distinctive 
feature of the photoresponse in separated LLs is the pres- 
ence of windows in lo c within which \lo — nu) c \ > 2r for 
any integer n and hence both intra- and inter-LL single- 
photon transitions are impossible, i.e., v E i> £+ul = for 
any e. The resulting gaps in the dependence of the pho- 
toresponse on_^BJo£_a_grven radia tion frequency were ob- 
served by iDorozhkin et al. | (|2005h . see Fig.rjJ At suffi- 
ciently high microwave power, the photoresponse in the 
gaps becomes visible due to multiphoton effects (Sec. lVIl) . 
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d. Sensitivity to different types of disorder. The displace- 



FIG. 12 Magnetoresistivity at a radiation frequency of (a) 
40 GHz and (b) 121 GHz as a function of B, compared to the 
magnetoresistivity in the absence of radiation. The shaded 
boxes around the points at which (a) cj = cj c /2 and (b) 
id = 3lj c /2 mark the range of B within which the magne- 
tores istivity with and witho ut radiation coincide. Adapted 
from lDorozhkin et al\ (120051 ). 



b. Polarization dependence. While the phase of MIRO in 
cr dls and cr ln is essentially the same, their polarization 
and temperature dependences are qualitatively different. 
For linear polarization with s± = e^'^/v^ character- 
ized by the angle ip between the directions of the mi- 
crowave and dc fields, cr dls contains the anisotropic term 
Re (£+£*) oc 2Res+ = cos2?/>, while cr in oc V does not 
depend on tp. The po larization dependence of cr d ls is illus- 
trated in Fig. 13 of (|Vavilov and Aleinerl . [200l . Exper- 
imental results on the polarization dependence are dis- 
cussed in Sec. lVII.A.6l 



c. Screening of the microwave field. The dimensionless 
field £± and the polarization parameters s± in Eq. (|T5[) 
represent the strength and the polarization of the total 
(screened) microwave field. As discussed in Sec. III.C31 
(see also Sec. IVILA.6l) . in high mobility samples the ac- 
tive circular component of the incoming electromagnetic 
wave is strongly suppressed near the CR due to a nearly 
complete reflection. Expressed in terms of the field E, 
and the polarization of the incoming wave, Eq. 
acquires the form 

Jo), P (o) 



(0) 



uj(uj ± uj c + ij) 



(92) 



where 7 [Eq. (|57|> . see also footnote [TU] is assumed to be 



much larger than 



Since in typical experiments 7 



is of order cu > t _1 , the CR in V [Eq. <[7I|)] is strongly 
broadened by the screening. Note also that the polariza- 
tion of the screened field is generally different from that 
of the field in the incident wave and depends on B. 



ment contribution cr is hig hly sensitive to the de- 
tails of the disorder poten tial ( Auerbach and Pal l2007t 
K ho das and Vavilo vL l2008lh which are difficult to extract 
from standard transport measurements. In the mixed- 
disorder model (Sec. III. Al) . believed to provide an ad- 
equate description of high mobility 2DEGs, the partial 
contributions (J7TJ) to the scattering rate are given by 



1 



1 



Tsh T q , sm 1 + X n 



(93) 



where the first term describes scattering off the short- 
range component of the random potential (modeled here 
as white-noise disorder) and the second term describes 
small- angle scattering off the smooth random potential 
created by remote donors. The characteristic scattering 
angle in the latter case is y/x = (2kpd)^ 1 -C 1. Since 
in high mobility structures the quantum scattering rate 
Tq -1 , given by Tq" 1 , is much larger than the transport 
scattering rate t _1 = t^ 1 — rj~ , the quantum scattering 
rate is dominated by the contribution of the long-range 
component: t^ 1 3> . That is, the LL broadening is 
largely determined by the smooth disorder, whereas the 
relative weight of the short- and long-range components 
in 1/r ~ 1/T s h + x/ T q,sm may be arbitrary. Moreover, 
the displacement contributi on g dls is proportional to th e 



rate r" 1 given by Eq. ([80]) (|Khodas and Vavilovi [2008h . 



12 X 2 



(94) 



which contains one more power of the small parameter 
Y in front of compared to t -1 (jVavilov and Aleinen . 
2004). As a result, even a small amount of short-range 
scatterers may give the main contribution to <7 dls , as il- 
lustrated in Fig. [K3"h . 
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FIG. 13 (a) Dependence of t/t* on the relative weight r/r s h 
of the short-range component of disorder in the transport 
scattering rate for fixed r/r q = 100, 30, 10. (b)-(c) Depen- 
dence of cr in and cr dis on T for (b) T q = 10 K and (c) T q = 3 K. 
The curves for er m are calculated for Tj n = 2r at T — 1 K; the 
dashed curves for <r dls — for strong short-range disorder with 
t/2t* — 1, t/t s i, = 0.6, T* = 2K; the dash-dotted curves 
for 



t/t q = 
(2009b) 



for purely long-range disorder wit h t/2t» = 0.1, 
50, and ~ 6 K.) Adapted from iDmitriev et al\ 
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e. Temperature dependence. The displacement contribu- 
tion [Eq. (jTT|) ] is T-independent (apart from the T de- 
pendence of the DOS, which becomes relevant at high 
T for both mechanisms, see Sec. lIII.B.l.fj) . while the in- 
elastic contribution [Eq. (l82j) ]. proportional to the inelas- 
tic scattering time T ln , grows with lowering T. A de- 
tai led analysis of the ine lastic relaxation was performed 
bv lDmitriev et ali ([20051 ) . It was found that, for relevant 
T, the amplitude of the oscillations in the distribution 
function is controlled by electron-electron collisions. In 
overlapping LLs, the ansatz / = / T + f{e) d e fj with 
periodic (p(s) = ip(e + cu c ) reduces the linearized collision 
integral (l68|) to 



StinU} = - 



ir 2 T 2 + e 2 dfj 



2 de 
+ tp(e')-tp(e'-E)]) e , E , 



(A(E) [tp(e) -tp(e + E) 
(95) 



where the angular brackets denote averaging over e 1 and 
E within the period uj c . For a harmonic modulation of 
the distribution function with ip(e) oc cos(27re/w c + 9), 
the collision integral then acquires the form Sti n oc 
(A(E){1—cos(2ttE /uj c ) ]) E , where the last factor strongly 
suppresses the energy-relaxation rate compared to th e 
out-scattering rate. As shown bv lDmitriev et al. 
the relaxation time approximation is justified for this col- 
lision integral and the inelastic scattering time that en- 
ters Eq. (152")) is given in the relevant domain oj c <C T <C 
Wc(w c t) 1/2 by 



j 2ui Tcc '- e ' ife—u fe+ui) 



where 



l (e,T) 



7T 2 T 2 



■In 



2v F /a B 
w c (cj c t) 1 /2 



(96) 



(97) 



For T > a, Eq. $96$) gives r in ~ 0.822r ee (0, T) cx 
T~ 2 (Fig.[T3b,c), while for T < u the result is r in = 
(tt 2 T/2lu) t oo (0, T). It follows that the T~ 2 scaling of a in 
for T ^ w crosses over to a T _1 scaling for T <C 
In separated LLs, the relaxation time approximation is 
not accurate parametrically; still, for an estimate — up 
to a factor of order unity — one can use Eq. (|52"|) with 

Tin ~ (r/w c )T ee (0,T). 



by the damping factor 5 — systematically yield shorter 
^SdH) < T ^ ^ g di scussec j i n Sec. lII.C.2| the reason for 
the enhanced damping is that the SdH oscillations are 
sensitive to small, with an amplitude of the order of lu c , 
macroscopic inhomogeneities of the chemical potential, 
whereas the magnetooscillations that survive at T lo c 
are robust with respect to them. 

The experimentally observed MIRO (Sec. lIII.5| show 
at sufficiently high T an exponential suppr ession of the 
oscillation amplitude a s T is increased ()Hatke et ali 
l2009ct IWiedmann et ali l2010bl) . This effect can be ex- 
plained in terms of a T-dependent ren ormalization of 
the DOS by electron-electron interactions ( Chaplikl . 1 1 9 T it 
Dmi triev et al. |, l2009ri iRvzhii et ali [200J). Specifically, 



the quantum scattering rate r" 1 that enters the DOS 
v(e) — 1 — 25 cos(27r£ / ! lu c ) via the Dingle factor i5 should 
be substituted by 



f- 1 =r- 1 
q q 



_-i 



(98) 



where r cc 1 is given by Eq. (|97[) . Substitution of the mod- 
ified DOS into Eqs. ([r?rj|) - (|B5|) yields an additional factor 
exp(-T 2 /T q 2 ) in the amplitude of MIRO for T above the 
temperature T q defined by 27t/cj c t oc (0, T q ) = 1 [which 
gives T q ~ (wcEf) 1 / 2 up to a logarithmic factor], see 
Figs. [8] and [T3b.c. It is worth recalling that electron- 
electron scattering produces no additional exponential 
damping for the SdH oscillations, see footnote [5] 



g. Relative weight of the inelastic and displacement contribu- 
tions. The discussion in Sees. III.B.l.d-f shows that the 
relative weig ht of CT dis and cr in in the amplitude of MIRO 
strongly depends on T and the correlation properties of 
disorder. For an estimate of their relative importance, we 
omit the anisotropic part [Eq. (|77[) ] and the subleading 
for lj 3> uj c term proportional to 7?-3 [Eq. ([84]) ] in <r d 
and combine Eqs. ([77 |) -([79 ]l and (f52 ]> -<[gl j) into 



_dis 



0"D 



2t^ 
r 



,-dis 



Ad 2 V — sin . (99) 

2t* / UJ C L0 C 



Since r in - e F /T 2 fScc. lIII.B.l.c|) , Eq. ([99]) defines the 
temperature T* ~ (eft*) 1 / 2 Jt at which cr( dls ) = a m . In 
the case of smooth disorder [when r" 1 is given by the 
second term in Eq. ([9"4"])]. Eq. yields 



f. Exponential B- and T -damping. According to 
Eqs. ((Til), (j82jt- (|84| . the Dingle factor squared, 
5 2 = cxp(— 27r/o; c T q ), determines the damping of the 
MIRO in the limit of overlapping LLs for both the 
displacement and inelastic mechanisms. This same 
damping factor describes also the oscillations of the 
absorption coefficient for T 3> uj c (Sec. III.C3[) and, as 
will be seen in Sec. [V] HIRO and PIRO. Experimental 
data for these types of oscillations measured in the 
same sample support this prediction. On the other 
hand, the SdH measurements — if the data are fitted 



a™ _ 4r in (T)r* r in (T) 



r(dis) 



3t 



JL < lis 

Tsh T 



(100) 



For overlapping LLs, Eq. (|100j) gives T* ~ {ef/t^) 1 / 2 > 
T q ~ (wc^f) 1 ^ 2 , i.e., in the case of smooth disorder, 
the inelastic contribution dominates in the whole tem- 
perature range T < T q where MIRO can be observed 
(Fig. [T3b.cV By contrast, for the case of a strong short- 
range component of disorder (T sn ~ r ~ t*), one has 
T+ ~ Tqfa/uJcT 2 ) 1 / 2 , much smaller than T q for uj c t ^> 1. 
As a result, the cr ln -dominated T~ 2 dependence of the 
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MIRO amplitude for T < T+ crosses over as T is in- 
creased to a plateau (whose height is given by er dls ) in 
the interval T* < T < T q (Fig. [13b) before starting to 
fall off exponentially for T > T q . In the intermediate 
case of T* ~ T q , the range of T in which MIRO are T- 
independent shrinks to zero (Fig. Il3h). 

The experimental results, presented in Fig. [51 are in 
overall agreement with the above theory. Apart from the 
T dependence, the two mechanisms are different in the 
photoresponse at the "odd nodes" uj/oj c = n + 1/2, where 
cr in = while cr( dis ) oc TZ 3 f 0, see Eqs. ([77|l. (|75]l. and 
([54"]). Moreover, their polarization dependence is different 
(Sec. lIII.B.l.bj) . Experimental studies of these features 
[omitted in the estimate (|99|) ] can serve as an additional 
tool to quantify the relative magnitude of a dls and er ln . 



2. Saturation of the inelastic contribution at high radiation 
power and/or in a strong dc field 

When the inelastic contribution dominates the oscillatory 
photoresponse (Sec. |III.B.l.g[ ), the leading correction to 
<7 ln linear in V and independent of E [Eq. (1521 ] is given by 
the terms of higher order in V and E 2 in the distribution 
function. These are governed by the kinetic equation 
([55]). a generalized — to include higher powers of V and 
E 2 — s olution of which in th e regime of overlapping LLs 
reads ( Dmitriev et all 12005) 



of = o— -j— sin T r in , V, C) , 

ZTT OS UJ r 



27TU) 



r/T in +Vsm 2 ^ 
where the dimensionless parameter 

C = 7re dc = 2ireER c /uj c 



+ C 2 /2 



(101) 
(102) 

(103) 



characterizes the strength of the dc field. The current 
induced by the inelastic mechanism is then given by 



U = a u E [ 1 + 2S 2 - 2S 2 _F(r in , V, Q } 



(104) 



The dependence of the resulting magnetoresistivity on V 
and E is illustrated in Fig. [14] 

Let us consider first the limit E — > 0. At fixed u>/oj c , 
MIRO saturate with increasing microwave power: 

„ s 47TW 7TLJ „ 9 7TCJ T 

•F(T in ,"P,0) -)• cot — , 7? sin 2 — >— . (105) 

U! c U! c U! c Tin 

In the limit of high power V T/r m , the maxima and 
minima of MIRO at e ac = e ac [Eq. (|63[) ] shift towards the 
nearest nodes at e ac = n, with the amplitude of MIRO 
at the extrema being propotional to (rinP/r) 1 / 2 : 



1 



T in V 



1/2 



n 



T2tt 



T in V 



1/2 



(106) 



Note that the T dependence changes in the nonlinear 
regime. As follows from Eq. (|102[) , there remains — for 
arbitrarily strong V — a range of e ac around the nodes 
within which the photoresponse is linear in V (apart 
from heating by microwaves which can modify Tj n ); how- 
ever, this range shrinks with increasing V . In a broader 
vicinity of uj/uj c — n which includes the nearest max- 
ima and minima of MIRO, the result (|104[) applies — 



under the conditions r q <C oj, 



< T 



< 



Til, 



M p, which are 

met in typical experiments dDmitriev et all l2007al 120051 : 
lHatke et all l2011al iKhodas and VavilovL l2008h — for ar- 
bitrary microwave power and for arbitrary relation be- 
tween the components of mixed disorder ([93]) . The re- 
sulting behavior with V reproduces the observations in 
Fig. 13 and is also in agreement with othe r experimental 
results (jMani et al 1 12010b lYe et al. I. l200lh . 

For T sn ~ r, i.e., in the presence of a strong short- 
range component of disorder, one should for r- m < t 
take into account the displacement contribution and 
also multiphoton processes (the latter are important at 
V sin 2 nuj/uj c > 1, see Sec. lVI.B.4[) . This results in 
Eq. (1 1 061) with a f a ctor of order unity substituted fo r 
Tin /V (|Hatke et all l2011at IKhodas and Vavilovl 120081 ). 
If only the smooth component of disorder is present, 
Eqs. (|104|) . (|102|) are applicable and give the main con- 
tribution to jd in a broader range V sin 2 ttoj/uj c < t/t^; 
at higher V, both the excitation of higher angular and 
temporal harmonics of the distribution function and the 
multiphoton processes become important, which results 
in the emergence of a series of distinct strongly nonequi- 
librium regimes as ui/u i r . is varied [for more details see 
( Dmitriev et all l2007aT) ]. The higher harmonics of the 



distribution function also determine the Hall photore- 
sponse (Sec. HrfB.4T) . 

The above consideration of the nonlinear-in-P effects 
Eqs. (|102p . (|106p shows that the linear response conduc- 
tivity [Eq. (I104[) at E — > 0] can become negative de- 
spite 5 <C 1, see Fig. [14] As discussed in Sec. IIV.Bl this 
causes an electric instability leading to the formation of 
domains. The electric field E* in the domains can be 
obtained from the current-voltage characteristics (I104[) . 
as illustrated in the inset to Fig. [HI see Sec. HV.Cl for 
details. The behavior of the magnetooscillations at high 
microwave power and in a strong dc field is further dis- 
cussed in Sees. IVllVIl 



3. General approach: Quantum kinetic equation 

A comprehensive approach to noncquilibrium trans- 
port in high LLs was initially formu l ated for the 
case of smooth disorder ( Dmitriev et all l2007al 120051 : 
IVavilov and Aleinerl . l2004f) and l ater generalized to the 
case of generic [Eq. (I7T1)1 disorder (fbmitriev et al. ., 2009b; 
K ho das and Vavilo v, 2008; V avilov et aLl . l2007l) . The key 
steps of the derivation are highlighted below. 
(1) Transformation r — > r — £(<) to a moving coordinate 
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FIG. 14 Resistivity normalized to the dark Drude value 
(given, as we assume lj c t 3> 1, by jd/or>E) as a function 
of uj c /ui at fixed cjr q = 2n for different levels of microwave 
power (r m /T)V\u, c =o = 0.24 (bottom), 0.8, 2.4 (top). Inset: 
I—V characteristics at the minima (marked by the dots) of the 
u) /u>c dependence of the resistivity. The arrows show the value 
of the dc field E* in domains i n unit s of Eo = nuj c T /2r; n e_R c . 
Adapted from lDmitriev et all (|2005l 1 . 



frame, where obeys 



with i xy = —e 



dt 



Of 



dc 



(107) 



= 1, is similar to the transfor- 



mation to Floquet state s, used in a nu m ber of 
works (lAuerbach and Pal 120071: iKashubal l2006b| 



Lee and Leinaasl . 
20061: 



Park, 2004: 



2004; 



Lva pilin and Patrakov, 
Torres and Kunoldl. 12003 



Volkov and Takhtamirovl 120071) on the displacement 



mechanism of MIRO. This transformation eliminates 
homogeneous electric fields in the Hamiltonian of a 
clean electron system at the expense that the impurity 
potential, "dressed by the electric fields", becomes 
time-dependent. Scattering off the moving impurities 
can then change the energy of electrons. 
(2) Keldysh equations within the SCBA. In the moving 
frame, the Green's functions G a and the self-energies 
S Q (a — R,A,K refers to the retarded, advanced, 
and Keldysh components) are related within the SCBA 
(Sec. III.C.ll) as 



^21 — 



al 2 q 



-tg£ai 



iqr qol £) — iqr 



(108) 



where W q is the correlation function of the bare disorder 
potential, the subscript (21) denotes times t% and t\ on 
the Keldysh contour, and £21 = £(fo) — £(^i)- 
(3) Quasiclassical approximation. For a degenerate 
2DEG in high LLs, Eq. f|108[) can be reduced to a simpler 
quasiclassical equation in the "action-angle" representa- 



tion: 

EM = -i£ 2 i , (109) 

g^(ip) = ioj c ^2G^(h + k; <p) , (110) 

k 

where the operators n and <p are canonically conjugated, 
[n,ip] = —i. The eigenvalues of h and <p are the LL index 
and the angle coordinate of the momentum, respectively. 
The effects of disorder and external fields are encoded in 
the integral operator 1C, 



dip' e ifeF ^ riv ~ rl v')€ 21 



2tt 



(111) 



The distribution function /, denned by G R -G A - G K = 
2(G R f — fG A ), commutes with <p. Accordingly, the oper- 
ator ip, which enters Eq. (|109[) and the impurity collision 
integral 



iSt im {f} = ± R f-ft A + (t 



K 



E R )/2, (112) 



can be treated as a c-number. The result is the quantum 
kinetic equation 



h uj c 8 v ) /21 
St; m {/}21 



St in {/} 2 1 = St im {/} 



(113) 



dU 



£21 (523/31 - haaii) 



hlK-23923 + /23^3l531 



(114) 



Here t = (ti+t^)/^ is the "center-of-mass" time. The in- 
elastic collision integral St; n {/}2i accounts for electron- 
electron scattering and for the coupling to a thermal 
(phonon) bath. 

Except for specific cases (Sec. lVI.A.3l) in the regime of 
separated LLs, the effect of external fields on the spec- 
trum can be neglected, meaning g R - t2 m ^ ne moving 
frame remains the same as in the static frame and does 
not depend on t and <j). It follows that the direct current 



j = 2ev F Jdeu(e){n ip f{e,ip 1 t))-2evQe F dti(t), (115) 

is determined by the first angular harmonic of the 
Wigner-transformed distribution function /(e, ip, t). The 
bar denotes time averaging over the period of the mi- 
crowave field in the steady state. 

The approach formulated above [Eqs. (|111[) . (|113[) 
(|115p ] validates the description of the inelastic and dis- 
placement mechanisms of MIRO in Sec. ini.B~Tl It has 
been applied to a wealth of nonequilibrium phenomena 
that we address in the remaining part of the review. 



4. Quadrupole and photovoltaic contributions to the 
photoconductivity 



Care ful study of the quantu m kinetic equation (|113[) 
shows ( Dmitriev et all l2007al) that there are in total 
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four different contributions to the photocurrent at the 
leading order EE 2 . Diagrams (A)-(D) in Fig. [T~| repre- 
sent four ways to obtain perturbatively the first angu- 
lar harmonic /io of the distribution function f(e, <p, t) = 
Y,f^n eluv+mult , which defines the current (TTIo| . start- 
ing from the isotropic dark distribution fj. In the dis- 
placement contribution (A), /io oc (cj c d ¥ ,) _1 Stim{.f~} re- 
sults directly from the action of the collision operator 
Sti m oc EE 2 . In the inelastic contribution (B), /io oc 
(w c 5 ¥ ,) _1 St im {(5/ 00 }, where St im oc E and Sf 00 oc r in E 2 
is the microwave-induced correction to the isotropic part 
of /. Apart from <5/oo, the action of Stj m oc E 2 on 
/~ results in excitation of the second angular harmon- 
ics /2o- The linear dc response in the resulting state, 
/io oc (uj c d v )~ 1 St im {f2o} with St; m oc E, produces 
the "quadrupole" contribution (C) to the direct cur- 
rent. Finally, in the "photovoltaic" mechanism (D), a 
combined action of the microwave and dc fields Stj m oc 
EEu produces nonzero temporal harmonics /io and /i2- 
The ac response in the resulting state, fio oc (d t + 
w c 9 v )~ 1 St im {/io,/i2} with Stim oc E u , also gives rise 
to the direct current. 
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FIG. 15 Diagrams representing four distinct quantum contri- 
buti ons (A)-(D) to the phot ocurrent at order EE^. Adapted 
from lDmitriev et al\ i|2007ah . 



Calculations performed in ( Dmitriev et al . 2007af ) 

(for smooth disorder) and ( Dmitriev et all 2009bl ) (for 



generic disorder) show that all four mechanisms (A)-(D) 
produce microwave-induced corrections to the direct cur- 
rent j — Jdark = VphE. The full photoconductivity tensor 
E 2 consists of components of different symmetry 



d s + d a h s + h a 
h s h a d s - d a 



(116) 



oscillating with oj/uj c . The Hall part is solely due to 
mechanisms (C) and (D). In particular, the quadrupole 
mechanism (C) yields an unusual symmetric off-diagonal 
term h s violating the Onsager symmetry, while the anti- 
symmetric h a is produced by the photovoltaic mechanism 
(D). The photovoltaic mechanism also yields corrections 
to both isotropic (d s ) and anisotropic (d a ) diagonal com- 
ponents. It is worth mentioning that, unlike contribu- 
tions (A)-(C), the photovoltaic contribution (D) cannot 
be described in terms of the scattering rates similar to 
(|66[) . In particular, it involves not only the DOS given 



by the imaginary part of but also quantum correc- 
tions to the real part of the Green's function which do 
not appear in the golden-rule approximation. 

In the case of smooth disorder, the largest contribu- 
tion comes from the inelastic mechanism (B), which con- 
tributes only to d s . In this case, d s 3> d a ~ h s ~ h a , 
meaning the contributions of subleading mechanisms (A), 
(C), and (D) are comparable in magnitude. A suffi- 
ciently large amount of short-range scatterers enhances 
contribution (A), see Sec. lIII.B.l.dl while the magni- 
tude of contributions (C) and (D) remains essentially 
the same, d s ~ d a ^> h s ~ h a . Thus, the de- 
scription of MIRO in Scc. IIII.B~Tl which accounts only 
for the inelastic and displacement contributions, is jus- 
tified as long as the diagonal part of cj( ph ) is con- 
cerned. Mechanisms (C) and (D) yield oscillations in 
the Hall component (in the terms h„ 



which 

were also observed experimentally dMani et all 2004a; 
IStudenikin et all 12004 ' iWiedmann et aU l2011bh . While 
the amplitu de of the oscillations in th e Hall component 
observed by IWiedmann et al. ( 2011b ) is in good agree- 



ment with theory, iMani et al 

(2004) reported considerably stronger oscillations 



( 2004a[ ); lstudemkin et all 
In 



fact, a similar issue is known for SdH oscillations whose 
measured amplitude in the Hall component of the resis- 
tivity often considerably exceeds the theoretical predic- 
tions. It is likely that the physical mechanism of the 
enhancement of the Hall component is the same in both 
cases. 



5. Classical mechanisms of MIRO 

Most experimental results on MIRO reported so far 
have been consistently explained within the quantum ki- 
netic approach de scribed above. These quantum effects 
are d irectly linked (jDmitriev et al. L l2003tlFedorvch" et all 
|2010D to quantum magnetooscillations (w) in the dy- 
namic conductivity fSec. III.C~~) . On the other hand, 
the quasiclassical memory effects, discussed in Sec. III.B~fl 
also produce large u) /w c -oscillations in the ac conductiv- 
ity, <j^(uj) [Eqs. ~23D, d~~~D], see the comparison of the 
classical and quantum contributions in Fig. ~0 In turn, 
the oscillatory ct( c '(ll>) translates into w/o; c -oscillations 
of the electronic temperature T e , while the oscillating T e 
manifests itself in dc transport via the renormalization 
of the elastic scattering rate by electron-electron interac- 
tions. The resulting oscillato ry contribution to the pho- 
toconductivity has the form (jDmitriev et al. 1. 120041) 



(0) 



e— ph 



-v 



CD 



oj 2 A<j^(uj) 
e F T cr D (w) 



(117) 



where r° n ph is the electron-phonon inelastic relaxation 
time (which controls the heating of the 2DEG) and Aa^ 
is given by Eqs. (f~3[) and (j24[) . The phase of the mag- 
netooscillations in A(7pk(0) is opposite to the phase in 
the dynamical conductivity Aa^ c \to) and is shifted by 
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7r/4 with respect to the quantum oscillations [Eqs. ([82]) 
and ([77)) ]. At low T, the magnitude of the oscillations in 
(|117p may become comparable to that in Eq. ([82]) since 
the ratio Act( c )(w)/<7d(w) can be of order unity (Fig. [5]), 
while the small parameter uj/sf ~ 10 _1 — 10~ 2 is easily 
compensated by the large ratio t?~ p jr\ n of the electron- 
phonon and electron-electron inelastic relaxation times 
[electron-electron scattering provides for relaxation of the 
fast oscillations in the energy distribution of electrons in 
Eq. I[82p but conserves the total energy of the electron 
gas]. However, unlike the quantum oscillations of the 
photoconductivity, the oscillations described by Eq. (|117[) 

saturate with increasing V well before Acr^ (0)/<7d be- 
comes of order unity, and therefore, cannot cause ZRS. 
Namely, in the regime of strong heating (T e — T > T), 

A<7pk (0)/<7d does not exceed T e jep 1. 

It is important to mention that all classical corrections 
to the Drude formula induced by the homogeneous exter- 
nal ac and dc fields vanish to zero in the Boltzmann equa- 
tion framework in the case of a parabolic dispers ion rela- 
tion and er-independent r ([Dmitriev et all 120041) . Apart 
from the interaction-induced effects discussed above, fi- 
nite nonlinear corrections may result either from the 
energy dependence of r [which typically occurs on the 
scale of £f, thus yielding a n additional small fa ctor 
T/ep compared to Eq. (|117p ([Dmitriev et al ... 2004)] or 
from a weak nonparabolicity [which leads to the ap- 
pearance of an even larger energy s cale of the order 
of the bandgap in the den ominator ([Joas et al I 12004 
iKoulakov and Raikhl . I2003T )] . By contrast, the Landau 
quantization provides the much smaller scale w c in the 
dependence on e [Eqs. ([78)) and ([83]) ]. which explains the 
leading role of the quantum mechanisms of MIRO. 



Two recent works addressed classical ef- 
fects close to the boun dary of a 2DEG. 

IChepelianskii and Shepelvanskvl ( 20091 ) studied the 
influence of the microwave field on the phase space 
portrait of electron dynamics near the sample edge. 
The conclusion was that at uj/uj c = n + 1/4 the mi- 
crowave radiation tends to trap electrons in trajectories 
propagating near the edge, while at lo/uj c — n the 
particles a r e mor e efficiently kicked out into the bulk. 
iMikhailovl ([201 if ) discussed an effective electrostatic 
(pondcromotive) potential for electrons (whose cy- 
clotron dynamics is modeled using Newton's equation 
with a friction term), created by a s trong inhomogeneity 
( Mikhailov and Savostianoval l2007t) of the microwave 
field in the near-contact region. Depending on the ratio 
Ld/uJc, the ponderomotive force repels electrons from 
or attracts them to the near-contact area. Since the 
above works did not develop a systematic theoretical 
analysis of the contributions of the boundary effects 
to the resistivity of the 2D system, it is difficult to 
compare these contributions to those predicted by the 
bulk theory fSec. HII.B.l[) . In any case, in the systems 
in which MIRO are observed transport is dominated by 
the bulk contribution, so that even a strong modification 



of edge transport may modify the total current only 
weakly. Indeed, experimentally, no dependence of the 
MIRO amplitude on the sample dimension or geometry, 
characteristic of the edge effects, has been reported 
so far. Moreover, edge transport obviously plays no 
role in Corbino geometry. Furthermore, recent obser- 
vations of MIRO usi n g various contactles s techniques 
(|Andreev et all 1201 it iBvkov et all l2010bl) have ruled 
out the contact-related phenomena as a generic cause of 
MIRO. 



C. Related microwave-induced phenomena 



1. Microwave-induced photovoltaic effects 



Re cent experiments ( Bvkovl l2008eJ: iDorozhkin et all 
2009). performed using an asymmetric contact con- 
figuration, discovered that microwaves can induce os- 
cillatory current and voltage signals very similar to 
MIRO even in the absence of external dc driving, see 
Fig. [TBI These phot ovoltaic effects were attributed 



(Dmitriev et al ., 2009a) ;o the violation of the Einstein 
relation between the conductivity and the diffusion coef- 
ficient, induced by the microwave illumination. Consider 
an inhomogeneous system characterized by the electro- 
static potential </>(r) and the electron density n e (r), both 
smoothly varying in space. In equilibrium, the electron 
system is characterized by the constant electrochemical 
potential r\ = e0(r) + n e (r)/x^ dark ' , while a generic weak 



r (dark) 



V0- 



k)^(dark)\7^ At 27T 2 T/W C > 1, 



perturbation leads to the current j( dark ) 

e£) (dark) Vne = _ e £)(d 

the dark compressibility ;v/ dark ) = 2fo, so that the Ein- 
stein relation in equilibrium reads rj^ dark - ) = 2e 2 ^o-D^ dark - ) 
where £>( dark ) = (D 2 (e)) £ R 2 c /2T, see Eq. ([76]). 
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FIG. 16 (a) Photovoltaic current and (b) voltage at a fixed 
radiation frequency / = 50, 112 and 168 GHz versus u> c /^ in 
the absence of dc bias applied t o the sample. The curv es are 
offset for clarity. Adapted from IDorozhkin et al\ l)2009t ). 



Analysis of Eqs. (|55]) - (I551) , adapted to the in homoge- 
neous conditions, shows ([Dmitriev et al. , 2009a) that the 
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direct current under microwave illumination, 



j = — <TrN<t> ~ ev§DXJr\ 



(118) 



T dis 



necessarily contains an extra "anomalous term" 
which violates the Einstein law. The term er^ = < 
includes the inelastic contribution to MIRO (|82|) and 
the most important part cr dls cx 1Z± of the displace- 
ment contribution ([77|) . The other part rj dls oc IZ3 enters 
the nonequilibrium diffusion coefficient, D = £)( dark ) + 
tJ3 ls /2e 2 fo. The total conductivity a, which defines the 
direct current j = aE in a homogeneous system, is given 
by (j = a. q + 2e 2 uoD. 

The current (I118jl is nonzero even at Vr/ = 0, provided 
there is a nonzero built-in electric field (created, e.g. by 
an asymmet ric contact configuration or by local gates) , as 
observed in dBvkovl l2008al l2010t iDorozhkin etall . I2009L 



l2011bl : IWillett et al 



2004). The simplest model, pro- 



posed bv lDmitriev et al. ( 2009al ). assumes that the built- 
in field U c /L is provided by two contacts with the differ- 
ence of work functions U c , attached to a 2DEG stripe of 
width L. The I-V characteristic for this setup reads 



j = <7,,U C /L + aV/L 



(119) 



where the voltage between the contacts eV = rj\ x -o — 
v\x=l- As seen from Eq. (| 1 1 9[> . there is a finite oscillatory 
photocurrent j at zero bias voltage V = 0, as well as a 
finite photovoltage V (asymmetric with respect to V = 0) 
at j = 0, as shown in Fig.[TBJ The corre sponding field and 
densit y distributions are calculated in ( Dorozhkin et all . 
l2011ah . 



2. Microwave-induced compressibility oscillations 

The violation of the Einstein relation under microwave 
illumination has a number of other important conse- 
quences. In particular, the generalized compressibility 
Xq = n g /(—e(f)q), defined as a static density response 
n q e lqr to a weak local perturbation <\> q e %qr of the screened 
electrostatic potential, strongly deviates from its equi- 
librium value ^( dark ) = 2vq. Assuming a stationary 
nonequilibrium state with the homogeneous component 
of the current j — 0, one obtains from Eq. (|118[) : 



a v /e 2 D = a/e 2 D 



Note that, in contrast, the quantity Xg^ = n ql l l q is n °t 

modified: Xq — 2fo- Using the generalized compress- 
ibility, one can define the nonequlibrium screening length 



(120) 



A = 



eD 
2na 



(121) 



which replaces the equilibrium Thomas-Fermi screening 
length Ao = e/27re 2 x < - dark ' ) in all electrostatic problems. 

At q- 1 > l in = {Dt^) 1 / 2 but qR c < 1, the local ap- 
proximation for the current (|118l) is still justified. At 



the same time, the effect of spatial variations of the 
elect ric field on the distri bution function becomes essen- 
tial (IVavilov et aL ^ 2004j ). Summarizing the re sults of 
Dmi triev et al. I (l2009al) and IVavilov et all t0Q4) for the 
leading inelastic contribution, the generalized compress- 
ibility at qR c <gC 1 is given by 



(0) 
Xq 

2v 



= l-6 2 T(T in ,T,() 



1 + q 2 lf 



(122) 



meaning the amplitude of the oscillations in Xq reduces 
by a factor of 2 at ql in ^> 1 compared to the limit q — >• 0. 
These results suggest that ZRS corresponds to a plateau 
in the compressibility, as illustrated in Fig. [17] for the 
case qli n 3* 1- Therefore, local measurements of the com- 
pressibility [e.g., usi ng the techniques that utilize single- 
electron transistors (|Ilani et <zZ.L l200Ct 1200 lh ] may provide 
a real space snapshot of the domain structure in ZRS. Ex- 
perimental work in this direction is currently underway. 




FIG. 17 Microwave-induced oscillations of the compressibility 
as a function of Ld c /uj at fixed cjr q = 2n and (n n /T)'P\ U j c= o = 
1. In ZRS, the electric field inside domains E* fixes the com- 
pressibility at the level shown by a dashed line. Inset: the 
dependence of the compressibility o n the electric f ield E in- 
side the domain wall. Adapted from lVavilov et a/1 l|2004T ). 



3. Effects of a parallel magnetic field 



lYang et all (2006) observed the effect of a parallel 
magnetic field B\\ on MIRO: the oscillations were sup- 
pressed as B|| is increased and virtually disappeared at 
B|| ~ 0.5 T. This observation came as a surprise be- 
cause, according to the theory presented above, MIRO is 
an orbital effect and should not be essentially sensitive 
to the spin degree of freedom. Moreover, the transverse 
motion of 2D electrons (that woul d be a ffected by the 
parallel field) is frozen. In (|Manil . |2005| ). MIRO were 
indeed practically independent of B n up to about 1 T , 
i.e. in the same range of B» as in [Yang et all (|2006l ). 



The most likely explanation of the experimental data 
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of lYang et~al. is the effect of Bu on the quan- 

tum scattering rate l/r q . Indeed, the parallel magnetic 
field does affect the structure of the wave function across 
the 2DEG plane, which can in general modify the impu- 
rity scattering rate. This is a nonuniversal effect, as it 
depends on microscopic details of the nanostructure in 
which the 2DEG is formed, as well as on the character 
of imperfectio n s. It is plausible that in the QW used 
bv lYang et all (2006) this effect may be stronger than in 
a typical single-interface structure. I n fact, a closer in- 
spection of the experimental data of lYang et al. I (12006b 
shows that the dark resistivity (and thus the transport 
scattering rate 1 /r) was enhanced by the parallel field by 
a factor of about 3. While in general the enhancement 
factors for 1/r and l/r q may be different (because 1/t 
and 1/Vq may be controlled by different components of 
disorder), this clearly shows that Bu did strongly enhance 
disorder-i nduced scattering. T he effects of B\\ were also 
studied bv lHatke et all (|201 leT ) in the nonlinear (with re- 
spect to the dc field) transport regime, see Sec. IV.AI It 
was confirmed that Bu suppresses quantum oscillations 
via the enhancement of l/r q . 



4. MIRO in multi-subband structures 

The microwave-induced effects in magnetotrans- 
port have also been intensively studied in struc- 
tures wit h several occupied subbands, both experi- 
mentally (iBvkov et all 120081 l2010dl; iGusev et all , l201lt 



IWiedmann et al 



. l2010bUd I2011U 120081 12009^ an d the- 
oreticallv (jRaichevl I2008T IWiedmann am2010bll2008l 
l2009d) . Already the fir st experiments in a dou ble Q W in 



the Hall-bar geometry ( Wiedmann et aZ.l . l2008f) as well as 
in Va n der Pauw and Corbino geometries ( Bvkov et all 
12008b demonstrated strong interference between MIRO 
and MISO (the latter are discussed in Sec. III.C.4|) . 
Specifically, it was observed that (i) MIRO are strongly 
enhanced (suppressed) in the maxima (minima) of MISO 
and that (ii) MISO in the minima of MIRO are in- 
verted by sufficiently strong microwave radiation. Sim- 
i lar effects were recently observed also in triple QWs 
(IWiedmann et aUl2009cj). 



IWiedmann et all ( 20081) performed a direct compari- 
son of the experimental results wit h theory that combines 



the inelastic mechanism of MIRO ( Dmitriev et all 20051) 



and the intersubband scattering effects ( Raichevl . |2008|) 



The theory reproduced with a good accuracy the compli- 
cated interference pattern, as well as the nontrivial de- 
pendence of the combined oscillations on the microwave 
power, frequency, and t emperature, observed in the ex- 
perim ent. In all cases (IWiedmann et all . l2010bl I2008L 
2009c]), the compa rison confirmed t he T 2 -dependence of 
MI RO amplitude dDmitriev et all 120051 ) up to T ~ 4 K. 
In (jMamani et all 20081 ). the same T~ 2 -dependence of 
the inelastic scattering time reproduced the exponen- 
tial T-damping of MISO, governed by the interaction- 
induced broadening of LLs [the high-T damping of 



MIRO (|Dmitriev efaH . l2009li Irlatke et all, l2009cl )— of 
the same origin — is discussed in Sec. lIII.B.l.ej . 

In the case of a balanced double QW, with negli- 
gible interlayer correlations of the scattering potential 
and identical scattering rates, the theoretical result of 
IRaichevI (|2008l) reads [cf. Eq. (fTMj) and Eqs. (|58> (|qTJ1)] 



3d 
<J D E 



2vrA 



, (123) 



where T is given by Eq. (|102[) . This expression clearly 
shows features (i) and (ii) mentioned above. Indeed, 
MISO are inverted for T > 1. Further, if 25 2 {J 7 - 1) ex- 
ceeds unity, the conductivity (and, therefore, resistivity) 
at the MISO maximum becomes negative, which leads to 
the emergence of ZRS (Sec. IIVI) . ZR S of th i s kind have 
been recently observed by IWiedmann et all ( 2010d ). see 

Fig. mi 




FIG. 18 Measured resistance (solid lines) of a two-subband 
electron system without (no mw) and under microwave irra- 
diation at / = 143 GHz as a function of B. The inverted 
by microwaves MISO peak at B ~ 0.27 T exhibits vanishing 
resistance. Dashed line: theory. Inset: symmetric (S) and 
antisymmetric (AS) wavefunctions for two lo west subbands 
in a wide QW. From I Wiedmann et all (|2010d ). 



5. MIRO in spatially modulated systems 



The discussion in Sees. IIILC.ll and IIII.C. 21 shows that 
the violation of the Einstein relation strongly mod- 
ifies the transport and thermodynamic properties of 
the microwave-illuminated 2DEG already in the pres- 
ence of a weak smooth inhomogeneity. The oppo- 
site limit of a stron g sho r t-peri o d spatial modula- 
tion was studied in dDietell. 12 006: Die tel et all. 1200. 
Joas et al. . 20051: Kennett et all 120051 : iRobinson etal. 
2004 iTorres and Kunoldl 120061) 



Dietell (|2006l ): iDietel et all |2005l ): IJoas et all (|2005l ) 

examined the transport properties of the 2DEG with 
a unidirectional static modulation [this system has 
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been intensively studied before in the context of 
Weis s oscillations and t r anspo r t anisotropies (iBeenakkei 



19891: 



1989 : iMirlin and W dlflc. 1998; 



_ .Gerhardts et al... 

Weiss et all Il989t IZhang and Gerhardtsl Il990t) ~ In the 



presence of a periodic potential V(x) — V cos Qx the de- 
generacy of LLs is lifted and LL bands appear with the 
dispersion relation 



lk = cu c (n + 1/2) + VJo(QR c ) cos(QfcZl) 



(124) 



where the Bessel function Jo is an approximation that 
is valid for high LLs. The DOS of the clean modu- 
lated system has square-root singularities at the band 
edges. In the x-direction, the photoconductivity can be 
calculated using the golden-rule approach formulated in 
Sec. IIII.B.ll The results for the displacement and in- 
elastic mechanisms are given by Eqs. ([77]) and ([82} with 
the oscillating factors (1751) . (f79")l . and (|83l) calculated us- 
ing the DOS corresponding to Eq. (jl24|) . The oscillating 
factors acquire additional sign changes, related to singu- 
lar DOS, which can be detected in the experiment. In 
the perpendicular direction, the linear (in the dc field) 
photoresponse diverges for the displacement mechanism. 
The divergence is cut off by inelastic processes. As a re- 
sult, both the displacement and inelastic contributions to 
the transverse photoresponse are proportional to Tj n and 

have a similar magnitude. 

iKennett et al. | (I2005T) : iRobinson et~al\ (|2004h stud- 
ied the magnetoresistivity in the presence of surface 
acoustic waves (SAW ) . This perturbation combines 
( Levinson et ali . Il998l ) both the unidirectional spatial 



modulation with the wavevector Q and the ac excita- 
tion at frequency w = sQ, where s is the sound ve- 
locity. Apart from an ac analogue of the Weiss oscil- 
lations, which results in an anisot ropic positive MR os - 
cillating with both QR C and u)/u> c , IKennett et al. I (j2005ft : 
IRobinson et all ( 2004f ) calculated the associated quan- 
tum contribution for the inelastic mechanism of MIRO, 
which potentially gives rise to SAW- induced ZRS. The 
quantum contribution is given by Eqs. (|104|) and (|102l) . 
where the microwave power V is replaced by properly 
normalized SAW power absorbed by the 2DEG. Since 
the absorbed SAW power is proportional to the dynamic 
conductivity of the classical system at the wavevector Q, 
ct^q, which oscillates with QR C , p xx (Fig. [19]) shows ad- 
ditional Q-Rc-oscillations on top of relatively slow uj /lu c - 
oscillations identical to MIRO. 

lYuan et al\ (20065) investigated experimentally the mi- 
crowave photoresistance in the presence of a triangular 
antidot superlattice. This study revealed the conven- 
tional MIRO with a superimposed magnetoplasmon peak 
on top of a series of narrow geometrical resonances posi- 
tioned at 2R C = 7 n a, where a is the lattice period and "f n 
are numbers of order unity corresponding to various com- 
mensurate orbits. Essentially no interference between 
MIRO and the commensurability oscillations was ob- 
served. Note, however, that a strong interplay of the two 
types of oscillations was predicted ( Torres and Kunoldl . 
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FIG. 19 Calculated resistivity magnetooscillations in the 
presence of a surface acoustic wave for lot = 7 and vf/s = 63 
(solid line) and ojt — 60 and vf/s = 84 (dashed), both c urves 
for Lur q = 6 and wr in = 302. From IKennett et al\ (|2005f ). 



2006) for a short-period superlattice with a of the order 
of l B < Ra- 



ti. Microwave-induced _B-periodic oscillations 



iKukushkin et al I (I2004D discovered a different kind of 
photoresistance and photovoltage magnetooscillations in 
the region oj < uj c . These i?-periodic oscillations were 
attributed to the interference of edge magnetoplasmons 
(EMP) emitted from different potential probes along the 
Hall-bar edge. The condition for constructive interfer- 
ence of the EMP injected from two contacts separated 
by the distance L reads qL — usL/v — 2irN, where the 
EMP velocity v is proportional to the Hall conductiv- 
ity &xy oc n s /B. This interpretation e xplains the period 
AB oc n s /u)L found in the experiment. IKukushkin et aJ\ 
(2004) established the following properties of the B- 
periodic oscillations: (i) the amplitude of the resistivity 
oscillations shows a threshold behavior as a function of V\ 
(ii) the photovoltage scales linearly with V in the limit of 
small V and saturates with increasing V; (iii) the satura- 
tion power for the photovoltage is close to the threshold 
power for the resistivity oscillations; (iv) the resistivity 
and photovoltage oscillations were phase shifted by ir/4 
with respect to each other; (v) the threshold V for the 
resistivity oscillations is an order of magnitude lower for 
the microwave electric field perpendicular to the Hall- 
bar edge (which supports the EMP scenario). The phe- 
nomenon showed only a weak dependence on T in the 
broad interval T ~ 1 — 10 K. The photovoltage reduced 
only by one order of magnitude at T as high as 70 K. 
Systematic theory of the -B-periodic oscillations remains 
t o be developed. 

IKukushkin et ~al. I (l2005h also reported a successful op- 
eration of the microwave spectrometer based on the B- 
periodic photovoltaic oscillations. The B-tunable selec- 
tive detection was demonstrated in the frequency range 
20 GHz < / < 150 GHz and for temperatures up to 
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FIG. 20 B-periodic magnetoresistance (solid line) and pho- 
tovoltage (dashed) oscillations at / = 40 GHz measured in a 
Hall bar with potentia l probes separated by 0.5 mm. From 
iKukushkin et all l|2004 ). 



T ~ 80 K (with a foreseen possibility of extension to 
the THz frequency range). Implementation of plasmon 
resonances in such a device allows one to overcome the 
temperature limitation T < ui of the conventional selec- 
tive detectors based on electronic transitions. The ma- 
jor advantage of magnetoplasmons is that they can be 
exploited at arbitrary ujt provided w c r ^> 1, while the 
operation of the conventional plasmonic B = 0-detectors 
in field-effect transistors is limited to the region ujt > 1. 

In a recent experiment IStone et all ( 20071 ). where 
MIRO, ZRS, and B-periodic oscillations were observed 
simultaneously, the period AB of the B-pcriodic os- 
cillations was reported to scale as AB oc n s /ui and 
be independent of the distance L between the con- 
tacts. This re s ult i s in contradiction to that by 
IKukushkin et all (2004): the controversy has remained 
unresolved. Several oth e r experiments, in particular by 
iDorozhkin et all (|2007al) : lYuan et all (|2006f) . reported no 
interference of the magnetoplasmon effects and MIRO 
with each other. This agrees with the calculation by 
IVolkov and Takhtamirovl (|2007h . according to which the 
magnetoplasmons can strongly affect MIRO only in the 
limit of well separated LLs — which is out of the range 
explored in the above experiments. 



the way to zero forming the zero-resistance states (ZRS), 
see Fig. 1211 With appropriate microwave intensity, tem- 
perature, and sample quality, ZRS can span magnetic 
field ranges corresponding to several tens in filling fac- 
tors. However, unlike the QH effect, vanishing of di- 
agonal resistance in microwave-irradiated 2DEG is not 
accompanied by Hall quantization. 
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FIG. 21 Longitudinal [Ru(B), left axis] and Hall [Rh(B), 
right axis] magnetoresistance under microwave irradiation. 
Longitudinal magnetoresistance R(B) witho u t irra diation is 
also shown. Adap ted from (a) [Mani et all (|2002D and (b) 
IZudov et all (|2003l ). Parameters: (a) microwave frequency 
/ = 103.5 GHz, temperature T = 1.3 K, electron density 
n e ~ 3 x 10 11 cm" 2 , and mobility fi ~ 1.5 X 10 7 cm 2 /Vs; 



IV. RADIATION-INDUCED ZERO-RESISTANCE STATES 
(ZRS) 



(b) / = 57 GHz, T z 
H ~ 2.5 x 10 7 cm 2 /Vs. 



1.0 K, 



3.5 x 10 1 



and 



A. ZRS: Experimental discovery and basic properties 

Zero resistance is a rare occurrence in condensed mat- 
ter physics, usually signaling a novel s tate of matter, such 
as superconducti v ity an d QH effects (jvon Klitzing et all , 
119801 : iTsui et all 1 19821) . see Sec.lIV.B.2l f or few further 
examples. Experim ents by iMani et all ( 2002t ) and by 
IZudov et al. I (120031) on microwave-irradiated very-high 
mobility (// > 10 7 cm 2 /V-s) 2D EG rev ealed that the 
lower order minima of MIRO (Sec. lIII.A")) can extend all 



Over the past decade ZRS were observ ed over a wide 
frequ ency range, from as low as 9 GHz (IWillett et all . 



l2004h to as high as 254 GHz (|Smet et all . 120051 ). Ex 



pcriments in ultra-high mobility (/x > 2 • 10 cm /V-s) 
2DEG have also revealed ZRS stemming from the min- 
ima of fractional MIRO (Sec . lVI.ATT) . namely near e ac = 
3/2, 1/2 and 2/3 (|Zudov et all l2006b| ). as well as a ZRS 
presumably appearing due to a fre quency mixing unde r 
bichromatic microwave irradiation ( Zudov et all 2006b). 
Very recently, ZRS were also observed to emerge from the 
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MISO maxima in double QWs (|Wiedmann et aLl . l2010dh 
see Fig. [IS] Shortly after the discovery of ZRS, ex- 
periments bv lYang et all ( 20031 ) in Corbino-disk shaped 
2DEG demonstrated the existence of corresponding zero- 
conductance states ( ZCS) (Fig.l22p. Later studies by 
Bykov et all (|2010bl ) reported observation of ZCS in 
capacitively-coupled 2DEG. As discussed in Sec. lVII.B.U 
ZCS were also realized in a no n-degenerate 2D electron 
syste m on liquid 3 He surface ( Konstantinov and Konol 
120101) . 
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FIG. 22 Magnetoconductance G XX (B) under microwave irra- 
diation of frequency / = 57 GHz. These data were obtained 
at T = 0.65 K on a Corbino shaped 2DEG with inner (outer) 
diameter 0.5 (3.0) mm, density n e ~ 3.55 x 10 11 cm" 2 , and 
mobil ity n ~ 1.28 x 10 7 cm 2 /Vs. Adapted from lYang et all 
(120031 ) . 



the nonlinear resistance is zero (lAndreev et al. L 120031: 
lAuerbach et al. 1. 120051: iFinkler and Hafperinl . l2009h . To 
date direct confirmation of absolute negative resistivity 
and the domain model has been limited. First exper- 
imental support of the domain picture was provided by 
IWillett et al. I (120041) who observed large voltages between 
an internal and an extern al contact to 2DEG in the ZRS 
regime. Experiments by IZudov et al. | (|2006al) employ- 
ing radiation of two distinct frequencies have found that, 
away from ZRS, the bichromatic photoresistance can be 
approximated by a simple superposition of monochro- 
matic photoresistances, e.g. 8p u:iu:2 ~ a.\8p u:i + ct28pu 2 - 
Theoretica lly t he bi ch romatic pho t orespo nse was ad- 
dressed in iLeil (120061) : E ei and Liul ()2006al ). Assuming 
that such a superposition also holds in the regime where 
one of the frequencies gives rise to ZRS, e.g. underly- 
ing microscopic resistance is negative p + 8p u>1 < 0, but 
p + 8p UlU>2 > 0, one can easily reconstruct negative resis- 
tance which otherwise is masked by instability (Fig. |2~4"]). 
In other words, the experiment uses one of the frequencies 
to probe the absolute negative resistance corresponding 
to a ZRS induced by another frequency. 



Very recently, iDorozhkin et al. I (l2011bl ) have studied 
temporal dynamics of Hall voltages between internal 
probes and observed random telegraph signals in the ZRS 
regime. These signals were interpreted in terms of spon- 
taneous switching between two nearly degenerate config- 
urations of spontaneously formed domains (Fig. l2"5)l . 



Temperature depend ence of the ZRS has been exam- 
ined by several groups (Mani et all . 120021 : IWillett et all . 
12004 IZudov et al. , 2003). All studies found that ZRS dis- 
appear with increasing temperature (see Fig.l2"3"]). trans- 
forming into MIRO minima where the resistance approx- 
imately followed the Arrhenius law p u cx exp(— A/T). 
The "energy gaps" A cx B extracted in this way exceed 
relevant T ~ u> by an order of magnitude. 



G 




FIG. 23 Temperature evolution of ZRS and MIRO obtained 
under / = 85 GHz ill umination on the same sample as in 
Fig.[2TJi. Adapted from lMani efaH (|2002l V 



As discussed in Sec. IIV.Bl a homogeneous state with 
negative resistivity is electrically unstable and is expected 
to break into domains with local current density at which 
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FIG. 24 Magnetoresistances Ri(B) at /i = 31 GHz (a), 
R 2 (B) at f 2 = 47 GHz (b), and R 1+2 at both /i and f 2 (c). 
Maximum-maximum, minimum-minimum and maximum- 
minimum overlaps are marked by f, J,, and £, respectively. 
The upper boundary of the shaded area represents the aver- 
age of monochromatic resistances Ri and R 2 . Dotted lines 
in (a) and (b) represents reconstructed negative resistance. 
These data were obtained on a 2DEG (~ 5 mm x 5 mm) 
with density n e — 3.6 x 10 11 cm" 2 and mobi lity p ~ 2.0 x 10 7 
cm 2 /Vs. Adapted from lZudov et all (|2006aT ). 
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FIG. 25 Photovoltages (for B = 95 mT, T = 0.5 K, and / = 
48.1 GHz) across different pairs of adjacent contacts shown in 
inset (middle column). The voltages at t > (no radiation) 
correspond to the zero reference level. Closed loops in insets 
indicate the Hall current flo w between the contacts i n two 
different configurations. From iDorozhkin et al\ (|2011bT h 



B. Phenomenological approach: instabilities of 
non-equilibrium 2DEG in strong magnetic field and the 
emergence of domains 

As discussed in Sec. IIII1 for a sufficiently strong ra- 
diation power (and other parameters in the appropriate 
range), the linear resistivity of the microwave-illuminated 
2DEG becomes negative. Such a state is unstable, which 
results in formation of sponta neous current and v oltage 
breaking rotational symmetry ( Andreev et all [20031) ; the 
corresponding theory is presented in Sec. IIV.B.1I In 
Sec. IIV.B.2I we will review some previously obtained re- 
sults on related systems. 



1. ZRS effective theory: Spontaneous symmetry breaking 



The starting point of the theory of lAndreev et al\ 
(2003) is a model with a constant Hall resistivity p H 
and a non-linear current dependence of the dissipative 
component E^ of the electric field: 



E = PH j xz + ±-E d (\j\). 



(125) 



It is assumed that the function E d (j) has a negative 
derivative at j = and a zero at j = j c . The theory 
is phenomenological, i.e. it assumes that there is cer- 
tain microscopic mechanism that provides such a current- 



voltage characteristics. A distinct feature of the consid- 
ered strong magnetic field regime is a duality between 
current and voltage. Specifically, Eq. (|125j) is appropriate 
in the Hall-bar geometry, when only one (say, j x ) compo- 
nent of the current is present; the two terms in Eq. (| 125[) 
then give the Hall (E y ) and the longitudinal (E x ) field 
respectively. In contrast, for the Corbino-disc geometry 
one field component is zero (say, E x = 0) and one is in- 
terested in j(E y ). Since the Hall component is the dom- 
inant one in the conductivity tensor, the dependence of 
the field dependence of the dissipative component jd(E) 
is obtained from E d (j) in Eq. (|125D by exchanging cur- 



rent and voltage (jBergeret et all . 1200 



E 

j = a H zxE+— Jd (\E\) 



(126) 



with an — {ph)~ 1 and jd{E) = ajjEd((THE). In other 
words, the dissipative current-voltage characteristics is 
S-shaped in the first case and N-shaped in the second 
case, see Fig. [211 




FIG. 26 S-shaped [(a), Hall bar setup, Eq. (Tl25|) ] and re- 
shaped [(b), Corbino setup, Eq. (|126p ] current-voltage char- 
acteristics. 



Equation (|125|) is supplemented by the continuity 
equation, 



Further. 



dn e /dt + V • j = 0. 



E = -V0(r), 



(127) 



(128) 



where (f>(r) is the electrostatic potential; its variations are 
related to those of density according to 



S(j)(r) = / d 2 r'U{v,v')8n e (v'). 



(129) 



where L/(r, r') is the Coulomb interaction (that may be 
screened in the presence of an external gate) . To explore 
the stability of the system, one considers a small devi- 
ation in density Sn e (r, t). Using Eqs. (|126p - (|129[) and 
linearizing in 5n e , one gets 



dSn e /dt — Vod Vf7(5n e , 



(130) 



where ad is the dissipative differential conductivity 
(o- d ) a p = d{jd) a /dE l3 satisfying {a d ) a p = {<J d )p a . Note 
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that the Hall component of the conductivity dropped 
out of Eq. (|130p because of its antisymmetric character. 
The stability requires that the real part of all eigenval- 
ues of the operator on the r.h.s. of Eq. f|130[) is non- 
positive. Transforming it into the momentum space and 
using U{q) > 0, the stability condition reduces to the 
requirement that eigenvalues of ad are positive, which 
yields 



djd/dE > 0; 
Jd/E > 0. 



(131) 
(132) 



In other words, the stability requires that both absolute 
and differential conductivity are non-negative. As is clear 
from the above derivation, Eqs. (|131[) and (|132|) represent 
the stability to longitudinal and transverse fluctuations 
of the electric field with respect to its direction in the 
state under consideration. 

Therefore, when the linear-response conductivity is 
negative the system is unstable near j = E = 0: fluc- 
tuations will grow until the stability point with current 
j c and electric field E c — pnjc is reached, see Fig.l2"rJl 
This implies that the system will break into domains 
with spontaneous currents and fields of these magnitudes. 
The simplest possible domain structure consists of two 
domains as shown in Fig. [27] for the case of Hall-bar 
and Corbino-disc setups (see Sec. IIV.DI for a discussion 
of more complicated domain structures). The values of 
spontaneous current and field correspond to dissipative 
current jd = (i.e. zero absolute dissipative conduc- 
tivity) in the Corbino disc setup and to dissipative field 
Ed = (and thus zero absolute dissipative resistivity) 
in the Hall bar geometry. When a finite current (for the 
Hall bar) or voltage (Corbino disc) is applied, the domain 
wall shifts to accommodate it. The system remains in the 
ZRS (or zero-conductance state in the Corbino geometry) 
until the applied current (field) exceeds j c (respectively, 

E r ). 



(|2004l ) who included in the r.h.s. of Eq. (fT26)) the dif- 
fusive term — e_DVn e and assumed that, contrary to 
the conductivity, the diffusion tensor is not affected by 
microwav e radiation. [This assu mption was later jus- 
tified by iDmitriev et al. I l|2009alb ] This introduces in 
the problem the nonequilibrium screening length that 
governs the width of doma in walls, see Sec. IIII.CT2l 
Volkov and Pavlovskiil fj2004h assume however a 3D re- 
lation between the potential and the density [i.e. they 
assume V(q) ~ 1/<Z 2 momentum dependence of the 
Coulomb interaction rath er than V(q) ~ 1/q app ropri- 
ate for 2D geometry, see ( Dorozhkin et al\ , l2011al) ] and, 
as a consequence, obtain the 3D screening length = 
(eD /A-K\a d SD "'\) 1 /' 2 . The resulting equations are used to 
find the spatial profile of the domain wall. The central 
result of IVolkov and Pavlovskiil (|2004l ) is the residual re- 
sistivity of a finite sample in the ZRS phase (with a single 
domain wall assumed), which is found to be negative and 
exponentially small: p rcs cx — exp{— const L/X^jj}. Here 
L is the sample size in the direction transverse to the 
domain wall. 

iDorozhkin et~al. I (l2011aft established the electrical sta- 
bility condition in a finite sample with realistic 2D 
Coulomb interaction. Instead of conditions (|131l) and 
(I132[) (which in the linear regime read ad > 0), one gets 
ad > —eD/2L, or A -1 > tt/L in terms of the 2D nonequi- 
librium screening length (|121[) . In the diffusion-stabilized 
regime of negative conductivity , —eD/2L < ad < 0, the- 
ory of IDorozhkin et al\ ( 2011al ) predicts the emergence 
of two regions with opposite directions of electric field 
(this requires broken inversion symmetry, for instance, 
due to asymmetric contact configuration). The ampli- 
tude of these fields increases when the system approaches 
the instability threshold. This effect is a precursor of the 
domain structure in the regime ad < —eD/2L. 



2. Earlier results on related problems 





FIG. 27 Simplest domain structures for (a) Hall-bar and (b) 
Corbino-disc geometries. The current density j c and the elec- 
tric field E c — pnjc in the domains are fixed by the condition 
(a) E d = or (b) j d = 0, see Fig.^ Corbino-disc [Hall-bar] 
system adapts to an applied bias V = E c (Li — L2) [applied 
current I — j c (Li — L2)] remaining in zero-conductance [zero- 
resistance] state. 

A furt her step forwar d in th e study of the ZRS was 
done by iBergeret et all (|2003l) ; IVolkov and Pavlovskiil 



The instability of a semiconductor with dc nega- 
tive c onductivity was first pointed out by IZakharovl 
(1960). Subsequently many authors performed de- 
tailed analysis of such instabilities in various sys- 
tems with negative differential or absolute conductiv- 
ity. The corresponding literature is very extensive, and 
we restrict ourselves to quoting several early articles 
dBonch-Bruevich and Koganl 1965 ; Elesin and Manvkinl . 



119671: iRidlev 

review ( Volkov and Kogan . 



19631 Volkov and K oganl 



19671) . the 

th e books 

dBonch-Bruevich et all \l97& iPozheii Il98ll: IScholll . 



196% , a nd 



l200lh . The reader can find there an overview of vari- 
ous mechanisms of emergence of negative conductivity 
(characterized by N-shape or S-shape current-voltage- 
characteristics), derivation of stability conditions, and an 
analysis of domain formation as a result of instabilities. 
In this contex t, particularly well known is the Gunn effect 
(|Gunnl . Il963f) : an N-type current- voltage-characteristics 
leads to an instability resulting in formation of moving 
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domains and in microwave generation. The correspond- 
ing semiconductor device (Gunn diode) has found various 
application in high-frequency electronics. 

The abso l ute ne gative conductivity was observed by 
Banis et al. I jl972ft [see also Sec. VI.3 of the book by 



Pozhelal "( 198lh ] ■ This was achieved by exposing GaAs 



samples exhibiting a negative differential conductance 
to a microwave r a diation. The explanation given by 
iBanis et al\ (fl972l) : iPozhelal (Il98lh assumes that the mi- 
crowave field can be treated adiabatically and simply 
leads to an (oscillating in time) shift of the operation 
point of the device. After the time averaging of the corre- 
sponding differential conductivity this leads to a negative 
dc conductivity. It seems, however, that this explana- 
tion (that would be perfectly correct for a low-frequency 
field) does not fully catch the physics of the experiment. 
Indeed, the frequency of the microwave field (of order 
10 GHz) does not seem to be small compared to inverse 
characteristic times in the device, so that the adiabaticity 
assumption is not met. Furthermore, the authors report 
that, in the absence of an external circuit, a large sponta- 
neous dc voltage developed on the sample (whose polarity 
could be stabilized by a small symmetry breaking pertur- 
bation) . The amplitude of the voltage was determined by 
the zero of the current-voltage characteristics, I(V) = 0. 
This spontaneous formation of dc field domains resulting 
from the microwave-induced absolute negative conduc- 
ti vity indicates certa in analogy between the experiment 
of IBanis et al\ ( 1972h and the systems constituting the 
subject of this review. 

It was experim entally discovered ( Basun et all 1 19831 : 
iLiao et all I1980D that intense laser illumination gener- 
ates strong electric fields (detected via the splitting of 
luminescence lines) in ruby crystals. A pheno menological 
theory of this phenomenon was develo ped by IDvakonovl 
(|1984T ); IDvakonov and Furmanl (|l984f) . It is based on 
an assumption that the absolute conductivity in ruby 
becomes negative for a sufficiently strong illumination 
power. The resulting instability leads to formation of 
electric field domains, thus explaining the experimental 
observations. There is a clear analogy between the theory 
of ZRS as resulting from regions of negative resistivity in 
MIRO and this theory of ruby domains. There are, how- 
ever, also essential differences. First, this is the presence 
of strong magnetic field in the MIRO problem. Second, 
contrary to an isotropic 2DEG, ruby is anisotropic: the 
electric field and the current are directed along the C3 
axis. This modifies the analysis of the stability (only 
longitudinal fluctuations need to be considered) and ren- 
ders the part of the current- voltage-characteristics with 
I/V < but dl/dV > stable. To our knowledge, no 
microscopic theory of negative conductivity in ruby has 
been developed. 

The impact of a strong transverse magnetic field on 
this class of phenomena has bee n also appreciated ap- 
proxim ately half a century ago. iKazarinov and Skobovl 
(119631) showed that a quantizing transverse magnetic 
field may lead under non-equilibrium conditions to strong 



non-linearities in current-voltage-characteristics, in par- 
ticular, to a negative differential resistance. The fact that 
in the case of large Hall angle (w c r ^ 1) current-voltage- 
charactertics have a dual shape in the Hall bar (j y = 0) 
and Corbino disk (E v = 0) geometries (in particular, 
an S-shaped characteristics of the Hall bare device gives 
rise to an N-shape d characteri s tics for the Corbino disk ) 
was po i nted o ut bv lBassI (|l965l) ; iBogomolov et all (|l967l ) . 
iKoganl (|!968h studied an instability in a system with neg- 
ative diffe rential r esista n ce at intermediate Hall a ngles 
(bj c T - 1). lElesinl (1 1969h : iGladun and Ryzhiil dl970T) and 
particularly iRvzhiiT i 19701) : Rvzhii et al. (|l98*6l ) discussed 
the photoconductivity of a 2D gas in a quantizing mag- 
netic field and strong electric field (see Sec. MI. B) and 
reached a conclusion of possibility to get an absolute neg- 
ative conductivity in this class of systems. 

Another broad class of semiconductor systems in which 
similar phenomena have been intensively studied are 
single and multiple QW heterostructures. In particle, 
a great deal of work was devoted to resonant trans- 
mission via double-barrier structure s, where a a nega- 
tive differential resista nce was found (jSollner et aUI 19831 
iTsu and Esakl Il973l) . More recently, it was predicted 
that under laser illumination the system may show 
absol ute negative resistance (|Dakhnovskii and Metiul . 
1995). Further, the absolute negative con d uctiv- 
ity was theoretically predict ed ( Ignatov et all \l99& 
iPavlovich and Epshteint Il976l) and experimentally ob- 
served ( Keav et aLll995T ) in semiconductor super lattices 
under THz radiation; another theo retical work consider s 
a superlattice in a magnetic field ( Cannon et aLl . l200(il) . 
The theoretical predictions include generation of sponta- 
neous voltages and currents resulting from instabilities. 
A negative a bsolute conductiv i ty was also found in a re- 
lated model ( Hartmann et all Il997l ) of a particle in a 
periodic lattice driven by an ac field and coupled to a 
dissipative bath. 

Absolute negative conductance under strongly non- 
equilibrium conditions emerges also in other typ es of 
systems. In particular, lAronov and Spivakl ( 19751 ) pre- 
dicted this effect to happen in a Josephson junction of 
two superconductors with different gaps one of which is 
subjected to light illuminat ion. This prediction was ex- 
perim entally confirmed in iGershenzon and Faleil ( 19861 
1988). In view of the spatially local character of the 
junction, the negative conductance does not lead in this 
case to an instability (contrary to extended systems that 
tend to break into domains in such a situations). Fi- 
nally, an absolutely negative mobility was obtained in a 
purely classical model of a Brownian particle s ubjected 
to non-equilibrium noise (jEichhorn et al. 

1, mil- 
it is clear from the above discussion that many of the 
aspects of the ZRS problem have earlier appeared in re- 
lated problems. On the other hand, to the best of our 
knowledge, the combination of the key features of the 
ZRS problem, namely, (i) strong magnetic field (large 
Hall angle), (ii) instability induced by an absolute nega- 
tive dissipative conductivity, and (iii) 2D isotropy of the 
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problem, has not appeared in any other context. 



C. Microscopic theory; determination of currents and fields 
in domains 

We combine now the phenomenological ZRS theory of 
Sec. IIV.B.1I with microscopic calculations of photoresis- 
tivity. As was discussed in Sec. the linear resistivity 
in the presence of microwaves may become become neg- 
ative around its minima. Specifically, let us consider the 
inelastic mechanism that is dominant for sufficiently low 
temperatures and assume first the regime of overlapping 
LLs. According to Sec. IIII.B.21 the linear resistivity is 
negative when the dimensionless microwave power sat- 
isfies V > V* > 0, w ith the threshold value given by 
(iDmitriev et all 12004 l2005h 



V = — 4<r — sin sin 2 — 

UJr. UJr 



(133) 



see Fig. [TU The spontaneous field in the ZRS do 
is found to be 



Er — En 



V 



En = 



ireR r V 2t; 



mains 



(134) 



For the marked MIRO minima in Fig. Qj] the correspon- 
dent values E* / Eq are pointed by arrows in the inset. 

In the regime of separated LLs with width 2r = 
2(2w c /7rr q ) 1 / 2 and assuming a small-angle impurity scat- 
tering, the threshold power is 



T 2 t/ujuj c t- u 



and the spontaneous field is of the order of 



1/2 



eRr. 



(135) 



(136) 



Within the displacement mechanism (relevant for 
higher temperatures) and for a small-angle impurity 
scattering the corre s pond ing results were obtained by 
IVavilov and Aleinerl ( 2004[) . In this case the character- 
istic value of the spontaneous field is given by Eq. (|136[) 
in both limits of overlapping and separated LLs. 

It should be pointed out that the form (r/rq) 1 / 2 of 
the factor in front of uj c /eR c in Eq. (|136|) depends on 
the character of disorder. For a model of mixed dis- 
order (1§3"|) and for overlapping LLs this factor is re- 
duced and becomes of order unity when the weight of 
the short-range component becomes sufficiently large , 
see (jDmitriev eTaH l2009bl: iKhodas and Vavilovl I2008D 
and Sec. IIII.Bl A re lated result was obtained by 



lAuerbach and Pail ( 20071 ). who studied the displacement 
mechanism for mixed disorder (with short-range compo- 
nent determining the transport rate) in the regime of 
separated LLs and found ££ c ~ T/eR c (if one sets the 
correlation length of short-range disorder to be ~ Xf)- 



To our knowledge, the inelastic mechanism has not been 
studied for the case of separated LLs and mixed disorder. 

These predictions can be confronted with the exper- 
iment. Focussing on the inelastic mechanism and on 
the regime of overlapping LLs and assuming character- 
istic values of parameters uj/2tt ~ 50 GHz, T ~ 1 K 
and rr 1 ~ 10 mK, r" 1 ~ 0.3 K, r = 10 mK and 



v F = 2 ■ 10 7 cm/s, one gets (|Dmitriev etaH |2004 [2005h 
the threshold microwave intensity required for for the 
emergence of the ZRS P* < 1 mW/cm 2 , in conformity 
with the experiments. Further, for these parameter val- 
ues and for V — V* ~ V* (i.e. for the microwave power 
exceeding its threshold by roughly factor of two) the esti- 
mated dc electric field in the domains, Eq. (I134p. is found 
to be E c ~ 1 V/cm. In the experiment of Willett et all 
( 20041 ) the voltage drop between an internal and an ex- 
ternal contact (separated by 200 /^m) generated by the 
radiation in the absence of the drive current was of the 
order of 5 mV for ui c /2tt ~ 20 GHz. Assuming the sim- 
plest domain geometry (i.e no additional domain walls 
between the contacts), this yields E c ~ 0.25 V/cm. On 
the other hand, the above theoretical estimate yields 
E c ~ 0.15 V/cm [taking into account the w 2 dependence 
of E c following from Eq. (|134l) ]. so that the experimental 
value is somewhat larger . In a more recent experiment 
( Dorozhkin et all l2011bl ) spontaneous fields of order of 
E c ~ 0.15 V/cm were reported for uj c /2tt ~ 50 GHz, 
which is in this case several times smaller than the the- 
oretical estimate. In general, the agreement between the 
theory and the experiment appears to be very reason- 
able if one takes into account some deviations of param- 
eters (relaxation rates, ratio of the microwave power to 
its threshold value) from those used in the theoretical 
estimate. 



D. Effective theories of the phase transition into ZRS and 
of ZRS dynamics 

The transition into Z RS belongs to the class of dynam - 
ical phase transitions ([Hohenberg and Halperi n|, 119771) . 
These phenomena are in general governed by non- 
linear differential equations (of hydrodynamic type) with 
stochastic (noise) terms and encompass spontaneous pat- 
tern formation in a great variety of systems driv en away 
from equilibrium (jCross and Hohenberg, 1993). Con- 
trary to conventional (thermodynamic) phase transi- 
tions, such phenomena in general are not characterized 
by a free energy functional. This invalidates, in partic- 
ular, the Mermin- Wagner theorem, opening the way to 
spontaneous breaking of a continuous symmetry in 2D 
systems, which is of direct relevance to the ZRS prob- 
lem. 

lAuerbach et al. l (|2005l) : lFinkier" et al. | (|2006l ) performed 
an analysis of possible domain patterns in the ZRS phase. 
In spirit of Sec. IIV.B1 they started with a model with a 
constant Hall conductivity cth and a non-linear filed de- 
pendence of the dissipative current jd, Eqs. (|126[) - (|129l) . 
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Their crucial observation is that for this problem one can 
define a Lyapunov functional 



G[<f>] = J d 2 rg(E(r)), g(E) = J dE'j d (E'), 



(137) 



with the key property dG/dt < 0. This implies that 
minima of the Lyapunov functional are stable steady 
states. To explore the ZRS phase (negativ e linear dissipa- 
tive co nductivity), lAuerbach et al\ ( 20051 ) : IFinkler et all 
( 2006f ) expand g(E) near the zero E c of the current- 
voltage characteristics jd(E): 



5 (E) = g(E c 



-(E-E C ) 2 + K \V-V\ 



(138) 



The last term here [penalizing large field gradients and 
determining the width of domain walls, l^w ~ (k/ctc) 1 / 2 ] 
is added on the phenomenological ground. On the mi- 
croscopic level, a finite width of the domain wall emerges 
when one takes into account the diffusive contribution 
— el)Vrt e to the current j. However, the local form of 
the corresponding term in Eq. (jl 38[) correspond to the 
relation between E and n e characteristic for 3D rather 
than for 2D systems, see the discussion in the end of 
Sec. IIV.B.11 It remains to be seen, to what extent this 
may affect the results. 

Minimization of the Lyapunov functional yields the 
simple domain structures in the Corbino and Hall bar 
geometries, see Figs. [27] and l2"5k . Introducing further a 
weak and smoot h disor d er field Edi a (r) , with [ Enis I <C E c , 
lAuerbach efal\ ([20051 ): IFinkler etal\ (|2006l ) generalize 
Eq. (|138p [denoted below as ffo(E)] to 



g(E, E dis ) = 5o (E) - tri (£)EE dis + Q(E 2 i; 



(139) 



The term proportional to Edi s induces a correction to 
the stability condition, with the result that the current 
density at a domain wall is j = — (7i(i? c )Edis + 0{E^- is ). 
For a short correlation length of the disorder a com- 
parison of the domain wall and bulk c ontribution to 
the L yapunov functional in the spir it of llmrv and Mai 
(|l975l) shows (lAuerbach et all 120051) ;hat no new do- 
mains is formed. On the other hand, disorder with a 
sufficiently large correlation length £dis favors the break- 
down of the syste m in multiple domains with a size set 
by £d i s - Further, lAuerbach et al\ (|2005l ): IFinkler et a!\ 
(2006) used this approach to estimate the effect of the 
potential (/>di s (r) (regular or random) that induces mul- 
tiple domain walls on the ZRS conductivity. They ob- 
tained a — Co\ (\Eidis\) / E! c , with a numerical prefactor 
C ~ 1 depending on the specific form (or statistics) of 

t he potential. 

IFinkler and Halperinl (|2009h found that even in a ho- 
mogeneous system creating additional domains may be 
favorable. They pointed out that the Lyapunov "en- 
ergy" e dw of the domain wall depends on the angle 
4> between the wall and the field in the domains. If 
£dw(7r/4)/edw(7r/2) < l/2\/2, additional domains as 
shown in Fig. [25b will emerge. The dependence e dw (</>) 





FIG. 28 (a) Proposed domain structure near metallic contacts 
(shaded grey regions) in a rectangular sample; (b) possible 
structure with an additional pair of domains. Arrows show 
the d irection of the electr ic field E c , see also Fig.|27[ Ad apted 
from IFinkler et aH (120061 ) ; IFinkler and Halperinl (|2009T ) . 



is determined by the precise form of the Lyapunov func- 
tion g( E). Remarkably, for the sim plest choice of this 
function IFinkler and Halperinl ( 20091) obtain for the ratio 
£dw(7r/4)/edw(7r/2) exactly the value l/2-\/2, which would 
imply degeneracy between the states with and without 
additional domains with 7r/4 domain walls. They con- 
jecture that a more realistic Lyapunov function might 
yield a sm aller value of this ratio , thus f avoring additional 
domains. IFinkler and Halperinl (|2009l) consider further 
the effect of spatial variation of the Hall conductivity 
an and show that in this situation non-stationary (peri- 
odic) solutions may emerge. The simplest example is the 
Corbino geometry with Hall conductivity changing lin- 
early along the y axis (radial coordinate), dan/dy = a. 
Then V • jjj = aE x . Combining this with jd = (char- 
acteristic for a homogeneous system away from domain 
walls), we get dn/dt = — V • j = — aE x (r). In the sim- 
plest domain geometry (two domains with electric field 
in the radial direction, Fig. Wf\ E x — and we do not 
get any time dependence. However, in the presence of 
additional 7r/4 domains, Fig. [28b, the field there is along 
the x axis, which impl i es tha t such domains should move. 
IFinkler and Halperinl ( 2009f ) performed also a numerical 
modelling of the problem in the torus geometry and with 
a spatially varying Hall conductivity and indeed found 
non-stationary, time-periodic solutions, in a certain range 
of the parameter a controlling a typical gradient of <jh- 

What is the character of the transition into the ZRS 
phase? For the case of equilibrium phase transitions the 
answer to such a question is usually based on the Landau 
(mean-field) theory of phase transitions complemented 
by a renormalization group (RG) analysis. A counterpart 
of this approach applicable to dynamical critical phe- 
nomena (including those in strongly non-equilibrium sys- 
tems) is known as the dynamical renormalization gr oup 
(jForster et all Il977t [Hohenberg and HalpeTml Il977| ). A 
powerful framework for its technica l implementation 
is the Martin-Siggia - Rose formalism (iDominicisl. 1 19761 : 
Dominicis and Pelitil , 1 19781: iJanssenl . Il976t iMartin et all 
19731) closely related to the Keldysh formalism. This 



approach allows one to cast the evolution governed by 
a stochastic non-linear equation into a Lagrangian form 
convenient for the implementation of the RG procedure. 

lAlicea et ali ( 2005! ) made a first step in application of 
these ideas to the problem of ZRS transition. They have 
started with formulating a general equation describing 
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the long-scale, long-time dynamics of a microwave-driven 
system: 



u\i\ 2 j 



^d 2 a 

-jj,V(t> + w c z x j - r j 
+?7iV 2 j + 7/ 2 V(V • j) - %V 4 j - ry 4 V 3 (V ■ j) 
-Mi ■ V)j - f 2 V(j 2 ) - !/ 3 (V- j)j 

+ 7 20j + 7 30z x j + C + • • ■ • (140) 



Here the /x and w c terms represent the electric-field and 
Lorentz forces, the r and u terms describe the non-linear 
resistivity, the term with <9 2 (on the left-hand side) origi- 
nates from the frequency dispersion of the resistivity, the 
r\i terms characterize its momentum dispersion, the Vi 
terms are convective non-linear contributions, and the 7, 
terms account for the density dependence of transport 
coefficients. Finally, C is the Langevin noise source with 
a correlation function 

(Ur,t)^(r',t')) =2g6 a p6(r-r')6(t-t'). (141) 



Equation (I140[) is supplemented by the continuity equa- 
tion (|127l) and t he relation (11291) between the potential 
and the density. lAlicea et al. Tiqoi) postulate Eq. (fT30)) 
on symmetry grounds: it includes leading terms of the 
expansion in gradients, time derivative, and amplitudes 
of the current and th e potential. It shoul d be pointed out 
that the assumption ([Alicea et all 120051) of local relation 
between j and <j) appears to be an oversimplification, since 
there are contributions to dtj(r, t) that depend in a local 
way on the density n e (r) and since Eq. (II 29[) linking the 
density to the potential is non-local. It remains to be 
seen how sensitive is the result to this assumpion. 

Equation (I140[) describes a transition from a conven- 
tional resistive state at r > to the ZRS at r < 0. Within 
the transition picture of Sec. IIV.B1 the transition is con- 
tinuous (second order), and the spontaneous current in 
the symmetry-broken phase is j = W\r\/u, i.e. j ~ \r\^ 
with the critical index f3 = 1/2. Clearly, this is just the 
Landau mean-field description of the transition, and one 
needs to fin d out how fluctuati ons affect these results. 
This is what lAlicea et Ml (l2005h do; their main findings 
are: 

(i) The symmetry-broken state is stable with respect 
to current-density fluctuations (there is no infrared diver- 
gence contrary to usual thermodynamic transitions where 
the Mermin- Wagner theorem implies a logarithmic di- 
vergence in 2D leading to destruction of the long-range 
order) . 

(ii) For the model with short-range interaction, U(r — 
r') = C~ 1 5(r — r') (that would correspond to a system 
with a screening gate) and with a symmetry with respect 
to <\> —> 0+const the gaussian fixed point of the mean-field 
theory is stable, with non-linear terms being marginally 
irrelevant. As a result, the mean-field exponents hold, 
implying the scaling j ~ \r\" with j3 = 1/2 on the ZRS 
side of the transition, as well as j ~ \E\ S with S = 1/3 
and the overdamped "diffusion mode" ~ l/(—iu) + Dq 4 ) 



at criticality, with logarithmic corrections to scaling. The 
invariance with respect to — > cp + const forbids in 
Eq. (|140j) those terms that explicitly depend on the value 
of 0(r) (rather than on its gradients). In the physical sys- 
tem such invariance does approximately hold; however, it 
gets violated when variations of density (linked to those 
of 4> via Eq. (|129|) become comparable to the total density 
of the electron gas. 

(iii) The long-range Coulomb interactio n drives the 
system away from the gaussian fixed point. lAlicea et al\ 
(|2005l ) analyze the RG equations in the model with l/fc c 
interaction (in momentum space) that has the upper crit- 
ical dimension d = 2 + e and find no stable fixed points in 
2D. This leads them to the conclusion that the transition 
becomes first order in the presence of long-range interac- 
tion. Presumably, this should also hold in the physical 
case of 1/k interaction, i.e. e = 1. 

(iv) The terms not included in the — > + const 
model drive the transition first order as well (whether 
with short-range or with long-range interaction). These 
results are obtained by the analysis of the RG flow d\jc~e 
dimensions, where the upper critical dimension due is 
found to be due = 4 and due = 7 for short- and long- 
range interactions, respectively. In both cases no stable 
fixed points is found, implying that the transition gets 
first order. This conclusion is then extrapolated to the 
physical dimension d = 2. 

(v) lAlicea et all ( 20051) also point out a similarity be- 
tween the ZRS problem and a phase transition in the 
2D bird-flocking model. The latter is essentially a dy- 
namical generalization of the 2D XY-model: birds move 
with a constant velocity and at each time step every 
of them picks up a new direction governed by average 
yelocity of surrounding birds, with some l evel o f noise 
dGregoire arid Chatel 120041: iToner and TuL Il995l Il998t 
IVicsek et all 119951 ) . For a sufficiently weak noise the ro- 
tational symmetry gets spontaneously broken: the bird 
flock acquire s a m acroscopic collective velocity. While 
IVicsek et al. (1995) foun d a second-order phase tra nsi- 
tion, a more recent work ( Gregoire and Chate . 20041) ob- 



tains a first-order transition. (Alicea et al 



2005) find 



that the ZRS and the bird-flocking model are almost 
equivalent, up to the Lorentz-force term (present in the 
former but not in the latter) , so that the first-order tran- 
sitions in both of them would be mutually consistent. 



V. HALL FIELD- AND PHONON-INDUCED 
RESISTANCE OSCILLATIONS (HIRO AND PIRO) 

A. HIRO: Experimental discovery and basic properties 



A decade ago lYang et al ] j2002a) discovered promi- 
nent l/_B-oscillations of the differential resistivity r = 
p + I(dp/dl) in a 2DEG with moderate mobility [i ~ 
10 6 cm 2 /Vs subject to a constant direct current /. These 
oscillations, termed the Hall field-induced resistance os- 
cillations (HIRO), can alternatively be observed at fixed 
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B and varying 7 . The latter approach allows one to ob- 
serve more oscillations since the amplitude of the DOS 
modulation (fixed by the value oj c r q ) does not change 
during the measurement. Examples of both realiza- 
tions of HIRO in a high mo bility (p ~ 10 7 cm 2 /Vs) 
2DEG are shown in Fig.HS lYang et all (|2002al ) pro- 




0.5 1.0 

B (kG) 



FIG. 29 (a) Differential magnetoresistivity r showing four 
HIRO maxima marked by integers measured by sweeping B at 
fixed 7 = 80 pA. (b) Correction to the differential resistivity 
5r versus ed c obtained by sweeping I at fixed B = 557 G. The 
data were obtained at T ~ 1.5 K in a 100 )im wide Hall bar 
sample with n e ~ 3.7 x 10 11 cm" 2 a nd p ~ 1.2 X 10 7 cm 2 /Vs. 
Adapted from lZhang et all (|2007bh . 

posed that HIRO stem from the commensurability be- 
tween the cyclotron diameter 2R C and the spatial sep- 
aration Ay = u c /eE between the Hall held-tilted LLs. 
Here, E = pnj is the Hall electric held, j = I/w is 
the current density, and w is the Hall bar width. This 
yields oscillations with edc = 27? C /Ay = 2eER c /uj c with 



-1/2 



on the elec- 



a characteristic dependen ce edc o 
tron density, confirmed bv lYang et all In addi- 



tion to single-layer 2D ele ctron systems, HIRO were also 
observed in double QW s (iBvkovl l2008d: iMamani et all , 
2009j; IWiedmann et al. |, l2011ah and in 2D ho le systems 
based on carbon-doped GaAs/AlGaAs QWs ( Dai et all , 
2009). HIRO were reported to retain their dc charac- 
ter under ac excitat ion with frequencies up to 100 kHz 
(|Bvkov et aZ.I . l2005bh . 



At 27re ( j c > 1, the oscillatory part Ar of HIRO 



Fig. is well described by 



Ar/ p ~ aS 2 cos(27redc 



(142) 



In agr e ement with other experiments (jHatke et all , l2010l 
l2009bLl2011el) . the HIRO maxima and minima in Fig. 1291 
occur at integer and half-integer values of ed Cl respec- 
tively. The amplitude a, obtained by IVavilov et all 
(2007) for edc 3> 1 in the form 



16t/tTT„ 



(143) 



is proportional to the backscattering rate 1/tv, see 
Eq. (|150[) below. HIRO can thus be understood in terms 
of resonant inter-LL transitions resulting from backscat- 
tering off impurities, i.e., HIRO require a sufficient 
amount of short-range scatterers for their observation. 
Dingle-plot analysis using Eqs. (|142[) . (|143|) allows one to 



extract both the quantum lifetime r q and the backscat- 
tering time tv ~ T sh . For a 2DEG with p ~ 10 7 cm 2 /Vs, 
such analysis at low T produced r q ~ 20 ps and r n ~ 5r 



(lHatke et al. 



lHatke et al 



2009b). 



( 2009bl ) demonstrated that, similar to 
MIRO, the decay of HIRO at high T is primarily due 
to the LL broadening induced by electron-electron in- 
teractions, which gives S 2 = e ~ 27r/ " cT q in Eq. (I142p with 
V T q( r ) - l/r q oc T 2 , see Secs. HIEAimi.B.l.fl At low T, 
the experimentally observed T dependence of l/r*(T) 
deviates from the T 2 law, which can be related to the 
heating of 2DEG by the direct current. The heating may 
also explain the decay of HIRO with edc at constant 73, 
see Fig. [29b. According to Eq. (jl42[) . the decay likely re- 
s ults from a suppressi on of 5 at 7 increases. 

lHatke et all (|2011eh found that the HIR O amplitude 
is strongly suppr essed [similar to MIRO ( Manil . 120051 : 
lYang ef aUl2006h (see Sec. IIII.C3p ] by the in-plane mag- 
netic field 7?|| ~ 1 T, while the HIRO period depends 
only on the perpendicular component. The experiment 
demonstrates that the suppression should be attributed 
to a Bu -induced enhancement of the quantum scattering 



rate 1/r* - 1/r* 



Bn=0 



oc B 2 . 



In the regime of strongly 



developed SdH oscillations, Bu starts to manifests itself 
also via the Zeeman splitting, see Sec. IVII.A.5I 



B. Theory of nonlinear dc transport 



Similar to the mechanisms of MIRO (Sec. IIII.B[> . two 
most important contributions to the nonlinear electric 
current in the absence of microwaves come from the di- 
rect effect of the dc field on impurity scattering (displace- 
ment mechanism) and from nonequilibrium changes of 
the isotropic part of the distribution function (inelastic 
mechanism). For the case of smoot h disorder, these ef- 
fects were systematically studied by IVavilov and Aleirier 



(120041) (displacement mechanism) and by Dmitricv et al. 
(|2005l ) (inelastic). Comprehensive analysis for a generic 
model of disorder [Eq. (f7T|) ] fo r high T ^> u) r and over- 



lapping LLs was performed in (|Vavilov et all . [20071 ). 

Below we overview the main results of these works us- 
ing the formulation of Sec. IIII.BI (which, in turn, fol- 
lows from the general quantum kinetic approach outlined 
in Sec. IIII.C[> . In the absence of the microwave field, 
Eqs. ([H5 ]) -([rJ5 1) reduce to 



j d = 2u e 



(144) 



where e' = e + eEAX vv i and the angle brackets denote 
averaging over tp and tp'. The isotropic part / of the 
distribution function obeys 



(fe ~ /J) /71a = (T-^(fe ~ fe>))^ 



(145) 



We restrict further analysis to 2ir 2 T/ui c 3> 1 (ther- 
mally suppressed SdH oscillations) and w c r q <C 1 (over- 
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lapping LLs). In t his limit, Eqs . (11441) and (|145[) yield 
the main result of ( Vavilov et all |2007|) , 



3d 



1 - 2£ 2 r 7 "(C) + 2S 2 r^T(r ia , C) , (146) 



<JdE ' ^' ( 

where 7(C), for t~ , parametrized by Eq. (|71|) . reads 

7(C)=E ^ lj '(0, C = ^dc = ^^, (147) 

— <n uj c 

with J n (C) being the Bessel functions. The second (dis- 
placement) term in Eq. (|146p comes from Eq. (|144p when 
substituting the equilibrium distribution for f e . The 
last (inelastic) term is proportional to the amplitude 



-2C 7 , (C) 

of the oscillations of f e in Eq. (|145[) , 



U =fi +<5sin 



. (2-ke 
in 



-F(r in ,C)ga £ /J 



(148) 



(149) 



In the large-C limit, Eq. (|146p is significantly simplified: 



3d 



a D E + 2S 2 a D E 



At 



sin 2C, 



(150) 



where the backscattering rate r^ 1 = t~ 1 e~xp(inir) 
(the actual condition on £ depends on the type of disor- 
der, see below). The nonlinear current is then dominated 
by the displacement term and shows oscillations propor- 
tional to sin 2^, which yields the cos2£ dependence of 
the differential resistivity observed in the experiment, cf. 

Eq. 032]). 

To explain the HIRO phase, let us look at Eq. (|144j) 
for T 3> ll>cC ~ eER c and neglect the oscillations of f £ 
by assuming f e = fj. Then, Eq. (|144[) is expressible as 
an energy averaging over the period lj c : 



u / {Ax vtp ,y I 



2v e 2 E 



£ / <p<p' 



(151) 



In overlapping LLs, (v £ v e +u>)e = 1 + 2S 2 cos(2ttll>/lu c ). 
The angular integrations in Eq. (|151l) within the 
stationary-phase approximation in the limit E — > 00 
yield Eq. f|150[) . The main contribution to HIRO comes 
from the backscattering processes with ip ~ ±7r/2 and 
\ip — tp'\ — 7r. These correspond to the maximum possible 
shift AX VV > ~ 2R C along the electric field. The average 
(v e v e +2eER c }e is maximized at integer e^c = 2eER c /uj c , 
which results in HIRO maxima in the differential resis- 
tivity at integer £dc, cf. Eq. (|142[) . 

The amplitude of the resistivity oscillations in 
Eq. (|150p does not depend on E and T, while experi- 
ments show considerable suppression for strong current 
and/or high T. As discussed in Sec. IV.Al this suppres- 
sion can be explained as a result of additional broadening 
of LLs by electron-electron interactions, combined with 



the effect of heating of the 2DEG by the dc field which 
is controlled by electron-phonon in teractions. 

To analyze the obtained results, IVavilov et ali (|2007l ) 
introduced the mixed disorder model (|93p (used in 
Sec. IIII.Bl to describe MIRO). The smooth (character- 
ized by 7~q~g m ) and sharp (t^ 1 ) components of disorder 
yield two separate contributions to 



7(0 = 



Tsh 



(152) 



The contribution of the sharp component changes on 
a scale of £ ~ 1 and shows strong oscillations for 
£ > 1 . The se cond term describes the n onlinear ef- 
fects, studied by IVavilov and Aleinerl ( 20041 ). that result 
from small-angle scattering by the long-range component 
(x 1 ^ 2 "C 1 gives a typical value of the scattering angle, 
see Sec. MLB]) . These become relevant in a much stronger 
dc field, £ ~ X -1 / 2 ^> 1, and do not contribute to HIRO. 
At C < 1, Eqs. (fT46]) and (fl"48]) reduce to 



3d 



a u E 



1 + 25 2 



3r 



4r* 
2r in ( 2 /r 
l + r in C 2 /2r ' 



S 2 C 2 -2^(r in ,C) ,(153) 



(154) 



where 1/r* ~ 3/r sh + 12x 2 /"r qiSm and 1/r ~ l/r sh + 
x/ r q,sm- Note that, in this regime, the nonlinear cor- 
rection to the Drude formula coming from HIRO can be 
equivalently obtained by taking the limit of uj — > in 
the MIRO terms in Eqs. ([ZZl), ([82]). and ([84]). At or- 
der E 3 in jd ! the ratio of the inelastic and displacement 
contributions in the case of smooth disorder is of order 
T- m T*/r 2 ~ T in /r q > 1 [similar to MIRO, Eq. ([TOO])]. 
That i s, the inelastic term, obtained for smooth disor- 
der bv lDmitriev et all (|2005l ). dominates and Eqs. (|153p 
and (|154p are valid for all £ <C x The suppression 
of jd/ E in a relatively weak dc field at £ ~ (t/"^) 1 / 2 is 
the strongest effect in the nonlinear dc response for the 
case of smooth disorder. 

It is important to mention the dual role the electric 
field plays in Eqs. ([148]) and ([154]) . On the one hand, 
it creates the oscillations (|149p in the energy distribu- 
tion of electrons; on the other, it also opens an addi- 
tional channel of inelastic relaxation, thus controlling the 
magnitude of the oscillations. Indeed, in the presence of 
the dc field, an electron can change its kinetic energy 
via elastic collisions with impurities. At £ » (T/Tin) 1 ^ 2 , 
the resulting "spectral diffusion" becomes more efficient 
than the inelastic relaxation due to electron-electron col- 
lisions. This leads to the saturation of the oscillations in 
Eqs. (|154p and (|149p and makes the nonlinear response 
independent of T for £ ^> (t/t^) 1 / 2 . Due to the spec- 
tral diffusion, J-"(Ti n ,£) remains smaller than or of the 
order of unity for arbitrary £. Note that the oscillatory 
term in Eq. (|149p contains additionally two small factors 
d and uj c d E fJ ~ —lu c /T and, therefore, remains small at 
any £, which makes the expansion (I149P legitimate. Un- 
like Eq. I|153p . which contains the "pure" displacement 
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and inelastic contributions, the last term in Eq. (jl46[) is, 
strictly speaking, a result of the interplay of the inelastic 
and displacement mechanisms in the strongly nonlinear 
response. 

Let us now consider the case x T sh "C Tq )S m *C T s h and 
£ <C X _1 ^ 2 j when the short-range component of disorder 
determines the transport scattering rate, i.e., T s h ~ t, 
while the nonlinear corrections generated by the smooth 
component are negligible. Equations (|146[) and (|148|l 
then reduce to 



3d 



1-25 2 [J 2 (C)]"- 



166 2 J 2 (()J 2 (Q 
r/r in + 1 - J 2 (C) 



(155) 



In this limit, r qjSm ~ r q enters jd only via the Dingle 
factor S — exp(— 7r/w c r q ). The above expression captures 
both the high-field limit (|150|) for HIRO and its essential 
modifications at intermediate (C ~ 1) and small (C ~ 
\J t /t\xi) electric fields. Unless T m <C r (which may occur 
at elevated T), the inelastic and displacement terms in 
Eq. (|155[) are equally important at £ ~ 1. Provided r s h ~ 
r, the smooth component only slightly modifies Eq. (|155[) 
in the relev ant range of C -C Y " 1 ^ 2 , see Eq. (3.8) and 
Figs. 2,3 in (jVavilov et a ll . 120071 ) . 



C. Nonlinear resistivity: Inelastic effects 

In accordance with the theoretical predictions, the 
most pronounced inelastic effects in the nonlinear re- 
sistivity were observed in 2DEGs with high density 
10 12 cm" 2 and modera t e mobility ji ~ 10 6 cm 2 /Vs 



(IBvkov et all l2005bl l2007t I Kalmanovitz et 



Mamani et aR 120091 : IVitkalovL l2009t IZhang et all . 



MMflp^i i r I ■ ■ i; i tj I Q 17 i 17 

2007al ). Indeed, according to Eqs. (fl~53"j) and ([TBI)) , the 



2008; 



2009 



differential resistivity at order £ 2 re 



r = p B (l + 2S 2 )-12S 2 p B 



3t 
16t* 



(156) 



The maximum value of the factor 3r/16r^ in the in 
the displacement term is 9/16 (which corresponds to 
t = T s h, see Fig.[T3"a). while the factor Ti n Jr in the 
inelastic contribution can be estimated as Ef/T 2 t ~ 
130(1K/T) 2 (for fi = 10 6 cm 2 /Vs and n e = 10 1 2 cm 



For r sh ~ 5t found from HIRO by lHatke et all (|2009bf ) 
(see Sec. lV.A[) . 3t/T6t* ~ 0.15 and the ratio of the 
two factors is 16T*r in /3T 2 ~ 9 x 10 2 (1K/T) 2 . This 
justifies the anal ysis of the low - field n onlinear resistiv - 
ity performed bvlMamani et all (|2009t) : IVitkalovl (|2009l) : 
IZhang et all (|2009l . l2007afT solely in terms of the inelas- 
tic mechanism [and also justifies the retaining of the 
high-order terms in Eq. (|154| )]. All these experiments 
reproduced the theoretical predictions in the range of 
applicability of the theory 17 and reported the values of 



17 The analytical results by l|Dmitriev et al 1. 120051: IM amani et ali . 
2009) are applicable for 2tt 2 T/lj c 3> 1 and overlapping 



ri n oc T 2 close to th ose calculated by (IDmitriev et all . 
2005}). In p articular. iMamani et all i|2009) generalized 



the theory ( Dmitriev et al . 2005( ) to the two-subband 
case and found that the MISO peaks in the nonlinear re- 
sponse arc inverted at ( = (2T/3r in ) 1 / 2 , which was used 
to accurately determine r; n . The effect is similar to the 
interplay of MISO and MIRO discussed in Sec. IIII.C4l 
In the nonlinear dc resp onse, interaction of MISO and 



MIRO was observed by IWiedmann et all ( 2011al ) and 
explained therein in terms of the inelastic mechanism. 
The data (Wiedm ann et all . l2011al ) indicate also a two- 
subband counterpart of the nonlinear mixing of HIRO 
and MIRO described in Sec. lVLBl 

In the ultra- high mobility 2DEG with n e . — 3.7 x 



by 



lO^cm-^ a n d n ~ 1.2 x 10 7 cm 2 /Vs used 
(|Zhang et all . l2007bf) . the estimated ratio T- m jr ~ 
Sf/T 2 t ~ 4(1K/T) 2 is much smaller than in the sam- 
ples discussed above and, as a result, both contributions 
in Eq. Q156p are relevant. In this case, the Dingle analysis 
of the low-fiel d differential resistivity at fixed T — 1.5 K 
(|ZudovL 1201 ll) yielded a value of 2.15 for the expression 
in the brackets in the nonlinear ter m of E q. (1156|) . Using 
Eq. (SI]), for r = 20r q and r sh ~ 5r (|Hatke et all . l2009bl ) 



one obtains 3r/16r* ~ 0.15 which divides almost equally 
between the sharp and smooth components of disorder, 
see Fig.[T3k. It follows that the inelastic term Ti n /r ~ 2 
at T — 1.5 K, which agrees with the estimate above. This 
example shows that measurements of the nonlinear mag- 
netoresistivity in high LLs in both limits ( > 1 (HIRO) 
and C <C 1 [Eq. (|156p ] provide a method to determine 
various scattering rates — in particular, T-dependent r q 
and ri n , as well as and 7V — which bring valuable in- 
formation about disorder and interactions in 2D electron 
systems. 



D. Zero differential resistance states (ZdRS) 



Expe riments by IBvkov et al. | (|2007l) : IZhang et all 
(2008b) revealed that the differential resistance can drop 



all the way to zero (Fig. [30]) . leading to the forma- 
tion of the zero-d i fferential resistance states (ZdRS). 
Bykov et al. | (|2007l ): iKalmanovitl et al. (2008) (in sam- 



ples with n e ~ 8 x 10 11 cm~ 2 and fi ~ 8 x 10 5 cm 2 /Vs) 
observed ZdRS emerging from the maxima of SdH oscil- 
lations. The transition to ZdRS was accompanied by a 
reproducible negative spike (Fig. 131)10 and by temporal 

fluctuations at higher dc bia s. 

By contrast, I Zhang et all (|2008bl ) (in a sample with 
n P ~ 4.8 x 10 11 cm " 2 and \x ~ 4.4 x 10 6 cm 2 /Vs) and 
lHatke et all (|2010f) (in samples with n e ~ 3.8 x 10 11 cm" 2 



LLs. In the regime of SdH oscillatio ns (2tt 2 T/lj c < 1) , 
which was also studied e x perimentally by IBvkov et ali (20071 ): 
IKalmanovitz et al. (2008'): IZhang et ali (2007a), such a descrip- 
tion is not parametrically justified and additional theoretical 
analysis is required. 
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FIG. 30 (a) Differential resistivity r versus current 7 for 
several values of T at fixed B = 7.64 kG, measured in a 5 (jm 
wide Hall bar sample. Adapted from lBvkov et all (|2007f ). (b) 
Similar dependence at B = 1.7 kG and B = 2.5 kG, measured 
at T ~ 1.5 K in a 100 /jm wide Hall bar sample. The first 
HIRO maximum is marked by f . The c urve for B = 1.7 kG 
is shifted by 2fi upwards. Adapted from lHatke et al\ (|2010h . 



and /j, ~ 1.0 x 10 7 cm 2 /Vs) demonstrated the possibility 
of ZdRS evolving from the principal minimum of HIRO 
in the regime of suppressed SdH oscillations. Neither 
overshoot to negative valu es nor temp o ral flu ctuations 
were detected (Fig.l3"0"b). lHatke et al\ (|2010| ) reported 
ZdRS over a continuous range of currents and in mag- 
netic fields extending well below the onset of SdH os- 
cillations. Similar to ZRS, ZdRS were found to disap- 
pear with increasing temperature (T > 2 K) and with 
increasing overlap between LLs (B < 1 kG) . The mini- 
mum current I\ (Fig. 130b) required to support ZdRS was 
found to be roughly B-indcpcndcnt. The maximum cur- 
rent li was found to increase roughly linearly with B, 
tracing the fundamental HIRO peak (cf. f in Fig. [30b). 
ZdRS wer e also observed to develop from the maxima 
of MI RO (IZhang et all l2007d) and MISO ([Bvkov et all 
l2010dl: iGusev et aUl201ltlWiedmann et a/ll2011al) . 

On the theoretical side, ZdRS originate from the parts 
of (local) I-V characteristics where the differential resis- 
tivity is negative, which violates the stability condition 
([j"3l~]l . Similar to the case of ZRS (Sec.H7|), this results 
in an instability leading to the formation of current and 
field domains. Unlike ZRS, however, the absolute resis- 
tivity remains positive, which leads to an essential dif- 
ference between the Corbino (N-shaped I-V characteris- 



tics) and Hall bar (S-shaped I-V characteristics) geome- 
tries. While in the Hall bar geometry a stationary domain 
configuration wit h the current flowin g along the domain 
walls is expected (jBykov et all 120071) , in the Corbino ge- 
ometry nonstationar y (mov ing) domains are predicted 
(|Vavilov and Aleinerl . \2004) similar to the Gunn effect 
Gunnl . Il963h~ 



E. PIRO: Experimental discovery and basic properties 

Resonant interaction of 2D electrons with lon- 
gitudinal optical ( LO) phonons was predicted 
( Gurevich and Firsovl . Il96ll ) and confirmed in mag- 
netotransport measurem e nts i n GaAs/AlGaAs het- 
erostructures (|Tsui et all [l980) long time ago. This 
interaction manifests itself as an enhancement of the lon- 
gitudinal resistivity whenever the LO-phonon frequency 

wlo — nuj c , n = 1,2,3, In GaAs, wlo ~ 10 13 s~ 1 

and observation of the LO phonon-induced oscillations 
requires hig h T > 10 2 K and strong B > 10 2 kG 
(|Tsui et alllimh . 
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FIG. 31 (a) Magnetoresistivity p as a function of B, show- 
ing eight PIRO maxima, four of which are marked by inte- 
gers, (b) PIRO order n versus The linear fit yields 
s ~ 3.4km/s. The data were obtained at T = 3.5 K in a 
100 fim wide Hall bar sample with n e ~ 3.8 x 10 11 cm" 2 a nd 
fj, ~ 1.2 X 10 7 cm 2 /Vs. Adapted from lHatke et al\ (|2009al ). 

A decade ago, another class of phonon-induced oscilla- 
tions, termed the phonon-induced resistance oscillations 
(PIRO), was discovered in the linear-respo nse resistivity 
of a 2D EG with mobility // - 10 6 cm 2 /Vs (|Zudov et all 
l2001b[ ). PIRO emerge at much lower T ~ 1 - 10 K and 
much lower B ~ 1 — 10 kG, and are understood in terms 
of resonant interaction of 2D electrons in high LLs with 
acoustic phonons which carry momentum 2kp and have a 
characteristic frequency = 2ItFS. Here s is the speed 
of sound and the out-of-plane component of the phonon 
momentum is neglected. Such interaction causes a cor- 
rection to the resistivity Ap ph which oscillates with the 
ratio 



e p h = ^tt/^c oc nl /2 . 



(157) 



The nV" 1 dependen ce was confirmed experimentally 
(|Zudov et a/lboOlbJl . 
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Electron backscattering by an acoustic phonon is most 
effective when w T = nu c , n = 1,2..., which maximizes 
the (thermally averaged) product of the initial and final 
densities of states. In the simplest model of 2D isotropic 
phonons, the oscillator y part of the res i stivity for e p h > 1 
and r>u c , cjjr reads ( Dmitriev et 



F. PIRO in a strong Hall field 



ApP h 2 9 2 Tt x2 



PD 



S cos(27re p h — 7r/4) 



(158) 



[see also Eg. (|161[) ] . The phase — 7r/4 was con- 
firmed experimentally by lHatke et al I (1201 ldT) m a va- 
riety of high mobility samples. Many other exper 



lments 



Hatke et al. 


, 2009a 


2008ri Zudov et al. 



v and Goranl. 
Yanj 



20091: iBvkov et all 12005. 



Yang 
120091, L 



et al 



2002b; iZhang et all 



2001bl ) reported PIRO max- 



ima at integer e p h- As discussed in Sec. IV.Gi the phase 
of PIRO (unlike that for MIRO and HIRO) is not ex- 
pected to be universal and is sensitive to the anisotropy 
of relevant phonon modes, crystallog raphic orienta t ion of 
the sample, wi dth of the QW, etc. ( Dmitriev et aUl20ld : 
lRaichevi . 120091) . 

In contrast to other \ow-B magnetoresistance oscilla- 
tions (SdH oscillations, MIRO, MISO, and HIRO), which 
are best observed at low T, PIRO are best resolved 
at T ~ ujir and get strongly suppressed at both low 

' 20091 lHatke et all l2009al: 
2001bh. At lowT, both 



and high T ( Bykov and Gora 
lYang et~al\. 



z 



2002b; Zudov et al 



emission and absorption of 2fc^-phonons is exponentially 
suppressed by the thermal factors for available electron 
and phonon states, see Sec. IV.Gl On the other hand — 
similar to MIRO, MISO, and HIRO— PIRO are insen- 
sitive to the temperature smearing of the Fermi surface 
which suppresses SdH oscillations at 2it 2 T/uj c 3> 1. A de- 
tailed study of the T dependence of PIRO revealed that 
higher-i? (lower-e p h) oscillations are best developed at 
progressively higher T: the temperature at which a given 
oscillation reaches its maximum amplitude was found 
to sca le with \/"B( Bvkov and Goranl l2009t lHatke et all 
l2009af ). This suggests that, at high T, the suppression of 
PIR O (lik e MIRO and HIRO, see Sees. iHTAllIII.B.Lfl 
and IV.A1 but unlike the SdH oscillations) is due to 
interaction-induced broadening of LLs. 

Similar to HIRO and MIRO, PIRO are most pro- 
nounced in high-quality 2DEGs. Early experiments in 
moderate mobility samples (/i ~ 10 6 cm 2 /Vs) required 
temperatures in the r ange about 5 — 20 K and revealed 
only a few oscillations ( Bykov et all [2005a; Z udov et all 
l2001bl ). More recent experiments in high mobility sam- 
ples (fj, ~ 10 7 cm 2 /Vs) showed up to eight oscilla- 
tions (Fig.|5T1) which remained visible down to T ~ 2 K 
(jHatke et all l2009aT) . 

Related phonon-induced oscilla tions in the magne - 
tothcrmopower were observed by IZhang et al] (|2004l) . 
Interplay of PIRO and MISO was investigated by 
IBvkov et ali ( 2010al ) in a single GaAs Q W with two pop- 
ulated subbands. 
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FIG. 32 Differential magnetoresistivity r as a function of B, 
normalized to the resistivity at B = 0, for different I (a) from 
2 to 50 /iA in steps of 4 fiA at T = 10 K and (b) from 30 fiA 
to 150 i-iA in steps of 10 fiA at T — 6K. The traces are offset 
for clarity. Integers I, II, and III at the top traces mark the 
order of the HIRO peaks. The vertical lines correspond to 
integer e p h = 1,2,3. The curves labeled by letters A, B, and 
C in part (a) are commented in the caption to Fig. 1331 The 
data were obtained in a 50 fim wide Hall bar sample with 
n e ~ 4.8 x 10 11 cm" 2 an d (i~4.4x 10 6 cm 2 /Vs. Adapted 
from I Zhang et all (|2008bTl . 



Recent studies (| Zhang et all l2008bl ) of PIRO in a 
strong Hall field established that the dc field leads to the 
evolution of the PIRO maxima into minima and back, 
see Fig.[32j It was also found that the strong dc field 
enables detection of PIRO at low T where the linear- 
response PIRO are exponentially suppressed. One more 
experimental finding was the observation of a pronounced 
resistance maximum at edc = £ph, where the Hall veloc- 
ity equals the speed of sound. All these results are re- 
produ ced within the theoretical model of lDmitriev et al\ 
(|2010ri . see Sec. lVXfl 



G. Microscopic theory of PIRO 



In the case of MIRO (Sec. IIII.Bj) or nonlinear dc trans- 
port (Sec. IV. B"|) , the inelastic effects that lead to oscilla- 
tions in the energy distribution of electrons play an es- 
sential and often dominant role. By contrast, PIRO are 
observed at elevated T, where a fast inelastic relaxation 
makes effects of this type much less important. There- 
fore, one can consider only effects that are similar to the 
displacement contribution to HIRO and MIRO. In terms 
of migration of the guiding centers of cyclotron orbits 
[Eq. (|65p - (|rJT)) ]. the phonon-assisted dissipative current 
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has the form of Eq. (j65j) with 



W X1 - 



dei / de 2 M eie2 5 (x 1 - x 2 + AX V>1V>2 ) 



where r^ s f^ 2 denotes the probability of spontaneous emis- 
sion of a phonon with frequency W12 and M u is the Planck 
distribution function. In the factor M. ee i [Eq. (|67p], one 
can substitute, in accordance with the above, fj for f e . 
Similar to Eq. (|66|) . the delta function containing Ae^' 
ensures energy conservation. 

We start with a simplified single-mode model which as- 
sumes interaction with 2D isotropic acoustic phonons via 
a deformation potential. We also assume that the out-of- 
plane component of the phonon momentum is negligible, 
i.e., the width b of the QW to which the 2DEG is confined 
is large, b^> kp. In this case, 



tot , N 

12 TW12 



(159) 



r 



(sp) _ „2 



g W12/2, g = mP /pfes" 



(160) 



where phonon frequency for quasielastic scattering is 
given by u)\2 = 2kpssm\(ipi — ip2)/2\, p is the mass 
density, and T> is the deformation-potential constant. In 
overlapping LLs and at high T > ■ u} r ,uj v , Eqs. (j6"5"j) . (|159l) . 
and (TT6Tjj) yield (jDmitriev et all . l201dh 



=^D^r« ? 2 r(l + 2<5 2 [J (27re ph )-J 2 (2^ ph )]) . (161) 

The linear-in-T dependence of is due to the fact that 
at T > w T the occupation number for relevant phonon 

modes is large, so that T^fl 2 Afu 12 ~ g 2 T/2, while the 
contribution of spontaneous emission is negligible. The 
current shows oscillations with e p h = w»/w c , controlled 
by commensurability between the cyclotron energy and 
the maximum possible — in the process of scattering — 
phonon ener gy. That is, in accor d with the original in- 
terpretation ( Zudov et all . l2001bh . PIRO originate from 
resonant inter-LL transitions caused by the backscatter- 
ing of electrons by acoustic phonons. The position of 
higher-order maxima of PIRO in Eq. (|161|) at e p h S> 1 is 
given by Eq. (|158jl . In narrow QWs (kpb ~ 1) or in wide 
QWs for e p h 3> (fci?6) 2 , one should take into account the 
3D character of acoustic phonons. Still in the isotropic 
approximation, the cos(27re p h — tt/4) behavior of PIRO 
[Eq. (|158[) ] changes then to cos(27re p h), while the ampli- 
tude of PIRO reduces by a factor of /2bkp (jRaichevI . 

While the isotropic single-branch model captures the 
essential physics of PIRO, in real structures the electron- 
phonon interaction is more complicated. In bulk GaAs, 
there are three anisotropic phonon branches and two 
mechanisms (via the deformation and piezoelectric po- 
tentials) of electron-phonon interaction. A comprehen- 
sive study of PIRO, using the general form of inter- 
action with bulk acoustic phonons, was performed by 



iRaichevi (2009). For higher harmonics of PIRO, the ef- 
fects of anisotropy were treated there analytically, pro- 
ducing three distinct contrib utions with different phases 
and periods (jRaichevI I2009T ) . Additional contributions 
can arise from interaction with interface phonon modes 
(jZudov et al!l2001bh. 

The analysis bv lRaichevI ( 20091 ) shows that PIRO in the 
ultra-hi gh mobility a n d mod erate electron density sample 
used bv lHatke et all ( 2009at ) are dominated by the bulk 
transverse acoustic mode propagating along the 2D plane 
and polarized perpendicular to the plane. The sound ve- 
locity of this mode sto — 3.40 km/s is indeed very close 
to s = 3.44 km/s ext racted from the period of PIRO by 
(jHatke et all . l2009al ). A subleading contribution of the 
longitudinal mode with a noticea bly higher velocity p ro- 
duces an extra peak observed bv lHatke et all ( 2009al ) at 
sufficiently high B, and also explains the beating pattern 
at lower B, caused by the interfere nce of these tw o con- 
tributions to PIRO [see Fig. 1 of (jRaichevI l2009h ]. By 
contrast, in the moderat e mobility and high e lectron den- 
sity sample studied by (B vkov et all . l2005aj) , PIRO are 
dominated by the higher-e nergy longitudinal m ode. In 
another recent experiment ( Zhang et all . l2008rj) , which 
studied the effects of a strong dc field on PIRO (dis- 
cussed below), the longitudinal mode also prevail, which 
explains the good a greement between the experiment 
(Zha ng~e£ all . T2008bl ) and the theory ( Dmitriev et all . 
|2010D based on the single- mode approximation (|160[) . 

The Hall field tilts LLs and changes the commensu- 
rability condition for electron-phonon scattering. This 
results in oscillations with e± = £dc ± £ph, where £d c 
is the parameter that controls HIRO (Sec. IV.Bp . The 
change of the period of the oscillations by the dc field 
can be viewed as a Doppler shift of the phonon modes 
in the frame moving with the Hall velocity vh — cE/B 
across the field. In the moving frame, the electric field is 
absent while the phonon dispersion becomes anisotropic, 
s s — vh cos[(cpi + ip2)/2] (this formulation is particu- 
larly useful in the case of a complicated dispersion rela- 
tion of phonons). For 1 <C £dc < £ph and 1 <C e p h — £dc, 
Eqs. flBSJ), (HMJ), and ([TE0|> yield for the oscillating part 

if 



of ^ 



• (p,osc) 
J_d 



Ad 2 



g Tt 



I sin 27 



cos 27re_ 
\Adc|e-| 



(162) 



The evolution of the differential resistivity r = dE/djd 
with va rying e t \c in this regi me was studied experimen- 
tally bv lZhang eCaU (j2008bh . The results of calculation 
according to Eqs. (|159[) and (|160p . illustrated in Fig. 1331 
quantitatively reproduce the experimental data in Fig. 1321 
without fitting parameters. Note that the theoretical 
plot does not include HIRO which become relevant at 
higher dc bias, see Fig. [321 but are negligible in the re- 
gion of interest, B > 2 kG and I < 40 p,A. I n qualitative 
agreem ent with the experimental results of lZhang et all 
(2008b), similar oscillations with e ph — e dc were also ob- 
tained numerically, using the balance-equation approach, 
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FIG. 33 Calculated oscillatory part 
resistivity dE/djd versus B, in units of pog'^Tr, for several 
values of current / from pA to 50 pA in steps of 5 /J.A. The 
traces are offset for clarity. The vertical lines correspond to 
integer e p h. The traces for (A) 10 ^tA, (B) 15 pA, and (C) 
30 fiA correspond approximately to the values of / at which 
the third, second, and first PIRO peaks, respectively, disap- 
pear before evolving into minima as I increases. The curves 
are very similar to those measured in the experiment: cf. 
the curves in Fig.r3"2l labeled a s A (1 /xA), B (14/xA), and C 
(28 //A). From IDmitriev et all <|2010h . 



bvlLeil (j2f)08l ). 

In a stronger dc field, £dc > £ph, corresponding 
to a supersoni c Hall velocity vh > s, the theory of 
iDmitriev et all ( 20101 ) predicts that the phase of the os- 
cillations changes in the second term of Eq. (|162l) from 
cos27re^ to sin27re_ for e^ c > e p h. Note that at vh > s, 
the effective sound velocity s — vh cos[(<^i + (p2)/2] in the 
moving frame becomes negative for certain transitions. 
Emission of a phonon with energy to in such transitions 
increases the electron kinetic energy by A > at the ex- 
pense of the electrostatic energy which is decreased by 
A + uj. 

The physical meaning of the "sound barrier" is clear 
from the low-temperature behavior of the phonon- 
assisted transport. At T -C u K} stimulated emission 
and absorption of phonons are exponentially suppressed 
by the Planck factor A/L„ — exp(-w 7r /T) in Eq. (|159l) . 
On the other hand, spontaneous emission in the scatter- 
ing processes reducing the kinetic energy is also strongly 
suppressed by the factor M eie2 . However, in the sponta- 
neous emission processes that increase the kinetic energy, 
the constraint imposed by the factor M ElE2 is not effec- 
tive. As a result, at £dc > £ph the phonon-assisted trans- 
port due to spontaneous phonon emission survives even 
at T = 0. Specifically, for 1 <C £- <C £dc the phonon- 
induced oscillations at T < are described by 



. (p.osc) 

3d ' 

(JdEt 



2S 2 



1/2 



2 3/2 



-g w c sin 27re_ 



(163) 



The evolution of the oscillations with temperature is il- 
lustrated in Fig. |M] (where also the smooth part of the 
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FIG. 34 (a)-(c): Oscillatory part of the phonon-assisted cur- 
rent j^ ,OBcS> in units of 2S 2 g 2 u) c rar>E, calculated for the case 
of (a) T = 5lj c , (b) T = 0.7o; c , and (c) T = 0.25oj c . (d): 
Smooth part of the p honon-assisted current at T — 0. From 
IDmitriev etall (|2010l ). 



phonon-assisted current at T = is shown). 

Nonlinear mixing of PIRO and MIRO was investigated 
by iRaichevI (|2010bt ) . This study predicted oscillations 
with ujyr ± u and discussed the feasibility of the ir exper- 
i menta l observation. Further, the analysis by IRaichevI 
(|2010b[ ) showed that heating of the 2DEG leads to addi- 
tional a-v/wc-oscillations with a phase shift with respect 
to the equilibrium PIRO. Separately, nonlinear mixing 
of PIRO and M I SO in t wo-subband syst ems was shown 
by iBvkov et al\ (|2010al) ; IRaichevI (|2010al) to lead to os- 
cillations at frequencies determined by ui^ and the sub- 
band splitting energy. A further example of frequency 
mixing for the case of MIRO and HIRO is discussed in 
Sec. lVl.H.2l 



VI. TRANSPORT IN STRONG AC AND DC FIELDS 

In the preceding sections, we have discussed two types 
of nonequilibrium magnetooscillations, MIRO (Sec. lIII[) 
and HIRO (Sec.|V|. The nature of MIRO has so far 
been considered at the level of single-photon processes. 
We now turn to mode-mixing phenomena; specifically, to 
multiphoton processes in MIRO and to nonlinear mixing 
of MIRO and HIRO. 



A. Fractional MIRO 

1. Experiments at high microwave power levels 

In addition to MIRO, whose extrema in the vicinity 
of integer e ac values are described by Eq. (f63|) and the 
oscillatory structure is shown in Fig.[Bl similar oscillatory 
features have been observed near e ac given by certain 
rational fractions, with maxima (+) and minima (— ) of 



45 



the photoresistance at 



— T <£ac , 
TO 
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where n and to are integers (with nonintcgcr n/m) and 
<ys ac > iSj in typical experiments, considerably smaller 
The oscil l atory feature at e 



than 1/4 (Fig. _ 

1/2 wa s obser v ed bvlDorozhkin ( 20031); Dorozhkin et al. 
2007bl l2005h : IWiedmann et all (|2009a[) : IWillett e~ 



[in retrospective, the first signature of the half- 
integer MIRO can be recognized already in early ex- 
perime nts on moderate-mobility samples ( Zudov et all . 
l2001al) . see a weak oscillation at B ~ 2kG in the 
/ = 45 GHz trace there] . Similar fractional MIRO 
(fMIRO) for other fractions in the n = 1 series e ar = 
1/to were reported for to = 3, 4 ( Dorozhkin et all . 



(Wiedmann et al. 



l2007bt IWiedmann et all l2009af) and to = 5,6,7,8 



2009af ). Fractions that have been 
observed in the fMIRO series with n = 2 (apart 
from the even denominator fractions coinciding with 
the member s of the n = 1 series) are represented 
by to = 3 (IDorozhkin et al. . 2007bt IWiedmann et all . 



l2009at IZudovl. 12004 IZudov et all l2006al) and to = 5. 7 



(jWiedmann et all . l2009al ). Note that the fMIRO series 
with n = 1 , 2 belong to the high-B side of the CR and 
thus do not overlap with the integer MIRO which oc- 
cur at e ac > 1. A third series of fMIRO among those 
observed is, in contrast to the variable-denominator se- 
ries above, characterized by fixed to = 2 and variable n. 
The reported oscillatory features in the e ac = n/2 series 
that fall in the intervals between integer e ac are (apart 
fro m e ac = 1/2 mentioned ab o ve) those wit h n = 3 and 
5 (|ZudovL 12004 IZudov et all . I2003L l2006al) . The half- 
integer fMIRO with n > 3 lie on the low-£? side of the 
CR and therefore require for their observation (other ex- 
perimental conditions being equal) higher mobilities com- 
pared to the fMIRO series with n = 1. 

Similar to MIRO, the most prominent minima from 
all three fMIRO series, namely those in the vicinity of 
£ ac = 1/2, 2/3, 3/2, were shown to evolve into ZRS with 
incre asing microwave pow er in ultra-high mobility sam- 
ples (jZudov et aUl2006al) . As seen in Fig. 135a . the for- 
mation of these fractional ZRS was accompanied by the 
diminishing and narrowing of the neighboring photore- 
sistance peaks, as well as by the overall suppression of 
the photoresistance at e ac < 1/2. A similar strong sup- 
pression of the resistance by radiation on the high- B side 
of the CR was examined by IDorozhkin et al. I (|2005h . 



2. Multiphoton absorption via virtual states 

The theoretical framework developed to describe 
fMIRO relates them to multiphoton processes in which 
m photons are absorbed in transition s between elec- 
tron s t ates separate d in energy by nujr- dDmitriev et ~al 



200714 iLei and Liul. [2003) l2006bl iPechenezhskii et al. 
20071: 1 
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20071 ; iTorres and Kunoldl . 120051 ) , as illustrated in the in- 



FIG. 35 (a) Magnetoresistance under microwave illumination 
at / = 27 GHz for different radiation intensities. The traces, 
offset vertically for clarity, are labeled according to the atten- 
uation levels. The downward arrows mark the minima near 
e ac = 3/2, 2/3, and 1/2, which develop into fractional ZRS at 
higher radiation intensities. The inset illustrates multiphoton 
processes responsible for these features. The data were ob- 
tained at T ~ 1.0 K in a sample with n e ~ 3.6 x 10 11 cm -2 
and n ~ 2 x 10 7 cm 2 /Vs. From I Zudov et all (|2006al ). (b) 
Similar data for / = 17 GHz taken at T ~ 1.4 K in a sample 
with n e ~ 2.7 x 1 11 cm " 2 and fi ~ 1.7 X 10 7 cm 2 /Vs. From 
IDorozhkin et all (|2007bD . 



set to Fig. l3T>h . One can distinguish two (generally, inter- 
fering with each other) scattering channels for multipho- 
ton transitions: sequential absorption of single photons 
via real (resonant) intermediate electron states and mul- 
tiphoton absorption via virtual intermediate states. For 
overlapping LLs, both types of multiphoton transitions 
yield in the limit of a strong ac field oscillatory features 
in the vicinity of e ac = n/m — with an amplitude pro- 
portional at T > muj to 8 2ra — which, for given 6, are 
dampe d by disorder much more str ongly than the integer 
MIRO ([Pechenezhskii et al. I. I2007D . For separated LLs, 
however, single-photon real transitions are forbidden if uj 
is not sufficiently large to bridge the gap between the LLs 
(see Fig. [12]), so that microwave absorption occurs then 
by means of multiphoton transitions via virtual states. 
This is, in particular, the case for two-photon transitions 
between neighboring LLs at w = cj c /2 3> 1/t q . 



4G 



In the limit of well-separated LLs, the mechanism of 
multiphoton absorption via virtual intermediate states is 
much similar to that of mul tiphoton i oniza tion in atomic 
or semiconductor systems ( Keldvshl . Il965() . In particu- 
lar, a perturbative expansion in the radiation power is 
justified if u) is large compared to the transition rate be- 
tween LLs the slope of which oscillates in phase with the 
ac field. The parameter that governs the validity of per- 
turbation thj3OTy_is_(7} l A")P_^l, where V is given by 
Eq. (72 ) (|Dmitriev et all , l2007al) . Here and in the rest 
of Sec. IVI.Al we over view the approach develo p ed by 
iDmitriev et~cd\ (|2007bl) and IPechenezhskii et all (|2007r ) 
to describe fMIRO for the case of smooth disorder. Wc 
expect that the picture that results from this approach 
remains qualitatively correct for more general models of 
disorder. Similar to the inte ger MIRO, fMIRO ind uced 
by the inelastic mechanism ( Dmitriev et all . l2007bl ) are 
much stronge r than those induced by the displacement 
mechanism 18 (|Dmitriev et all, l2007bt iLei and Liul . 2003. 
l2006bl: iTorres and Kunoldl . 120051 ). provided r in /r q > 1. 
Under this condition, the half-integer fMIRO near e ac = 
n/2 with n = 1,3, 5,... are described at order V 2 for 



well-separated LLs ( Dmitriev et aZ.I . [2007bh by 



— = ^^(2^(1 + ^) 



Pph 

PD 
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where the function IZ2 (w), odd in the detuning 2uj — nco c , 
is defined in Eq. (85), and = 2\£ + £_\ 2 /T 2 depends on 
the polarization of microwaves; specifically, $ — [{ui 2 — 
uj 2 )/(uj 2 + u! 2 )) 2 /2 and d = for the cases of linear and 
circular polarization, respectively. For n = 1, Eq. I|165|) 
describes the two-photon CR in the photoresistivity. 

The two-photon contribution of the displacement 
mechanism is smaller (in the parameter r q /Tj n <C 1) but 
possesses, in the case of linear polarization, anisotropy 
with respect to the mutual orientation of the dc and ac 
fields. At order V 2 , the depth of the modulation, as the 
orientation is changed, between the principal values of 
the photoresistivity tensor induced by the displacement 
mechanism is given by the fact or 1 ± V2~&/(1 + ■&) (fo r 
the full analytical expression, see lDmitriev et all . 12007b). 
The orientation dependence thus provides a measure of 
the role of the displacement mechanism in fMIRO. 



18 The method of "fo r ce- and energy- balance equations" ijLeil . l2O10l : 
ILei and Li 3, 120031. I2006bl 1 tacitly assumes an instant thermal 
equilibration in the frame moving with the drift velocity, i.e., 
the inelastic mechanism of photoconductivity is lost within 
this formalism "by const ruction" . Similarly, the argument by 
ITorres and Kunoldl (|2005| ) that the "inelastic processes can be 
safely ignored" if they are slower than the elastic ones is in 
contradiction with the fact that without them the stationary 
regime of photoconductiv ity cannot be established: in effect, 
ITorres and Kunold (2005) also implicitly assume an infinitely 
fast equilibration of the same kind. 



3. Multiphoton absorption via sidebands 

The multiphoton absorption via virtual states 
[Eq. (I165[) ] is the main mechanism of fMIRO in the limit 
of well-separated LLs. A subleading (in the parame- 
ter l/uj c T q <C 1) contribution to fMIRO comes from 
single-photon transitions via radiation-induced "side- 
band s" formed in the gaps of the DOS between the 
LLs (|Dmitriev et all , l2007bl ). The sideband-mechanism 
can be viewed as resulting from disord er-induced transi- 
tions between the quasienergy levels (|Zeldovichl Il967l ) 
that characterize the stationary states of a homoge- 
neous electron system under the action of a periodic- 
in-time perturbation. A systematic approach to study- 
ing the spectral properties of a 2DEG driven by the 
ac field in the presence of disorder is formalized in 
terms of a time-dependent self-energy of the Wigner- 
transfo r med retarded Green's functi on ( Dmitriev et all , 
l2007allbl ; IVavilov and Aleinerl . l2004f ). The notion of 
a time-dependent DOS was also introduced to de- 
scribe sidebands in a closely rela ted prob l em of the 



1968) by 



dynamical Franz-Keldysh effect (lYacobv , . 
iJauho and Johnsenl ( 19961) ; iJohnsen and Jauhol (| 19981 ) 
In the fMIRO problem, the nonequilibrium DOS of each 
LL acquires, to first order in l/w c r q , two static satellite 
peaks 



(166) 



where v>^ N \e) is given by Eq. (27), centered for u = ui c /2 
in the middle between LLs, and also satellites centered 
at the sam e energies and oscillat ing in time with a pe- 
riod 2n/u) ([Dmitriev et al I l2007bl ). Note that the DOS 
peaks (|166[) can be probed by resonant spectroscopy in 
the same way as if they were present at equilibrium. For 
the half-integer fMIRO, single-photon transitions to and 
from the static sidebands yield a contribution which, at 
order V 2 , has the same shape as in Eq. (|165l) and is of or- 



der w c r q times smaller (jDmitriev et all , 2007b|. On the 



other hand, the time-dependent sidebands give a contri- 
bution which is smaller than that in Eq. (|165[) by only 
a factor of order (wcTq) 1 / 2 and has a distinctly differ- 
ent shape; in particular, it is even in th e detuning from 
the fM IRO resonance (for details, see IDmitriev et all , 
l2007bD . 



4. Multiphoton absorption via real states 

Multiphoton absorption via virtual states or sidebands 
provides the mechanisms of fMIRO in the case of well- 
separated LLs. If, however, gaps between LLs are opened 
(i.e., w c T q > 7r/2, Sec. II.C.l), but the width of the 
gaps A g < w c /2, then for any given u> in the interval 
A g < uj < uj c — A g real single-photon transitions are 
possible both between neighboring LLs and within the 
same LL. For w from the above interval, an electron can 
thus be transferred between neighboring LLs via two se- 
quential single-photon transitions. More generally, for 
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(n — 1)lu c + A g < uj < noj c — A g , single-photon real tran- 
sitions are possible with a change of the LL index by both 
7i—l and n, where n = 1,2,3,.... For given n, there- 
fore, there exists a finite range of u) and lj c within which 
m 7^ n photons can bridge the gap between n LLs via 
real intermediate states, which leads to the emergence 
of an additional mechanism of fM IRO ( Dorozhkin et all . 
l2007bt iPechenezhskii et all . |2007|) . fMIRO induced by 
this mechanism of sequential single-photon absorption 
are specific to the crossover between the regimes of over- 
lapping and well-separated LLs and disappear altogether 
for 2A g > oj c . At the crossover, however, they are the 
main fMIRO channel at order V 2 if Ti n /r q ^> 1, with 
an amplitude proportional to (r- m V /t) 2 , similar to the 
second-order expansion in V of the integer MIRO contri- 
bution in Eq. (97). 



B. Nonlinear mixing of MIRO and HIRO 

The physics of high LLs out of equilibrium is further 
enriched by the possibility of simultaneously applying 
both ac and strong dc fields to make MIRO and HIRO 
superimpose onto each other. A nonlinear response of the 
2DEG to the fields provides a way of mixing the "frequen- 
cies" e ac and e<j c of the two types of magnetooscillations. 



to HIRO. 

As discussed in S ec.lV.D[ a m i nimum of HIRO may 
evolve into a ZdRS ( Hatke et al ., 2010). The nonlinear 
mixing of MIRO and HIRO gives rise to an interplay be- 
tween ZRS associated with MIRO and ZdRS. Namely, as 
£dc increases, ZRS that are present at zero dc field trans- 
form into maxima of r w , as shown in Fig. 136b , for the case 
of ZRS near e ac = 2 + 1/4. In turn, a "strong maximum" 
of MIRO may be converted by the strong dc field into 
a ZdRS. The ZdRS in a microwave-driven 2DEG, origi- 
nating from the maximum o f MIRO near e ac = 2 — 1/4, 
was observed bv lZhang eiail (|2008all2007cr) . In fact, this 
observation appears to be an example of a more general 
(not relying on the presence of microwaves) phenomenon 
which relates ZdRS to the maxima of magnetooscillations 
at zero dc field. The vanishing of r w of a similar nature — 
with ZdRS evolving from the maxima of oscillations — 
has also been observed in the absence of radiation: in 
particular, the emergence of ZdRS from the maxima of 
str ongly developed Sd H oscillations was demonstrated 
by iBvkov et alt (|2007t) . The ZdRS was also found to 



emerge from the maxima of MISO (IBvkov et all . f2010dt 



iGusev et al. l. l201ltlWiedmann et al. 



2011a ) 



2. Frequency-mixing for magnetooscillations 



1. Experiments on MIRO in a strong dc field 

The nonlinear mixing of MIRO and HIRO was demon- 
str ated in exper i ments on ul tra-high mob il ity stru ctures 
bv IZhang etaJl (j2007cT) and lHatke eTaH (|2008bh . Fig- 
ure [36^ shows the effect of a direct current I varied un- 
der continuous microwave irradiation on the differential 
resistivity r w : the MIRO maxima are seen to evolve into 
minima, and vice versa, with increasing /. One of the 
most remarkable features in Fig. 136b is the crossing points 
(marked by arrows) of the curves measured at different 
currents: the resistivity shows little change at e ac = n 
and n + 1/2, where n is an integer. 

More insight into this phenomenon is gained if one 
plots r u as a function of e^c for different fixed e ac , as 
in Fig. [55b . The amplitude of the oscillations of r w in 
Fig. l3"ob is seen to be much larger around e ac = n ± 1/4 
(near which e ac the extrema of MIRO are located) than 
around e ac = n,n + 1/2 (near which e ac MIRO exhibit 
zero response). Figure [35b also reveals a substantial dif- 
ference between integer and half-integer e ac : while at 
e ac = 2,3 the oscillations are suppressed but still clearly 
visible, they virtually disappear at e ac = 5/2,7/2. The 
phase of the oscillations of r w as a function of e^c also 
shows a strong dependence on e ac . Specifically, at integer 
e ac = 2, 3, 7\j behaves very similar to HIRO, with maxima 
occurring at ed c — n, Eq. (| 142|) . In contrast, the oscilla- 
tions of r w at the MIRO maximum for e ac ~ 3 — 1/4 and 
the MIRO minimum for e ac ~ 3 + 1/4 are phase-shifted 
by about a quarter cycle in opposite directions compared 



Analysis o f the experimental data in Fig. I3"5b led 
lHatke et all (|2008bl) to conclude that the local maxima 



/= 69 GHz 2 
T= 1.5 K 



*s 2 




FIG. 36 (a) Differential photoresistivity r u versus magnetic 
field at a radiation frequency of 69 GHz for different levels of 
dc excitation, from 7 = (thick black curve) to I = 20 fiA 
in steps of 1 fiA. Integers show the order of the MIRO peaks. 
The arrows mark the zero-response nodes that remain largely 
immune to dc excitation, (b) Part of the data in Fig. 136a . 
represented as the dependence of r w on edc at different fixed 
e ac ranging from 2 to 3 in steps of 0.05. The traces are offset 
in steps of 0.75 Q. The data were obtained at T = 1.5 K in a 
100 /im wide Hall bar sample with n e ~ 3.7 x 10 11 cm" 2 a nd 



/' ; 



1.2 x 10 7 cm 2 /Vs. Adapted from lHatke et all (|2008bl) . 
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of r u in the e ac -e dc plane occur at the points given by 

e ac + £dc ^ n , e ac - e dc ~ m - 1/2 , (167) 

where n and m are integers [t he first cond it ion w as 
established i n earli er work by IZhang et all ( 2007cl )]. 
lHatke et all (|2008bl ) also argued that Eqs. (fTFfj) point 
to the key role played by backscattering off disorder in 
shaping the oscillations in Fig. l36b (similar to HIRO in 
the limit of e dc ^ !)■ Within a qualitative picture pro- 
posed by lHatke et all ( 2008blh r w has a maximum when 
disorder-induced scattering, accompanied by shifts of the 
cyclotron orbit in real space by 2i? c , is maximized for 
shifts in the direction of the dc field (first condition) 
and minimized for shifts in the opposite direction (second 
condition) . 

As the calculation by iKhodas and Vavilovl (|2008| ) 
shows, the above picture is qualitatively accurate (to first 
order in V and for e ac > e dc ) in describing the mechanism 
of "frequency mixing" of MIRO and HIRO; however, it 
misses an important microwave-induced correction to the 
amplitude of HIRO which arises at the same order in V '. 
Specifically, to order 0(V) and at e dc » 1, the oscillatory 
correction to the resistivity Sr^ = — pu for overlapping 
LLs reads (|Khodas and Vavilovl 120081 ) 

5r u /p D = (16<5 2 /7re dc )(r/T T ) (C + C+ + C_) , 
C = (1 - W) e dc cos 27re dc , 

C±=V (e dc ± e ac ) cos 27r(e dc ± e ac ) . (168) 

The frequency-mixing terms C± are seen to be exactly 
canceled by the linear-in-'P correction to the amplitude 
of HIRO in Co at integer e ac (in which case the oscillatory 
correction Sr^ reduces to the plain HIRO). 

Moreover, Eqs. (jl 6T[) can be obtained from Eqs. (| 168[) 
only if the modified HIRO, described by the term Co, 
are neglected. In general, the location of the extrema 
of r w given by Eqs. (|168[) . and the whole picture of the 
oscillations for that matter, depend in an essential way 
on V . In particular, the most characteristic feature in 
Figs. [3(2 — the vanishing of the oscillations at half-integer 
£ac — can be obtained, within the theoretical framework 
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on which Eqs. (jl68[) are based, at special values of V only. 
Specifically, Eqs. (|168p yield for half- integer e ac : 



Srjpo = (16(5770 (t/tv)(1 - 47>)cos27re dc 



(169) 



The amplitude of the oscillations in Eq. (|169[) is seen to 
be reduced by microwaves. In fact, exten ding the calcu- 
lation of Sr^ for e dc ^> 1 t o arbit rary V ( Khodas et all , 
120101: IKhodas and Vavilovl 120081 ) gives for half-integer 
e ac Eq. (I169P with the factor (1 — 4P) being substituted 
by the Bessel function J {Ay/V) [see Eq. (IT70|) below]. 
The nodes of the oscillations of r w as a function of e dc 
at half-integer e ac are thus only reproduced at zeros 19 
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FIG. 37 Differential resistivity r w versus edc (bottom axis) 
and the direct current I (top axis) at (a) e ac ~ 1/2 and (b) 
e ac ~ 1. The data were obtained at / = 27 GHz and T — 1.5 K 
in a 100 /im wide Hall bar sample with n e ~ 3 . 8 x 10 11 cm -2 
and /i~1.3x 10 7 cm 2 /Vs. From lHatke et all (|2008aT ). 



of Jo (4 VP). Given the accuracy of the reported nodes 
at half-integer e ac , it appears that the data in Fig. 
were likely obtained in the very close vicinity of a zero of 



3. Fractional HIRO 

The nonlinear response of a 2DEG to the dc field was 
also investigated experimentally in the regime of frac- 
tional MIRO, in t he vicinity of the mo st prominent frac- 
tion e ac = 1/2, bv lHatke etail (|20Q8aD . Remarkably, the 
experiment revealed the emergence of fractional HIRO 
at e ac = 1/2, with "frequency-doubling" for the oscil- 
lations of r u as a function of e dc . As seen in Fig. [37j 
the overall behavior of r w at e ac = 1/2 in the interval 
< e dc < 1 closely replicates that at the CR in the in- 
terval < e dc < 2 — with the period of the oscillations in 
e dc at e ac = 1/2 being half that at e ac = 1 (which means 
that the period of the oscillations of r w as a function of 
the direct current is the same). 

A description of the fractiona l HIRO in a microwave- 
excited 2DEG was suggested bv iLeil (|2009l ) in terms of a 
postulated ansatz which assumes an infinitely fast ther- 
mal equilibration in th e frame mo ving with the drift ve- 
locity (see footnote 1). ILeil ( 20091 ) argued that the max- 
ima of r w should obey the condition e dc + ne ac = m with 
integer n and m, and the numerically obtained solution 
to the ansatz equation yielded a picture similar to the one 
in Fig. 1371 however, the dependence of the strength and 
shape of the oscillations on the major parameters, like V 
or the degree of overlap of LLs, remained unclarified. 

On the other hand, the observed frequency-doubling 
for HIRO may be a signal of the inelastic mechanism of 
nonequilibrium magnetooscillations. Indeed, in overlap- 



A reversal of the sign of the oscillations with e^c when V reaches 



a certain threshold, reported bv lZudov et all i2009h . may provide 
a valuable calibration point for the microwave intensity seen by 
the 2DEG, which is difficult to obtain from other experiments. 
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FIG. 38 Differential resistivity r u} (B) (a) for zero direct cur- 
rent I and (b) for I — 64 /iA. The vertical lines mark the po- 
sition of the CR and its second harmonic. Parameters of the 
sample, frequency / = 27 GHz, and temperat ure T = 1.5 K 
are th e same as in Fig. 1371 Adapted from iKhodas et al\ 
(l2010h . 



pingLLs, a controlled calculation bv lKhodas and Vavilovl 
(2JD08>) at order S 2 and at first order in V in the collision 
integral showed that, while the displacement mechanism 
is dominant in the limit of e^c 3> 1, the inelastic mech- 
anism for Tj n /r > 1 is equally important at e^c ~ 1 (see 
also Sec. V.C) and produces strong subharmonics in the 
dependence of r u on edc — in contrast to the displacement 
mechanism at this order in 8. Importantly, the strength 
of the subharmonics in the p resence of radiation depend s 
on V and oscillates with e ac (K hodas and Vavilovl . 120081) . 
In particular, for overlapping LLs, the subharmonic with 
a period in edc equal to 1/2 is as strong at e ac = 1/2, 
edc ~ 1) and V ~ 1 as the harmonic with period 1. 
In the case of separated LLs, the half-integer MIRO 
and half-integer HIRO are also likely to strengthen each 
other in analogy with the theory of two-photon absorp- 
tion (Sec. VI. A). Further theoretical work is necessary to 
study the exact form of the interplay of the two types of 
fractional oscillations. 



4. AC field-periodic oscillations 



In Sec. IVI.B.3l the focus was on the regime of a mod- 
erate microwave power {V < 1). The nonlinear mixing 
of HIRO and MIRO in the case of high-power radiation 
(V ^> 1) results in the emergence of oscillations of r u 
as a function of V (in addition t o the oscillations i n edc 
and e ac ). These were observed ( Khodas et <zZj . l2O10h as a 
series of oscillations grouped around the CR and decay- 
ing away from it, with much less pronounced but similar 



behavior also detected around e ac = 2. In contrast to the 
integer MIRO, which have one maximum and one mini- 
mum per unit interval of e ac (Fig. 138a ) , or the fractional 
MIRO with their oscillatory features around rational val- 
ues of e ac , these oscillations exhibit multiple maxima and 
minima around the CR (and its second harmonic) with 
a period which depends on ?, as shown in Fig. 



A theoretical descriptio n by IKhodas et al. I (|2010h : 
IKhodas and Vavilovl (|2008h of the reported oscillations 
in V gives for the oscillatory correction Sr^ at order S 2 
for e dc > max{l, VV}: 



edc cos(27redc) Jo (4 VP sin 7re a 



5r u _ 16S 2 t 
— 2e ac sin(27redc) cos(7re ac )v / 7- > J\ ( 4 VP sin7re 



(170) 



which generalizes Eqs. (|168|) to the case of large P, with 
the Bessel functions Jo,i oscillating with a period which 
depends on V and e ac . Oscillations in V of this form are 
a characteristic signature of multiphoton transitions be- 
tween stationary states of a periodically (in time) driven 
electron system. Namely, the absorption rate for TV- 
photon transitions between these states is proportional 20 
to Jj^[2-\/Psin(0/2)], where 6 is the angle between the 
initial and final direct ions of the ve l ocity (in the qu asi- 
classical formulation) (|Khodas et aZl l20l'd : iLel l2009h . In 
the limit edc ^ 1 , where the main c ontribution to HIRO 
comes from 9 ~ 7r, summing over TV (jKhodas et al l l2010D 
reproduces Eq. (|170| . The characteristic number of pho- 
tons N in the transitions that yield Eq. (|170[) for V ^> 1 
is 21 2\fV . This conclusion is further corroborated by the 
observation that r u from Eq. (|170l) in the limit V ^> 1 
oscillates, as a function of e ac in the vicinity of inte- 
ger e ar , with a perio d \jl\fV. The experimental data 
(Kho das et al ., 2010). shown in Fig. l39k . demonstrate 
that the oscillations of r u with e ac become more rapid 
with increasing microwave power, in close agreement with 
the result of Eq. (|170j) . also shown in Fig. [39b . 



20 A similar structure arises in a variety of closely related prob- 
lems, e.g., in study ing the multiphoton CR in three dimen- 
sions jSeelvl, I 1974) or t he dynamical Franz-Keldysh effect 
(Joh nson and Jauhclll99Sft . 

21 More specifically, for 1/P 1 / 2 <C |eac|, where e a c = £ac — n and 
n is the nearest integer, the sum over N, which is of the type 
E;v^Ar( 2 v / 7>)(edc + A^tac) cos[27r(e dc + 2Ve ao )] and contains the 
factor oscillating in TV with a period l/|e ac |, is determined by a 
close vicinity of the sharp peak in J^(2\/V) at N = 1\/V whose 
width is of order P 1 / 6 . The peak falls off exponentially on the 
side of larger N and as (2\/P — TV) -1 / 2 on the other side. For 
l/P 1 / 2 <C I eac I l/P 1 / 6 , the main contribution to Sr^ comes 
from a narrow interval of positive 2vP — N whose width is of 
order l/|e ac |. 
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FIG. 39 Measured (a) and calculated (b) differential resis- 
tivity versus e ac for / = 54 /iA and different microwave 
intensities. The traces are labeled according to (a) the at- 
tenuation levels and (b) the value of V ■ Parameters of the 
sample, frequency / = 27 GHz, and te mperature T = 1.5 K 
are the same as in Figs.[37]and[38] From lKhodas et al.l (|2010T ). 



VII. OUTLOOK 

A. Perspectives 

Here, we briefly outline several, in our view, particu- 
larly promising directions of ongoing and future research 
in the field of nonequilibrium magnetotransport phenom- 
ena in semiconductor structures. 



a very strong negative MR (also in the absence of mi- 
crowaves) attributed to the quasiclassical memory effects, 
see Sec. lII.B^2l which indicates a possible connection be- 
tween these phenomena. As a matter of fact, the memory 
effects that are responsible for the QCMR may also pro- 
duce spikes in the photoconductivity at the CR harmon- 
ics, see Sees. III. B.2HII.B.31 and lIII.B~5l However, it is un- 
clear why only even harmonics (and, most prominently, 
the one at n = 2) are enhanced. Further theoretical work 
is thus needed. 
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FIG. 40 Magnetoresistivity measured at f — 121.45 GHz and 
T = 0.3 K in a Hall bar sample with n e ~ 2.9 x 10 11 cm" 2 
and fi ~ 3.0 x 10 7 cm 2 /Vs. The dark resistivity, normalized 
to its value at B = 0, is shown in the inset. Adapted from 
iDai et al\ (j2010h . 



1. Domain structure of ZRS and ZdRS 

The present status of research in this area has been 
reviewed in Sec. lIVl The physics of the broken-symmetry 
states with spontaneously formed domains remains one of 
the central research directions in the field. Particularly 
important open questions include the character of the 
transition into ZRS, the corresponding critical behavior 
(in the case when the transition is continuous) beyond the 
mean-field picture, the residual resistance of a finite-size 
system, the effect of thermal and nonequilibrium noise 
on ZRS, as well as dynamics of the domain structure. 



3. Nonequilibrium phenomena in separated Landau levels 

Nonequilibrium phenomena in the regime of fully sepa- 
rated LLs deserve further study. On the theoretical side, 
the theory for mixed disorder, which is the type of disor- 
der characteristic of most of the experimentally relevant 
samples (see Sees. IIIII and PVT) . needs to be extended into 
this regime. On the experimental side, this corresponds 
to studying the nonequilibrium phenomena in stronger 
magnetic fields. It is further interesting to see how the 
phenomena discussed in this review evolve when entering 
the QH regime, see also Sec. IVII.A.71 



2. Giant photoresponse at the second harmonic of the 
cyclotron resonance 

Recent experiments on ultra-high mobility structures 
reported a narrow spike a t u c± 2oj r . superposed on the 
n = 2 maximum of MIRO (|Dai et aq[2010l ; lHatke et all 
l2011bl ). The height of the spike was several times (up 
to an order of magnit ude) large r than the height of 
other MIRO maxima. IDai et al. also observed 

an enhancement of higher-order even MIRO maxima 
(n = 4,6,8), see Fig. 001 It is important to note that 
these effects were observed in samples which also show 



4. Nonequilibrium transport at low temperature 

Most of the experiments and theories presented in this 
review were focussing on the temperature regime where 
SdH oscillations are suppressed. The interplay of the 
nonequilibrium physics and SdH oscillations at lower T 
is also of interest. A remarkable pheno menon of this 
kind t hat was studied both experimentally (iBvkoy et al 



2007 ; Kalmanovitz et Ml 120081: IStudenikin etal 



2011 



Zhan g et all l2009l l2007al) and theoretically (|Dmitriev 



20111 ) is the dc field-induce d inversion of t he pha se of SdH 
oscillations. In particular, iBvkov et all ( 2007|) reported 
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the formation of ZdRS that arise, as the strength of the 
dc field increases, from the maxima of SdH oscillations, 
see Sec. IV. Dl A similar inversion of the SdH oscillat ions 
by the microwave field was observed bv lManil (12004ft . 



5. Spin-related phenomena 

In the regime of strongly developed SdH oscillations, 
the s pin degree of f reedo m starts to play an important 
role. iRomero et all ( 20081 ) explored the effect of a paral- 
lel magnetic field B\\ on ZdRS in this regime (Sec. IIV.DP 
and found that the ZdRS are destroyed for sufficiently 
stron g _B|| . The result was explained (IRomero et all , 
l2008f ) within the kinetic equation formalism described 
in Sees. IIIIIVI in terms of a Zeeman splitting-induced 
modification of the osc illatory component of the DOS. 
IStudenikin et all ( 201lh further extended the study of 
the effect of Zeeman splitting on nonlinear magnetotrans- 
port to the regime of strong magnetic fields normal to the 
2DEG. Still, spin-related phenomena in nonequilibrium 
magnetotransport remain largely unexplored. This re- 
search direction overlaps with the investigations of the 
regime of separated LLs (Sec. IVII.A.3[) and of the QH 
regime fSec. IVII.A~7)) . 



6. Polarization dependence of MIRO 



As discussed in Sec. lIII.B.l.bl the relative orientation 
of the dc field and the linearly polarized microwave 
field affects considerably the MIRO amplitude within 
the framework of the displacement mechanism, while 
the inelastic contribution is insensitive t o the d irectio n 
of the linear po l arizatio n. iMani et all (|2002l . l2004ct) : 
IWiedmann et all ( 2011al ) reported equal MIRO signal 
strength for the microwave field polarized along and per- 
pendicular to the current direction, which is consistent 
with the dominance of the inelastic contribut ion. On the 
other hand, very recent experimental study ()Mani et all , 
l201ll ) reported a pronounced difference of the MIRO am- 
plitude for two directions of the linear polarization, in- 
dicating the relevance of the displacement contribution. 
Systematic experimental study of the dependence of this 
effect on the temperature and microwave power would 
permit a quantitative comparison with theoretical pre- 
di ctions (see Sec.|III.B.l.g[ ) . 

ISmet et all (2005) used an optical setup in order to 
study MIRO induced by circularly polarized microwaves. 
According to the theories presented in Sec. IIIIl the ampli- 
tude of MIRO depends on the direction of circular polar- 
ization through the factor l/[(w ± u> c ) 2 + 1/t 2 ], where 
the minus (plus) sign corr esponds to a ctive (inactive) 
polarization. Surprisingly, ISmet et~al. I (|2005l ) reported 
very similar results for two directions of circular polar- 
ization, with a difference being much smaller than ex- 
pected from this factor. To understand a possible reason 
for the failure to see a big effect after changing the direc- 



tion of circular polarization of the incoming radiation, it 
is important to take into account that the electric field in 
Eqs. dHJ), ([77]). and fl52J| is the total electric field acting 
on electrons in the sample. 

First of all, the total field is different from the 
source field because of screening by the 2DEG, see 
Sees. III.C.31lHI.B.l.cl This leads to a strong broadening 
of the CR: the above factors are modified to l/[(w±a; c ) 2 + 
7 2 ], where 7 = 2ne 2 n e /mcn r [Eq. (|57|) ] for 7 ^> 1/t. For 
a typical concentration n e — 3 • 10 11 cm" 2 and the ef- 
fective refractive index n r = (V12.8 + l)/2 ~ 2.3 for the 
case of a GaAs-based structure one finds 7/27r ~ 17 GHz, 
which is comparable to characteristic microwave frequen- 
cies. While this explains why no sharp resonance features 
is experimentally observed near uj — uj c , this is not suffi- 
cient to explain why the difference in the amplitude of the 
MIRO harmonics with n = 2, 3, 4 for two d irections of the 
circular polarization in ISmet et all ( 20051 ) was within a 
few percent. Indeed, since 7 is several times smaller than 
the frequencies uj/2tt — 100~250GHz used in this experi- 
ment, one would expect a much larger difference (a factor 
of about 9 for n — 2 and of about 4 for n = 3, etc.). This 
discrepancy can be most likely attributed to the fact that 
the metallic electrodes and other elements of the sample 
setup may strongly affect the polarization of the e lectric 
field. In particular, iMikhailov and Savostianoval (2006) 
showed, for a simple geometry, that metallic contacts 
may strongly reduce the degree of circular polarization 
of microwaves [th e actual geometry of the experiment 
( Smet et al .1. 120051) is much more complex, and a quanti- 
tative evaluation of the polarization of the actual electric 
field would require a numerical modeling of the whole 
setup.] 

In fact, a close inspecti on of the tr a nsmis sion data 
shown in Figs. 1 and 2 of ISmet et all ( 20051 ) provides 
additional support to this explanation. While the data 
do show, as expected, a CR dip around uj c — ui for the ac- 
tive polarization, the — l/(u)~ uj c ) 2 tail of the Lorentzian 
for the active polarization is not observed when \uj c \ be- 
comes smaller than uj by more than 7. Instead, the data 
for both polarizations of the incoming wave become prac- 
tically identical, with the transmission amplitude for the 
active polarization being even slightly higher. This sug- 
gests that, in this range of magnetic fields, the polariza- 
tion of the incoming radiation was essentially immaterial 
for the actual polarization of the electric field. Since this 
is exactly the range of magnetic fields where MIRO are 
measured, this observation is consistent with the fact that 
the MIRO signal was almost identical for both directions 
of the circular polarization of the incoming wave. 

It appears that the presence of nearby metallic con- 
tacts and, more generally, the effects of a finite sample 
geometry may strongly affect the actual polarization of 
microwaves, thus seriously complicating the study of the 
polarization dependence of MIRO. Devising an experi- 
mental setup that would be free from these drawbacks 
remains a challenge for future work. 
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7. Nonequilibrium magnetotransport in the integer and 
fractional quantum Hall regime 

It is natural to ask whether the phenomena described 
in this review have their counterparts in the QH regime, 
both on a QH plateau and near a QH transition. In 
this respect, it is worth recalling another nonequilibrium 
phenomenon — the breakdown of the QH effect at a high 
current d e nsity- — that has been extensively studied, see 
iNachtweil (jl999t) for a review. While this is an essen- 
tially different phenomenon, certain aspects bear simi- 
larity with the MIRO/HIRO and, especially, ZRS/ZdRS 
physics (the nonequilibrium phase transition between 
states with zero and finite resistivity in a transverse mag- 
netic field; the inhomogeneous spatial distribution of the 
electric field and current; the importance of relaxation 
processes; etc.). It remains to be seen whether there is a 
deeper connection between the two classes of phenomena. 

It would be further highly interesting to observe phe- 
nomena analogous to those reviewed here near half- 
filling of the lowest LL, where transport is associated 
with composite fermions (CFs) (jHeinonenl . 1 19981 : iJainl . 
120071 ). A possib i lity to observe MIRO of CFs was dis- 
cussed bv iParkl (|2004l) . There is, however, a serious 



difficulty related to the fact that CFs are much more 
strongly scattered by disorder than electrons. The dom- 
inant type of disorder for CFs is an effective random 
magnetic field produced by the random scalar poten- 
tial through the flux attachment ( Halperin et [l993l) . 
The damping of quantum magnetooscillations induced 
by the random magnetic field is much stronger than 
the one (oc e -^/ ^c T q) characte r istic of the regime of low 
magnetic fie l ds (lAronov et al. |, 119951: lEvers et all Il999t 
iMirlin et all Il996l 19981 ). This appears to be an essen- 
tial obstacle for the observation of pronounced MIRO or 
HIRO in this regime. On the other hand, while the DOS 
oscillations around v = 1/2 are very strongly damped, 
one does observe conductivity oscillations which are of a 
different origin: they are induced by quantum localiza- 
tion of CFs, leading eventually to the fr actional QH effect 
(|Evers et all Il999l: IMirlin et all Il998l ). An intriguing 



question, deserving theoretical and experimental investi- 
gation, is whether these oscillations would be sufficient to 
induce MIRO or HIRO under nonequilibrium conditions. 



B. Nonequilibrium magnetotransport phenomena in other 
2D systems 

1. Electrons on liquid helium 

Electrons on a surface of liquid helium represent a re- 
markable 2D electronic system with parameters strongly 
different fr om thos e of 2 D EGs in s e micon ductor struc- 
tures, see lAndreil (jl997T ); iGrimed (|l978l ) for reviews. 
First, because of a very low density, correlations play 
an essential role: in the zero-temperature limit this sys- 
tem would be a Wigner crystal. Second, the electron 
gas on helium becomes nondegenerate already for very 



low temperatures. Finally, relaxation times at Kelvin 
temperatures are several orders of magnitude longer 
than in semiconductor 2DEGs. Quasiclassical magne- 
totransport in this system has been studied both exper- 
i mentally and theoretically during the last two decades 
(iDvkman et a/J. 1 19971 111)93 iDvk man and Prvadkoil200^ ; 



iLea and DykmanT " 1998 ). The focus of a more re- 
cent activity was on the photoresponse of the electron 
fluid. A characteristic feature of the photoresponse is 
resonant absorption at a frequency of order 100 GHz 
(jKonstantinov et all 120091 120071 l2008allbh . correspond- 
ing to the transition from the lowest to the first excited 
subband (i.e., the first excited state for the motion in 
the direction transverse to the surface). It is worth men- 
tioning that the conductivity of electrons on the helium 
surface is measured in a contactless way, by capacitively 
coupling the 2D electron system to electrodes and mea- 
suring the complex conductance (admittance) of the cir- 
cuit. 

In very r ecent experiments (jKonstantinov and Konol . 
120091 . 12010h . the resonant intersubband excitation mea- 
surements were performed in a transverse magnetic field. 
Remarkably, the system showed 1/B-periodic magne- 
toconductance oscillations and zero-conductance states 
(the experiments were performed in the Corbino-disc ge- 
ometry), very similar to those found in the semicon- 
ductor 2DEG (Sees. [TlTl[lVj). The similarity includes 
the same periodicity and phase, as well as a strong^ 
tempe rature dependence. Further, iKonstantinov et al " . 



(j201_]|) reported spontaneous currents (in the absence 
of driving voltage) in the regime of zero-conductance 
states. An explanation of the s e oscillations put for- 
ward b y lKonstantinov and Konol ( 2009. 2 lot ) : [Monarkhal 
(|201lh is analogous to the MIRO theory (SecHIIl) and 
combines features of both the inelastic and displacement 
mechanisms. Specifically, it is proposed that the resonant 
microwave absorption leads to a strongly nonequilibrium 
occupation of the first excited subband, which is followed 
by a quasielastic transition from the lowest LL of the first 
excited subband to one of higher LLs of the lowest sub- 
band. 

The above observations appear truly remarkable as 
they show that the nonequilibrium phenomena in high 
LLs, discussed in the review, are also relevant to strongly 
interacting and nondegenerate 2D electron systems. On 
the other hand, a systematic theory of the nonequilibrium 
magnetotransport phenomena in the electron liquids on 
helium remains to be developed. It particular, one has to 
clarify the role of resonant absorption to the excited sub- 
band, the relative weight of the contributions of different 
mechanisms, relevant relaxation processes, etc. 



2. Graphene 

The ex perimental brea k throu gh in fabrication of 
graphene ( Novoselov et all |2004) has brought a new 
2D system with fascinating properties into condensed- 
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matter physics ( Geim and Novoselovl 120071 ). The hall- 
mark of graphene is the massless Dirac-fermion char- 
acter of charge carriers. It manifest itself, in particu- 
lar, in the unconventional QH effect bot h in graphene 
dNovoselov et all 120051: IZhang et aZ.l.l2005l) and graphene 



bilay ers (|McCann and Fal M, [2006; iNovoselov et all 
2006) . The Dirac spectrum implies also a linearly vanish- 
ing DOS at t he ne utrality point, i.e., the effective mass 
vanishes as vT^eT as a function of the electron concen- 
tration or, equivalently, of the gate voltage. As a result, 
the LLs are not equidistant, with the distance between 
them (cyclotron frequency) being strongly enhanced near 
the Dirac point. One of the manifestations of this is 
the o bservation of th e QH effect up to room tempera- 
ture ( Novoselov et all 12007). For a r ecent review of the 
electr onic properties of graphene see ( Castro Neto et all , 
2009). One of the most active directions in the graphene 
research is devoted to the optoelectronic properties of 
graphene, including a variety of prospective applications 
(photodetectors, l i ght emitters, solar cells, lasers, etc) 
(|Bonaccorso et q/.l . l2010t iMueller et all , l2010h . 

Spectroscopy studies of graphene in a transverse mag- 
netic field confi r med th e Dirac-fermion LL spec trum 
(jSadowski et all 120061 ). iNeugebauer et all (|2009h pre- 
sented a spectroscopic manifestation of the exception- 
ally high purity of graphene samples: the LL quantiza- 
tion was observed at B as low as 1 mT and T up to 
50 K. For a review o f LL s pectroscopy of graphene see 
(|Orlita and Potemskil |20Tfj) . 

The properties of graphene make it a very interest- 
ing and promising system for studying nonequilibrium 
physics in quantizing magnetic fields. Several differences 
compared to the semiconductor systems may be foreseen. 
First, according to the above, one can expect that in 
high-quality graphene samples with the chemical poten- 
tial located near the Dirac point the quantum phenom- 
ena will become pronounced at relatively high T and 
relatively low B. Second, the nonequidistant character 
of LLs may lead to a modification of the functional de- 
pendence of magnetooscillations on T and the strength 
of disorder. In fact , a first experim ent in this direction 
appeared recently: iTan et al. I (1201 ll) studied SdH oscil- 
lations in graphene under dc bias and observed a phase 
inversion analogous to that found in semiconductor struc- 
tures, see Sec. IVII.A.4I 
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